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1.—The spectrum of a nova, before and during maximum brightness, 
is not unlike that of a supergiant star of Class B, A or F, with the ab- 
sorption lines greatly shifted toward the violet. After maximum the 
spectrum changes rapidly with the development of wide, complex bright 
lines, undisplaced in the mean from their normal positions; complex struc- 
ture, corresponding to several different velocity shifts, also appears in the 
absorption lines, which become progressively fainter as the bright-line 
spectrum increases in intensity. 

It has long been realized that the violet shift of the absorption lines is a 
result of a rise of the absorbing material, and this upheaval of the star’s 
surface helps to account also for the brightening of the nova, since it is 
evident that the rapid increase in brilliance cannot be attributed solely 
to increase in surface brightness; the area of the radiating surface must be 
increasing as well. It has been suggested that a nova outburst begins with 
the expansion of the star as a whole, and that near maximum brightness 
the star ejects a shell of nebulous matter, and then contracts and grows 
fainter. Though this model interprets many features of the behavior of 
novae, it fails in several important points, and the present note attempts to 
outline another model that may be more satisfactory. 

2.—We suppose that the expansion is confined to the photospheric layers, 
which are blown away from the stellar surface at the time of the first out- 
burst. The star then becomes a sort of planetary nebula, in structure 
if not in physical conditions. In the early stages the shell should still be 
dense enough to radiate like a stellar photosphere, receiving diluted high- 
temperature radiation on its inner surface, and reradiating the transmitted 
radiation undiluted at a lower temperature. The observed progression 
of the nova spectrum from an absorption spectrum to a bright-line spec- 
trum, which soon develops the nebular lines and ultimately shows Wolf- 
Rayet characteristics, is in general harmony with this model; that the 
theory is also substantiated in considerable detail will be shown later in the 
present note. 
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3.—The profile of an absorption line in an expanding photosphere of 
constant (high) velocity depends upon two factors: the variation of 
intensity of the line over the stellar disc, and the degree of darkening of the 
photosphere toward the limb in the wave-lengths that border the line. In 
the case of a nova, where the velocity shifts are of the order of 1000 
km./sec., it is unnecessary to know the form of the original line profile; 
we assume that the profile is rectangular and that the intensity within the 
line is constant. In figure 1 are shown schematically the predicted profiles 
corresponding to four assumptions: (a) absorption line of constant inten- 
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FIGURE 1 


Schematic profiles of lines in the spectrum Schematic profiles given by a transparent 


of a nova expanding with constant high 
velocity. Bright lines are blackened. A 
dotted line indicates the undisplaced wave- 
length; violet is to the left. The four cases 
represented are: (a) uniform disc; (0) disc 
darkened at the limb; (c) lines dark at the 
center, bright at the limb; (d) emission 
from transparent expanding shell. 


expanding: shell, ejected directionally as 
shown in the upper right corner. Four pre- 
dicted profiles, seen as the star is viewed 
from the direction of the four arrows, are 
shown. Ejection is supposed to be more 
plentiful at the center of the jet. The dia- 
gram is two-dimensional; the fact that the 
actual case is three-dimensional does not 





change the effects essentially. 


sity on expanding disc of uniform brightness; (6) absorption line of constant 
intensity on expanding disc darkened toward the limb; (c) ‘“‘chromospheric’”’ 
emission present (line dark at the center of the disc, bright toward the 
limb); (d) shell effectively transparent, except perhaps to the radiation 
of the nuclear star. The last case, discussed by Beals,! gives a rectangular 
emission profile. 

The foregoing assumptions are too simple. Unsymmetrical ejection of 
material would obviously affect the observed profile. In figure 2 is shown 
the predicted profile for a star with an effectively transparent radiation 
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shell, ejecting matter symmetrically with respect to two hemispheres, and 
viewed from several angles. The superposition of successive outbursts 
would also complicate the profile. Further, the predictions of figures 1 and 
2 were made on the basis of a single constant velocity, but, as one of the 
writers has pointed out,? any symmetrical profile can be predicted by 
appropriate choice of velocity distribution, for an atmosphere expanding 
isotropically. 

4.—A detailed study, now being made by the writers, of the available 
material of the spectra of bright novae at Harvard and elsewhere, will 
ultimately provide data for comparing observed and predicted profiles. 
Cases (a) and (5) have not at present been discriminated, but a study of 
good premaximum spectra should provide the necessary data. Case (c), 
“‘chromospheric’”’ emission, is shown by the premaximum spectra of Nova 
Persei, Nova Aquilae, Nova Cygni and Nova Pictoris, and is doubtless 
a normal feature of the early spectra of all novae. Case (d), the rectangular 
profile, is not shown by any novae that we have at present examined. 
Some of the bright lines show a rounded form, and others, while sharply 
bounded at the red and violet, show a central dip; the two types coexist 
in the spectrum of Nova Aquilae, where the hydrogen lines show a saddle- 
shaped, and the nitrogen lines a rounded profile. Both of these types 
of profile could be interpreted equally well in terms of appropriate velocity 
distribution, or of directional ejection. The structure of individual bright 
lines in the spectrum of Nova Aquilae’ clearly points to directional ejection 
in jets, with respect to a definite plane which may be defined by the star’s 
equator. 

5.—The interaction of the greatly widened profiles of neighboring lines 
is of much interest. Emission profiles overlap in a purely additive way, 
producing a complex structural appearance. The relation between emis- 
sion and absorption lines seems to indicate that absorption takes place at 
the lower level of the stellar surface, while the emission lines are produced 
above them in the tenuous envelope. This conclusion is drawn from the 
fact that in the case of a close group of lines, all the emission components 
are present but only the violet absorption component. This observation 
and this interpretation also apply to the Wolf-Rayet stars. 

6.—Most of the interpretation of the spectra of novae that has hitherto 
been made has hinged on the observations of radial velocity, and hence has 
dealt primarily with the state of motion of the outer shell. Extensive 
analysis by Wright,‘ by Spencer Jones® and by Stratton® has led to the 
possibility of drawing conclusions also as to the composition and constitu- 
tion of the atmosphere. 

7.—Recent investigations of the spectra of the W (Wolf-Rayet) stars have 
resulted in the almost complete identification in them of the lines of carbon, 
nitrogen and oxygen up to and including the three-electron stage.” The 
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temperatures of these stars, estimated by the method of Zanstra,’ are of 
the order of 50,000° to 100,000°. The W stars, as has been shown by 
Menzel’ and by Beals,” are losing their atmospheres at a slow but steady 
rate with high velocity, so that they are surrounding themselves with 
nebulous envelopes, which are however of too high a density to permit the 
“forbidden’”’ nebular lines to appear. The novae present an obvious 
parallel, since they too are losing their atmospheres with high velocity, 
and since they ultimately show a characteristic W spectrum. 

8.—When a developing nova is regarded from this standpoint the 
successive spectral changes are seen to form a logical sequence. The 
absorption stage is followed by the successive appearance of normal bright 
lines and forbidden bright lines. The absorption spectrum changes, as 
the star passes maximum, in a way that indicates a rapid increase of excita- 
tion; for example, the appearance in Nova Aquilae" of the lines of N V 
indicates that the star has risen from perhaps 10,000° to an ionization 
temperature of about 100,000° in the course of a few months. The normal 
bright lines appear in order of excitation energy, the bright-line spectrum 
lagging appreciably behind the absorption spectrum of similar excitation. 
The forbidden bright lines appear also in order of excitation energy. The 
observed forbidden transitions usually occur for atoms whose ground states 
have the configurations p’, p* and p*. These configurations all produce 
two metastable terms, with three possible transitions between each other 
and the ground state, giving rise to what we have termed nebular, auro- 
ral and transauroral lines. For any one atom the auroral lines appear first, 
and are followed after a few days (Nova Aquilae) or a few weeks (Nova 
Pictoris) by the nebular lines. So regular is the sequence of appearance 
of lines of a given type that the date at which a line first appears furnishes 
an effective means of line identification in the nova spectrum. Reference 
should be made to an earlier discussion of forbidden lines in astrophysical 
sources. !? 

The nebular lines gradually increase in strength, and, in some novae, 
eventually surpass the auroral lines. They are almost always more promi- 
nent in the spectra of nebulae than the auroral lines. Bowen" has sug- 
gested that the difference of intensity of the nebular and auroral lines of 
O I in the spectra of the terrestrial aurora and the planetary nebulae may 
be a result of different relative lives for the metastable states. 

9.—The analysis of the spectrum of the developing nova along these 
lines promises to be an extensive investigation. As a general illustra- 
tion of the data, tables 1, 2 and 3 present the observations of prominent 
lines for three typical novae. The tables list the atoms whose spectra 
appear on the various dates, with the excitation potentials of their strongest 
lines. Since the material of the nova is probably not in thermodynamic 
equilibrium, the significant excitation potential is at present uncertain. 
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For convenience we have listed the cumulative potentials (in volts) counted 
from the ground state of the neutral atom up to the state of the atom from 
which the line in question arises. 

Any astrophysical study that employs marginal appearances is open to 
question, since the dates of the appearances are, no doubt, influenced by 
factors other than the excitation potential, for example, by abundance. 
The dates upon which the various lines were at maximum appear to be 
much more valuable indices of the excitation energies, a conclusion sup- 
ported by the meager evidence now available. The authors plan to in- 
vestigate this point further. 





TABLE 1 
SPECTRAL DEVELOPMENT OF NOVA PERSEI 
EMISSION 
DATE ABSORPTION - 
1901 PERMITTED AURORA! NEBUI.A? TRANSAURORAI. 
Feb. 22 MgI_ 2.7 
Ca FP.) Gal 
H 10.2 H 12.0 


Mg II 16.4 MgII 19.2 
Sis H..17.9) St T--20::9 


He I 20.9 
Al. TL 21.0 
CS: Tk Se:2 
N II 32.9 
OQ: Tt 29.2 
Al III 40.3 
He II 72.6 
Feb. 23 oe. TE - 626 
Th TE 1056 
Feb. 24 Mg I 5.1 
Ca II 9.2 
Fe II 13.2 
Feb. 25 He tf 2.00 1 7 ee. a | 2.4 8 I 2.3 
See | 2.7 
Feb. 28 ?N II 18.5 ?0 II 16.9 
Mar. 2 N_ II* S$ HE 18:3 
Mar. 19 He II 75.3 
NN Tt 778 
Apr. 12 O III 53.8 
Apr. 28t Ne IV 132 O - Tif 51.0 


Ne III 65.4 
* Reappearance. 
{1 No observations since April 13. These lines may have appeared earlier than 
April 28. 


10.—The excitation analysis of the bright-line spectra of novae furnishes 
a basis of comparison for novae that developed at very different rates. 
There are of course differences of detail, but all novae at the same excita- 
tion stage show approximately similar spectra. 
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TABLE 2 
SPECTRAL DEVELOPMENT OF Nova AQUILAE 
EMISSION 
DATE ABSORPTION — “ 
1918 PERMITTED AURORAL NEBULAR TRANSAURORAL 

June 8 H 10.2 H 12.0 

Mg II 16.4 

< ie > Ge ba 
June 9 Fe I 0.0 

at. 49 

Fe II 10.7 
June 10 Fe II 13.2 
June 12 He I 23.0 
June 14 N_ III 74.3 
June 15 (before) Ot 4.2 Sok 2.3 
June 15 He II 75.3 
June 16 N_ III 77.0 
June 18 N II 18.5 N II 16.4 

O II 16.9 

June 24 N_ IV 141.5 
July 3 N_ IV 145 Ne III 65.4 
July 5 O III 53.8 
July 8 © Tit 61.0 
July 12 Ne IV 132 


Aug. N V_ 225 


The most obvious difference between the bright-line spectra of different 
novae is found in the behavior of the lines of ionized iron. In only two 
novae (Nova Pictoris'* and Nova Aurigae'®) have the forbidden Fe II lines 
appeared during the outburst, though the celebrated ex-nova 7 Carinae 
has developed them with great prominence'® since 1889, and shows them 
still. Evidently the novae that develop the most slowly exhibit the 
forbidden Fe II spectrum, probably because only for such objects is the 
requisite low density attained before the temperature rises so high that all 
the Fe II is further ionized. The slow fall in brightness and the slow 
spectral changes of Nova Pictoris and Nova Aurigae are almost certainly 
associated with a slower rise of temperature than takes place for most novae. 
Possibly some of the fast developing novae may, in a similar way, fail to 
exhibit other forbidden lines of low excitation; further observations are 
required on this point. 

11.—In the final stages of the nova the forbidden lines usually disappear. 
This may be attributed in part to the general increase in ionization and to 
the decreased brightness of the star, but the persistence of the lines of H 
and He II shows that an abundant supply of energy is still available. 
Bowen" and Zanstra™ look with favor on the hypothesis that the forbidden 
lines in nebulae arise from collisional excitation rather than from recom- 
bination. In the late stages of a nova collisions would then be too infre- 
quent, owing to the low density, to excite an appreciable number of atoms. 
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TABLE 3 
SPECTRAL DEVELOPMENT OF Nova PICTOoRIS 
EMISSION 
DATE ABSORPTION - “ 
1925 PERMITTED AURORAL NEBULAR TRANSAURORAL 
May 25 Al I 0.0 
Cat 0.0 
Fe I 0.0 
Sis Ha | 0.8 
Mg I 2.7 
= a2 866 
Ce ff 6.4 
ae EF. 330 
oe WB 
H 10.2 
Fe II 10.6 
Mg II 16.4 
June 10 Caot 2.9 
Mn I 3.1 oa S 2.7 
Fe I 4.5 
Ca Te Be 
Se Tk SS 
Ti. 10.6 
H 12.0 
Fe II 13.2 
June 12 Nac ke 2 NM T*: 24° 8" 2.3 
See ue. eg 
Zr II 10.5 
Ce. Tt. 33:5 
He I* 23.0 
June 25 ?Si I 1.9 
Aug. 17 Fe it 10.7 Fe II 11.0 
Aug. 24 oe” S68 
Dec. 21 O III 53.8 
1926 
Jen, 7 He If 23.0 ©. EE 6r:0 
Ne III 65.4 
Mar. 15 He IIt 75.3 Ne IV 1382 
N IIIf 77.0 
Sept. 7 F III§ 56 


* Earliest record. 


t Reappearance. 
t These lines probably appeared much earlier; this is the first date on which they are 


sharp and well resolved. 
§ This is the date given by H. Spencer Jones. The Harvard plates seem to be 
systematically earlier than the Johannesburg plates and this line may have appeared 


before September 7. 


The mechanism of excitation of the bright lines will be discussed in a later 


paper. 
The bright-line spectra of novae can be interpreted in a general way, in 


terms chiefly of rapidly rising temperature, rapidly falling pressure and 
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rapidly increasing dilution of radiation. But the problem is made some- 
what more complex by the probable existence of differences of chemical 
composition. In comparing various novae, the irregularity in the distribu- 
tion of carbon, noted for the W stars," also appears. In the development 
of some novae the permitted lines of several’successive stages of the carbon 
atom have been noted; from others they are absent. Unfortunately there 
are no accessible forbidden lines of these stages of the carbon atom. The 
apparent high abundance of neon and fluorine in the atmospheres of the 
novae leads one to suppose that these atoms are cosmically very abundant; 
it should be noted that neon is conspicuous in the nebular spectra, while 
fluorine isnot. A search for the stellar absorption lines of these elements is 
to be desired; and a search for irregularities in cosmic abundance is now in 


progress. 
We wish to express our appreciation to Dr. J. H. Moore and Dr. W. H. 


Wright of the Lick Observatory for valuable suggestions made while this 
paper was in manuscript. 
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CHROMOSOME NUMBER AND MORPHOLOGY IN NICOTIANA 
VI. CHROMOSOME NUMBERS OF FORTY SPECIES' 


By THomMAS HARPER GOODSPEED 
DEPARTMENT OF BOTANY, UNIVERSITY OF CALIFORNIA 


Communicated May 29, 1933 


In the University of California Botanical Garden over 60 species and 
varieties of the genus Nicotiana have been under investigation, together 
with species of such related genera as Petunia, Fabiana, Salpiglossis and 
Nierembergia. Morphological, distributional and cytogenetic evidence 
indicates the desirability of recognizing among the races of Nicotiana 
under observation, the following 40 species. The list is arranged in order 
of date of description. 


glutinosa L. solanifolia Walp. sylvestris Speg. & Comes 
paniculata L. rotundifolia Lindl. Stocktoni Brandg. 

rustica L. corymbosa Remy Sanderae hort. 

tabacum L. Miersii Remy Debneyi Domin 

tomentosa Ruiz & Pavon Cavanillesii Dun. Raimondii Macbride 
undulata Ruiz & Pavon trigonophylla Dun. nesophila Johnston 
longiflora Cav. cordifolia Phil. tomentosiformis Goodspeed 
plumbaginifolia Viv. wigandioides Koch. & Fint. 4 unnamed Australian species 
bonariensis Lehm. megalosiphon van Heurck 

repanda Willd. & Miill. 

suaveolens Lehm. attenuata Torr. 

Langsdorffii Weinm. Bigelovii S. Wats. 

alata Link & Otto nudicaulis S. Wats. 

glauca Grah. caudigera Phil. 


acuminata (Grah.) Hook. pampasana O. Kuntze 


The haploid chromosome numbers characteristic of these 40 species are 
9, 10, 12, 16, 20, 22, 24 and 32, and in the following table the species are 
arranged under these various chromosome number classes. 

For the sake of completeness, chromosome numbers of species earlier 
published by the author and by others are included in table 1 along with 
those first reported in the above list. Our representatives of the species 
listed have received detailed taxonomic study to assure nomenclatorial 
accuracy. For reasons to be discussed elsewhere, it appears desirable to 
distinguish plumbaginifolia from longiflora, caudigera from acuminata, 
Raimondii from glauca and nesophila from Stocktoni. The status of 
wigandioides and tomentosiformis has already been commented upon 
(Goodspeed”*). 

As indicated in footnote 2 to table 1, separate discussion of the Nicotiana 
species peculiar to Australia proper, Lord Howe Island and Rarotonga 
will be forthcoming. It should be noted, however, that decision to recog- 
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TABLE 1 
CHROMOSOME NuMBERS OF Nicotiana SprEcrIEs! 

9 10 12 16 20 22 24 32 
alata longiflora acuminata suaveolens megalo- rotundi- Bigelovii unnamed 
bonariensis plumbagi- attenuata and 3 siphon? folia? Debneyi? species? 
Langsdorfhi nifolia caudigera unnamed nesophila 
Sanderae Cavanillesii species? nudicaulis 

cordifolia repanda 
corymbosa rustica 
glauca Stocktoni 
glutinosa tabacum 
Miersii 
pampasana 
paniculata 
Raimondii 
solanifolia 
sylvestris 
tomentosa 
tomentosiformis 
trigonophylla 
undulata 
wigandioides 
4 2 19 4 1 1 8 1 


1 Expressed as number of pairs of chromosomes and based upon counts in P. M. C. and in root-tips. 
2 Australasian species concerning the cytotaxonomy of which an article by Miss Helen-Mar Wheeler 
is in preparation. 


nize as new species 4 races allied to suaveolens and obtained from the wild 
state came only after thorough comparative study of all available her- 
barium material of types. Of special interest are two 32); races, one of 
which represents an instance of polyploidy within the species, being a 
replica of typical swaveolens on a large scale, while the other is apparently a 
distinct species. 

With the exception of N. sylvestris and paniculata, all the species 
listed above give evidence, from herbarium specimens or in the living 
state, of being more or less polymorphic. As a result, named varieties of 
many of them are recognized (cf. Millan‘), but often, particularly in 
botanical gardens, varieties, and in some cases minor ones, are given species 
designations. It might be noted in this connection that a considerable 
number of varieties of tabacum, rustica, Bigelovit, acuminata, longiflora 
and Cavanillesti examined show no variation in chromosome number. 

Other authorities for chromosome number determinations are cited by 
Tischler® and Gaiser.6 The only species not in the above list (table 1) 
mentioned by them are cerinthoides, noctiflora, petiolaris, Forgetiana, 
viscosa, Palmeri, angustifolia, Rusby: and clarionensis. ‘The species re- 
ferred to by others as Rusbyi is listed here as tomentosiformis (Goodspeed”*). 
The designation clarionensis was never published (cf. Johnston’), this 
species being referred to here as Stocktoni. De Vilmorin and Simonet*® 
reported cerinthoides, Forgetiana and noctiflora as possessing 9;; and ¢rt- 
gonophylla, solanifolia and petiolaris as 24;;. Seed of what was grown by 
them under these designations (except petiolaris) was kindly supplied by 
M. de Vilmorin and when grown here gave the following results: ‘“‘nocti- 
flora” = typical, blue-pollened Langsdorffii; ‘‘Forgetiana’”’ and ‘‘cerin- 
thoides’”” = apparently derivatives of Langsdorffii X alata; ‘‘solanifolia’’ 
and ‘“‘trigonophylla’”” = rustica varieties. N. petiolaris is a designation 
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usually applied to certain rustica or tabacum varieties. Gaiser® lists 
“angustifolia” (n = 10) as reported by Clausen.’ Although the taxonomic 
situation in regard to this species designation is confused, it is clear that 
the race referred to by Clausen is a longiflora variety. The following 
species for which Christoff’ gives chromosome number determinations 
are also listed by Gaiser: “viscosa (?)” (n = 24), ““Palmeri (?)” (n = 12), 
and ‘‘solanifolia Wolf—cardiophylla Ph.’’* (n = 12). There is no evidence 
as to the identity of Christoff’s ‘‘viscosa’’ since if it had been an attenuata 
variety, as he suggests, the chromosome number should have been n = 12 
rather than ~ = 24. Christoff refers ‘‘Palmeri’’ to trigonophylla. 

Some 50 interspecific hybrids in Nicotiana have been obtained here and 
have included three amphdiploid types. To digluta (n = 36), an amphi- 
diploid glutinosa X tabacum hybrid (cf. Clausen and Goodspeed"), can 
now be added the amphidiploid glutinosa tomentosa hybrid, diglutosa 
(n = 24), and disualovit (n = 40), an amphidiploid suaveolens X Bigelovii 
hybrid. If the chromosome numbers of these three “‘species’’ are added 
to the series given in table 1, the complete list is 9, 10, 12, 16, 20, 22, 24, 
32, 36 and 40 pairs. 

Among by-products of recent cytogenetic investigations in Nicotiana 
carried on here are numerous aneuploid types in various species—alaia, 
Langsdorffii, sylvestris, paniculata, suaveolens, tabacum and rustica. Tetra- 
ploid Langsdorffit, sylvestris and tabacum have been studied. Haplonts 
within all chromosome number classes doubtless occur; those at present 
known are in the 9;;, 12;; and 24;; groups, i.e., glutinosa, tabacum, nudi- 
caulis; and Langsdorffii (Kostoff!?). A variety of quantitative chromo- 
some alterations have been established in stable derivatives of tabacum 
and sylvestris as products of treatment with high frequency radiation. 
A number of sesquidiploid forms have also been obtained (cf. Webber"*). 

The chromosome numbers of related genera have in a few cases been 
studied. Petunia species have for some time been known to possess 77, 
and more recently a 9;; species has been reported by Ferguson and Cool- 
idge.'4 This last is of obvious interest since a group of Nicotiana species 
shows this number and because hybrids between members of it and Petunia 
have been reported (cf., however, Kostoff*). The count of 22;; for Salpi- 
glossis sinuata reported by de Vilmorin and Simonet® can be confirmed. 
For Fabiana imbricata and Nierembergia frutescens, the numbers, not before 
reported, are in each case 2m = 18. 

Evidence which has for some years been accumulating will make possible 


* Concerning Christoft’s “‘solantfolia Wolf—cardiophylla Ph.’’ there is only the matter 
of the authority for solanifolia to be noted. He apparently misinterpreted ‘“‘Wolf” 
for ‘‘Walp.” on the herbarium sheet from which the seed was obtained, and this mistake 
is referred to because it has been perpetuated by Gaiser® where, in addition, the author- 
ity for cardiophylia is stated to be “‘Rh.” instead of Ph. 
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a partially complete description of species origins and relationships in the 
genus Nicotiana, based upon morphological, cytogenetic and other 
studies. A few postulates which appear justified and are pertinent in the 
present connection may be mentioned: 1. Certain genetic groups 
characterized by similarity both in chromosome number and in certain basic 
morphological features have been recognized. The ‘alata group’’ includes 
alata, Langsdorffii, Sanderae and bonariensis (n = 9); longiflora and plum- 
baginifolia (n = 10); the ‘‘tomentosa group” (” = 12) consists of tomentosa, 
tomentosiformis, glutinosa and wigandioides; the ‘‘glauca group” (m = 12) 
of glauca, Raimondii, paniculata, solanifolia, cordifolia and undulata; 
the ‘‘corymbosa group” (” = 12) of corymbosa and Miersii; the ‘‘Cava- 
nillesii group” (n = 12) of Cavanillesti and pampasana; the ‘“‘acumi- 
nata group” (m = 12) of acuminata, caudigera and attenuata; the ‘‘repanda 
group” (m = 24) of repanda, Stocktoni and nesophila. The ‘‘suaveolens 
group,” consisting of suaveolens, Debneyi, rotundifolia and megalosiphon, 
is.a natural one despite the distinctions in chromosome number involved. 
2. Certain, if not all, of the 24;; species are of allopolyploid origin. Thus 
tabacum (n = 24) may be traced to derivatives of hybrids between the 
“tomentosa group” (m = 12) and sylvestris (n = 12); rustica to derivatives 
of hybrids within the “paniculata group” (m = 12), Bigeloviit obviously 
involves in its origin attenuata (n = 12) and probably trigonophylla (n=12), 
while the origin of the “repanda group”’ and nudicaulis is not so clear, due 
perhaps to the absence of living descendants of their progenitors. The 
occurrence in our cultures of an amphidiploid between two 12;; species— 
diglutosa—lends some weight to the postulated amphidiploid origins just 
mentioned. That higher chromosome races of Nicotiana can be derived 
is indicated by our digluta and disualovii and by the results of Kostoff,® 
Lammerts'*—'8 and Rybin.” 
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YOUNG'S MODULUS AND POISSON'S RATIO WITH 
REFERENCE TO GEOPHYSICAL APPLICATIONS* 


By W. A. ZISMAN 
DEPARTMENT OF GEOLOGY AND GEOGRAPHY, HARVARD UNIVERSITY 


Communicated May 20, 1933 


Scope of the Investigation.—The theory of elasticity furnishes the follow- 
ing relations among the constants describing the elastic behavior of an 
isotropic and homogeneous medium obeying Hooke’s law: 


mr. (1 — a) 
p (1 + a)(1 — 2c)’ 


where V is the velocity of a longitudinal elastic wave in an infinite volume 
of the material, v is the velocity of a transverse wave, E is Young’s modulus, 
a is Poisson’s ratio and p is the density of the medium. In rocks, V and v 
can be determined by seismological methods. Assuming perfect elasticity 
and isotropy for rocks, E and o can be calculated. Conversely, on the 
same assumption one can predict V and v, if he can measure E and o by 
statical methods in the laboratory, at the range of stresses involved in the 
propagation of earthquake waves. 

The validity of the underlying assumption can be tested. With the 
explosion-wave method V and v are measured in an exposed terrane and 
thence the elastic constants are calculated. These are then compared with 
the constants measured statically on specimens of rock from the same 
terrane and at the corresponding extremely low range of stress-differences. 
If the respective values do not agree, the rock is not ideally elastic or 
homogeneous, and corrections must be applied before the constants de- 
termined from V and v, measured seismologically, can be used for identi- 
fying the rock where hidden underground. In practice, the elastic co- 
efficient useful in such diagnoses is cubic compressibility. Hence the geo- 
physicist needs data concerning the statical compressibility of an extensive 
series of standard rocks at varying hydrostatic pressures and under the 
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minute additional stress corresponding to that involved in wave propaga- 
tion. Here published data are incomplete and furnish a poor foundation 
for geophysical research, especially for the outer layers of the earth. 

The wide range of published values of E and o for rocks may be illus- 
trated by the case of granite. 

The Annual Report of the Watertown Arsenal for 1895 gives values of 
E for this rock with limits between 3 X 10"! and 6.5 X 10! dynes per 
square centimeter; and of o with limits between 0.07 and 0.25. These 
represent statical tests of American granites when subjected to unilateral 
pressure in air. 

Adams and Coker! found for American granites E ranging from 4.0 X 
10" to 6.3 X 10", with o between 0.20 and 0.26. On page 9 of their report, 
they mention Nagaoka’s values of E for granites, ranging from 1.1 X 
10"! to 4.2 X 104. Adams and Coker, like the experimenter at the Water- 
town Arsenal, used a statical testing machine, whereas Nagaoka worked 
with the flexure method at loads relatively and absolutely small. 

Again, Bach? states that E for German granites is 2.4 101}. 

Similar discrepancies are found in the case of practically every other 
type of rock investigated, and are clearly explicable not only by real 
differences in the specimens of each of the rock species, but also by differ- 
ences in the methods of measurement. Quantitative study of both of 
these controls over results for E and oc, as well as new measurements of 
cubic compressibility, were suggested by the Harvard Committee on 
Geophysics. The conclusions from this study, so far as they relate to E 
and a, are described in the present paper. 

A pparatus.—The study of the elastic behavior of most rocks is rendered 
difficult by the low stresses necessary in order to avoid departure from 
Hooke’s law. It is therefore essential to use apparatus of the greatest 
possible sensitivity for the measurement of the strains. Instruments 
were consequently built to measure E and o by a compression test in 
which the applied stress rarely exceeded one thousand pounds per square 
inch. This meant that, if the measured values of the elastic constants 
were to be accurate to one per cent, the instruments must be capable 
of detecting strains of the order of magnitude of one part in ten million. 

The heterogeneity and coarseness of texture expected in most rocks 
made it wisest to employ samples of a size sufficient to guarantee results 
fairly typical of the respective terranes. The specimens were rock cores 
two inches in diameter and about ten inches in length. Each rock was 
ground with emery in order to form a good right cylinder. The ends of 
the cylinder were compressed in the jaws of a testing machine and the 
change in diameter and change of length between two points on the cyl- 
inder wall were measured as a function of the load. From the changes 
in length and diameter E and o are directly calculated. Cylinders were 
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cored from large blocks in directions mutually perpendicular to one an- 
other, permitting a study of the variation of the elastic properties with 
direction. The expression “‘hardest splitting way’’ used in table 4 means 
that the axis of the cylinder is in the plane of most difficult cleavage. 

The gauge constructed to measure E was a modified form of the well- 
known Bauschinger rocking-mirror strain apparatus. It measured the 
change in length of the rock between two points twelve centimeters apart. 
The sensitivity attained was such that vibrations in the laboratory often 
caused great trouble. Hence many of the experiments were made at 
night. With a telescope magnifying thirty-five times to observe the 
mirror deflections, an illuminated scale at a distance of 4.5 meters from 
the mirror and a rocking knife-edge but one millimeter wide, a change of 
length of the specimen of 10-7 cm. caused a millimeter deflection in the 
telescope. 

The reliability of this strain gauge was proved by comparing the result 
for a cylinder of soft steel with the value of E obtained by Professor G. W. 
Pierce on the same cylinder by a magnetostriction method of great ac- 
curacy. The two measurements agreed within one per cent. Another 
good check was obtained with a cylinder of stainless steel. The observed 
E of soft steel was 21.7 X 10! and of stainless steel 22.1 X 101! dynes/cm.?. 
It is a pleasure to thank Professor Pierce for his measurements. 

The apparatus for the measurement of o is an improved form of a well- 
known type, utilizing a rocking mirror and a mechanical lever magnifica- 
tion. A detailed description will be found in the ‘‘Review of Scientific 
Instruments,” June, 1933. 

The improvements incorporated in this gauge make it capable of measur- 
ing o to one per cent, and yet the instrument is rugged and easily manipu- 
lated. The sensitivity was such that with the same telescope and scale-to- 
mirror distance used for the Young’s modulus experiment it was possible 
to obtain a scale deflection of 1 mm. for an increase of 10~’ cm. in the 
diameter of the specimen. 

The behavior of the Poisson’s ratio gauge was checked by measuring 
the lateral strains in the same specimens of steel used to check the Young’s 
modulus apparatus. In three successive and independent runs on soft 
steel the values of 0.296, 0.303 and 0.304 were obtained. The last check 
was made four months after the first two, during which time most of the 
rocks were studied. Another check was obtained by computing the cubic 
compressibility, 8, of the steel, using the measured values of E and o and 
the relation from the theory of isotropic elastic bodies 


3 
= ae — 20). 


The result agreed within five per cent with P. W. Bridgman’s published 
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measurement of the cubic compressibility of soft steel under hydrostatic 
pressure. The agreement is fair, especially in view of the uncertainty that 
the rolled steel cylinder was truly isotropic. 

The compression on the specimen ends was obtained with a new Bald- 
win-Southwark ‘“‘Emery’’ testing machine of sixty thousand pounds 
maximum loading capacity. This was placed at our service through the 
courtesy of Professor A. Sauveur of the Harvard Department of Metal- 
lurgy. 

The rock densities given in table 4 were obtained from the dimensions 
and mass of the rock cylinders. 

Measurement of Young’s Modulus (E).—Elastic hysteresis. If the strain 
- cis plotted along one axis of 
coérdinates and the load is 
plotted along the other, the 
curve for decreasing loads 
never coincides with that for 
increasing loads. A_ typical 
hysteresis loop for surface 
Quincy granite, specimen No. 
1, is shown in figure 1. If a 
succession of such hysteresis 
loops is made on the same 
specimen, it is found that the 
width of the loops decreases 
rapidly in the first two cycles 
and then becomes fairly con- 
FIG. 1 stant. The method of obtain- 
A, oman (ing Yougs motile fom 

ee such a graph was as follows: 
To find E for a mean load of 
500 pounds, the strain was 
” - oo ‘measured successively for loads 
ene ree of 0, 250, 500, 750, 1000, 900, 
750, 500, 250 and 0 pounds. This cycle of loads was carried through 
about five times in the space of two hours. The graph for any one of the 
last three runs may be represented by figure 1. In general, if the value 
of E was desired, corresponding to a hysteresis cycle for which the maxi- 
mum load was F; and the minimum load was F,, the following procedure 
was adopted (see Fig. 2). A series of five or six hysteresis cycles was made, 
each requiring ten to fifteen minutes tocomplete. For each cycle the load 
was adjusted in steps commencing at F,, increasing in four to six steps 
to a maximum load F:, and decreasing in four to six steps until F, was 
reached again. 
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Permanent set was shown by the lack of coincidence of points O, and 
O,. Elastic after-effect was shown by a time change of the position of 
point O2. It was always true that O, approached closer to O; with passage 
of time, although O, rarely settled down to coincidence with 0. 

It was always found true that at least half the paths of the curves, say, 
from A, to Az and B, to Be, were practically straight lines which differed 
in slope by from one to five per cent. The mean slope of the lines A:A2 
and B,Bz was used to calculate Young’s modulus. This value of E was 
then identified in this report as that corresponding to the mean load 


: 
9 (Fi + F). 


Thus for the cycle shown in figure 1, with the load F; at zero, and the 
loads at Ai, Az and C, re- 
spectively, 250, 750 and 
1000 pounds, the mean of 
the slopes of A,A, and B,B, 
was used to compute the 


FIG. 2 





1 
value of E at the load 3 


(F, + Fe) = 500 pounds. 
This corresponds to a mean 
load of 11.2 kg./cm.?. 

The value of E obtained 
in this manner was found 
to be the same for all 
hysteresis loops after the 
second. G 4+ Fe Fo 

Variation of E with Mean loed # —> 
Stress—It was found that 
E varies greatly with the mean stress used. Thus the following data 
were obtained for Quincy surface granite, cylinder No. 2: 


—> 


al 
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TABLE 1 
Fi, F:, MEAN STRESS, . 
POUNDS POUNDS KG./CM.? DYNES/CM.? 
0 500 2.8 2.39 X 101 
0 1000 11.2 2.58 
500 1500 22.4 3.07 
1000 2000 33.6 3.64 
2000 3000 56.0 4.05 
3000 4000 78.4 4.73 


The results for a Rockport granite specimen with axis in plane of “easiest 
splitting” were: 
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TABLE 2 
MEAN STRESS, E, 
KG./cM.? DYNES/CM.? 
0 3.48 X 10! 
2.2 3.95 
22.4 4.10 
44.8 4.35 
67.2 4.60 
89.6 4.85 
112 5.10 
132 5.35 
157 5.60 
179 5.85 
202 6.10 


In table 2 the value of E for zero load was obtained by extrapolating 
on a graph of mean load against E. 

The great variation of E with different loads is the cause of much of the 
apparent disagreement among the results of previous experimenters. 
The high values obtained by Adams and Coker are due to the fact that 
they measured E by taking hysteresis loops for a range of loads from a 
minimum of 1000 pounds per square inch to a maximum of 9000 pounds 
per square inch. When their result for granite is compared with the one 
here reported for a mean stress of 4500 pounds per square inch, the agree- 
ment is quite satisfactory. The low values of Nagaoka, mentioned on a 
preceding page, are at least partly due to his use of much smaller mean 
stresses. 

The final results for E at two values of the mean stress are given in 
table 4. 

It was found that those rocks which showed the smallest change of E 
with stress had the least porosity. 

The range in values of Z for the rocks investigated is considerable— 
from about 1.5 X 10! for marble to 10.0 X 10" for trap rock. The rocks 
with the largest values of E show the smallest variation of E with load. 
This seems related to the porosity. Rocks free from cracks and pores 
should exhibit little change of E with stress. One would expect rock at 
great depths in the earth to be in such a state We cannot regard the 
value of E, derived on a surface specimen for a mean stress of 100 kg./cm.? 
with this static apparatus, as equal to that for the same type of rock at 
a depth in the earth where the rock pressure is 100 kg./cm.?, for at that 
depth the porosity may be much less. Again, the mean stress of 100 
kg./cm.? used in the testing machine is a stress applied only over the ends 
of the rock cylinder, whereas in the earth at great depths the pressure 
on the rock is more nearly hydrostatic. If such loads as are available 
with the testing machine could completely close the pores and cracks 
in the rocks, we might then obtain a limit ng value of E as the mean 
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stress used is increased in steps. The limiting value might then be identi- 
fied as that existing at depths in the earth great enough to insure true 
compactness of the rocks. Actually, all those tested, except Vinal Haven 
trap, still showed increase of E up to the highest mean stresses applied 
(3000 pounds per square inch). 

The data of Adams and Coker and of others should therefore be used 
only for the range of stresses that they employed. 

Anisotropy of Rocks —The variation of E with direction in the rocks 
was large in some cases. This was especially true of Vermont marble, 
layered gneiss and limestone. On the other hand, we find in Vinal Haven 
trap and in Quincy granite from a depth of 235 feet almost perfect isotropy. 
It should not be concluded, however, that only layered or bedded rocks 
show great anisotropy, for French Creek norite is decidedly anisotropic. 

Measurements of Poisson’s Ratio.—Elastic hysteresis. There was much 
less hysteresis in the stress-strain graphs for the lateral strains as com- 
pared with that obtained for the longitudinal strains. Often there was 
but a trace of hysteresis and evidence of elastic after-effect or permanent 
set was lacking. This difference in behavior is, for the most part, due to 
the fact that E was measured first, and afterward the same rock cylinders 
were used to measure o. Obviously, then, the rock had been thoroughly 
pressure-seasoned before the commencement of the determination of o. 

In measuring o, the same method as that for Young’s modulus, namely, 
repeating hysteresis cycles and obtaining the mean slope of the stress- 
strain cycle, was employed. 

Variation of Poisson’s Ratio with Mean Siress—It was found that 
Poisson’s ratio also varies greatly with the mean stress applied. Table 3 
gives the results for the cylinder of Rockport granite whose values of E 
as a function of mean stress are listed in table 2. 


TABLE 3 
MEAN LOAD 
IN POUNDS POISSON’S RATIO 
0 0.080 
1000 0.088 
2000 0.097 
3000 0.105 
4000 0.114 
5000 0.122 
6000 0.131 
7000 0.141 
8000 0.151 
9000 0.162 


Here the value of Poisson’s ratio for zero load was obtained by extrapo- 
lation on a graph of o as a function of mean load. 

This rapid variation of « with mean load clearly explains the differences 
in the results of various experimenters. Those using large mean stresses 
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obtained high values of c, while those using small mean stresses obtained 
much lower values. It is for this reason that Adams and Coker obtained 
o for granite not far from 0.25, whereas the Watertown Arsenal obtained 
values around 0.10. 

Summary of Results on Poisson’s Ratio.—Each rock studied was stressed 
under two ranges of loads; hence two values of o were found. These are 
listed in table 4, which also gives the respective densities. 

As in the case of the measurements of E, the rocks showing the greatest 
variation of o with mean stress are seen to be those having the lowest 
values of «. Also low values of E accompany low values of ¢. This allows 
one to predict a large compressibility for such rocks at low pressures. 

The remarks made on a preceding page on the geophysical significance 
of the variation in the values of E with mean stress also apply to these 
results on o. 

The variation of o with direction in the rock is parallel to that of the 
change of E. 


TABLE 4 


ae 
Younc’s Moputvs, Porsson’s RATIO AND DENSITY OF SOME AMERICAN ROCKS 
MEAN 


STRESS, E, DENSITY, 
MATERIAL SPECIMEN KG./cM.? DYNES/CM.? eo GRAMS/cCM.® 

Granite 1 3.1 Ae f eG | | cae 2.61 
Quincy surface 78.4 5.08 
Granite 2.8 2.39 
Quincy surface 2 78.4 4.73 rs 2.61 
Granite 11.2 2.13 
Quincy surface rf 56.0 3.98 ae 2.61 
Granite 11.2 2.38 
Quincy surface 8 56.0 4.65 toe 2.61 
Granite 11.2 5.18 
Quincy 100 ft. deep 1 56.0 5.82 ee ate 2.59 
Granite 11.2 4.95 
Quincy 100 ft. deep 2 56.0 5.66 pee 2.59 
Granite 11.2 4.85 
Quincy 100 ft. deep 3 56.0 6.77 ea 2.59 
Granite 11.2 3.48 
Quincy 235 ft. deep 4 56.0 4.70 iG 2.64 
Granite 11.2 3.45 0.092 
Quincy 235 ft. deep 5 56.0 4.63 0.130 2.64 
Granite “‘Hardest 

Way” 11.2 3.54 0.091 
Rockport 100 ft. deep 56.0 4.26 0.110 2.63 
Granite ‘‘Rasiest 

Way” 11.2 3.95 0.082 
Rockport 100 ft. deep 56.0 4.39 0.103 2. 
Norite 1 36.5 8.07 0.212 2.85 
Sudbury 2 36.5 7.96 0.224 2.85 
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TABLE 4 (continued) 
MEAN 


STRESS, E, DENSITY, 
MATERIAL SPECIMEN KG./CM.? DYNES/CM.?2 o GRAMS/cM.? 
3 36.5 8.05 tees 2.85 
4 36.5 9.36 ee 2.85 
Norite 11.2 7.27 0.154 
French Creek 1 56.0 7.67 0.160 3.05 
Norite 11.2 5.84 0.108 
French Creek 2 56.0 6.15 0.129 3.05 
Norite 11.2 6.72 0.141 
French Creek 3 56.0 7.03 0.146 3.05 
Gneiss 0.0 0.328 0.0283 
Pelham (across foliation) 1 9.0 0.43 0.0283 2.64 
Gneiss 11.2 1.42 
Pelham (parallel to foliation) 2 56.0 2.20 alles 2.64 
Gneiss 11.2 1.68 0.082 
Pelham (parallel to foliation) 3 56.0 2.55 0.139 2.64 
Diabase 11.2 10.20 0.257 
Vinal Haven 1 56.0 10.15 0.257 2.96 
Diabase 11.2 10.25 0.244 
Vinal Haven 2 56.0 10.04 0.244 2.96 
Diabase 11.2 10.17 0.261 
Vinal Haven 3 56.0 10.04 0.261 2.96 
Marble 11.2 3.43 0.134 
Vermont 1 56.0 4.60 0.180 3:71 
Marble 12 2.32 0.134 
Vermont (perp. to bedding) 2 56.0 3.99 0.180 2.71 
Marble 11.2 3.83 0.142 
Vermont (parallel to bedding) 3 56.0 4.95 0.209 2.71 
Limestone 11.2 3.37 0.148 
Pennsylvania 1 56.0 4.57 0.183 2.69 
Limestone 11.2 7.60 0.256 
Pennsylvania 2 56.0 7.50 0.286 2.69 
Limestone 11.2 6.06 0.168 
Pennsylvania 3 56.0 7.22 0.183 2.69 
Quartzitic sandstone 11.2 6.36 0.109 
Pennsylvania 1 56.0 6.65 0.104 2.64 
Soft steel 28.0 21.7 0.301 


Description of Rocks Studied 1n These Experimenis.—The following 
description of the rocks studied is based upon the results of petrographic 
examinations by H. J. Fraser, E. B. Dane, Jr. (Quincy granite, Rockport 
granite, French Creek norite, Sudbury norite, Vinal Haven diabase and 
Pelham gneiss), and C. S. Hurlbut (limestone, sandstone and Vermont 
marble). The thanks of the writer and of the Harvard Committee on 
Geophysics are due to these gentlemen. 
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Granite, Quincy, Massachusetts——The specimens employed are fairly 
typical of the main body of the Quincy granite, especially that part of it 
where field velocity determinations have been made by Leet and Ewing.’ 
All the specimens are moderately coarse, holocrystalline and of normal 
hypiodiomorphic-granular texture. Unusually large thin sections were 
made from fresh pieces taken at a surface ledge near the J. S. Swingle 
quarry in West Quincy, and from a depth of 235 feet in the quarry. The 
volume percentages of the constituent minerals in these two phases are 
given in the first two columns of table 5. 


TABLE 5 
VOLUME PERCENTAGES OF THE MINERALS 


QUINCY GRANITE, QUINCY GRANITE, ROCKPORT GRANITE 


SURFACE 235 FT. DEEP 100 FT. DEEP 

Quartz 12 26 28 
Feldspar: 

Microperthite 65 50 64 

Albite 16 13 Small amount 
Riebeckite 3 7 3 
Soda-amphibole 1 3 3 
Magnetite 2 1 
Titanite, kaolin sericite, leucoxene, 

limonite, calcite, etc. 1 1 2 


In the surface specimen from Quincy, the average diameter of the feld- 
spar and riebeckite grains is about three millimeters; the other specimen 
has a somewhat smaller grain. The average of Warren’s closely similar 
analyses of the Quincy granite, main phase, can'be used to give the chemical 
composition.* 

Granite, Rockport, Massachusetts——In chemical and mineralogical 
composition this granite is nearly identical with the Quincy granite, while 
the grain is a little coarser. Table 5, third column, gives the mineral 
composition of the rock. 

Norite, French Creek, Pennsylvania.—This rock, known to quarry men 
as “‘St. Peter’s trap,’’ is moderately fresh, coarse-grained, holocrystalline 
and hypidiomorphic-granular. The volume percentages of its constituents 
are given in table 6. 


TABLE 6 
NoriTE, FRENCH CREEK, VOLUME PERCENTAGES 
Plagioclase, Abs;Angs 51 
Augite (75 per cent hedenbergite) 32 
Hypersthene 15 
Magnetite : 0.5 
Biotite 0.5 


Secondary hornblende, leucoxene, limonite 
and quartz 1 
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The diameters of the feldspars reach 1.0 millimeter and average about 0.3 
of a millimeter; those of the augite reach 5.5 millimeters and average 
about 0.8 of a millimeter. The minerals show no sign of strain. No 
chemical analysis is available. 

Norite, Sudbury, Ontario.—On request, the officials of the International 
Nickel Company, Ltd., kindly presented a large number of diamond-drill 
cores from the norite at Sudbury, Ontario. Four of these, especially free 
from alteration, were selected for study. Their original depths below the 
surface are not known. The respective volume percentages of the constitu- 
ent minerals are given in table 7. 


TABLE 7 
NOoRITE, SUDBURY, VOLUME PERCENTAGES 
CORE 1 CORE 2 CORE 3 CORE 4 
Labradorite, AbysAnse “e 67 58 70 
Andesine, AbssAna, 74 is i = 
Potash feldspar 1 1 2 4 
Hypersthene (25 per cent MgSiOs) 15 8 18 10 
Augite 3 5 + 2.5 
Hornblende 2 5 8 3 
Biotite 1:6 6 1.5 3 
Quartz 3 8 8 6 
Magnetite, apatite, sericite, kaolin, 
carbonate 0.5 Nes 0.5 1.5 


All four specimens failed to show evidence of any internal strain. Cores 
1 and 4 are almost ideally fresh; the other two are slightly altered. The 
average of five analyses® of this norite illustrates fairly well the chemical 
composition. 

Diabase, Vinal Haven, Maine.—This fresh rock from the coast of Maine 
was chosen as nearly typifying fresh olivine diabase, one of the crystallized 
products from the quantitatively important basaltic magma. The volume 
percentages are given in table 8. 


TABLE 8 


DIABASE, VINAL HAVEN, VOLUME PERCENTAGES 


Labradorite, AbsoAneo 72 
Olivine 10 
Augite (75 per cent hedenbergite) 10 
Hypersthene 2 
Biotite 3 
Magnetite 2 
Uralite, antigorite, apatite, hornblende, 

quartz, limonite 1 


The texture is ophitic. The average lengths of the feldspars, olivines 
and augites are respectively 0.8 mm., 0.4 mm. and 0.6 mm. The thin 
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section shows absence of strain. No chemical analysis is available, but the 
composition must be near that of the average olivine diabase. 

Gneiss, Pelham, Massachusetts——Professor G. W. Bain of Amherst 
College was good enough to forward a cuboidal block of this fresh, notably 
foliated, but otherwise homogeneous granitic gneiss, which has an ex- 
tensive outcrop around Pelham, Massachusetts. Volume percentages 
of the minerals were measured in sections cut across the foliation and 
parallel to it. The results are given in table 9. 


TABLE 9 
GNEISS, PELHAM, MASSACHUSETTS 


PARALLEL TO FOLIATION ACROSS FOLIATION 

: VOL. PER CENT AVG. DIAM. VOL. PER CENT AVG. DIAM, 
Quartz, albite, orthoclase 73 0.2 mm, 84 0.3 mm. 
Microcline 13 0.5 8 0.3 
Biotite 12 0.7 6 0.25 
Apatite 0.5 0.5 
Epidote 1.0 1.0 
Carbonate 0.5 0.5 


As expected, the biotite, quartz and feldspars are elongated more or 
less systematically in the plane of the foliation. The minerals show 
sutured boundaries but not much intricate interlocking. 

The chemical composition is doubtless well represented by an analysis 
made of the gneiss at Ward’s quarry, Pelham.° 

Quarizitic Sandstone, Pennsylvania.—Blocks of this Cambrian rock 
from the Hardyston Quartzite formation, as well as those of the limestone 
next to be described, were obtained through the courtesy of Professor B. L. 
Miller of Lehigh University. The sandstone came from the north side of 
the mountain about midway between Allentown and Bethlehem, Pennsyl- 
vania—about 40° 37’ N. Lat., 75° 25’ W. Long. It is a strong, well-ce- 
mented rock approaching a typical quartzite in mode of fracture. The rock 
is comparatively massive, the bedding not being prominent in the speci- 
men studied. The volume percentages, maximum diameters and average 
diameters of the constituents are given in table 10. 

Nearly all of it shows strain shadows in polarized light. 


TABLE 10 
QUARTZITIC SANDSTONE, PENNSYLVANIA 

VOLUME, MAXIMUM AVERAGE 

PER CENT DIAMETER DIAMETER 
Quartz 80 3.0 mm. 0.3 mm. 
Microperthite 5 2:0 0.4 
Microcline 1 0.6 0.3 
Sericite 7 neRE. * 
Magnetite 7 0.5 0.07 
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Limestone (‘‘Cement Stone’), Pennsylvania.—This Ordovician rock also 
represents a formation occurring in Pennsylvania, sampled at the quarry 
of the Nazareth Cement Company, Nazareth, Pennsylvania—approxi- 
mately 40° 45.5’ N. Lat., 75° 18’ W. Long. The rock is dark-colored, 
carbonaceous and fossiliferous. The volume percentages of the minerals 
are: calcite, 80; carbonaceous material, 17; quartz, 2; feldspar, 1. 
The grain is extremely fine, the largest grains of calcite having diameters 
of 0.2mm. The abundant fossils are bryozoa. 

Marble, Vermont.—A large cube of this well-known white marble was 
kindly presented by the Vermont Marble Company of Proctor, Vermont. 
The rock is essentially a mass of irregular, interlocking grains of calcite, 
which almost everywhere in the thin section show pressure-twinning. 
The only other constituent, quartz, makes up no more than 0.5 per cent 
of the marble. The carbonate grains average 0.4 mm. in diameter and 
have a maximum of about 1 mm. 

Acknowledgments.—The author wishes to express his sincere thanks to 
Professor R. A. Daly and Professor P. W. Bridgman for their interest 
and valuable advice during the progress of the work. It is due to the 
courtesy of the Harvard University Department of Physics that the 
facilities of the new research laboratory were made available for these 
experiments. 


* Paper No. 3 published under the auspices of the Committee on Geophysical Re- 
search (P. W. Bridgman, R. A. Daly, L. C. Graton, D. H. McLaughlin and H. Shapley) 
and of the Division of Geological Sciences at Harvard University. 

1F. D. Adams and E. G. Coker, “Investigation into the Elastic Constants of Rocks,” 
Publication 46, Carnegie Institution of Washington, (1906). 

2C. Bach, “‘Elastizitat und Festigkeit,’’ 4th edition, 64 (1904). 

*L. D. Leet and W. M. Ewing, “Velocity of Elastic Waves in Granite,” Physics, 2, 
No. 3, March, 160-173 (1982). 

4C. H. Warren, Proc. Amer. Acad. Sct., 49, 227 (1918). 

5 A. P. Coleman, E. S. Moore and T. L. Walker, University of Toronto Studies, No. 
28, 18 (1929). 

6B. K. Emerson, Bull. 597, U. S. Geol. Survey, 1917, p. 249. 


i 
HM 
| 
| 
| 
4 
i 








666 GEOLOGY: W. A. ZISMAN Proc. N. A. S. 


COMPRESSIBILITY AND ANISOTROPY OF ROCKS AT AND 
NEAR THE EARTH’S SURFACE* 


By W. A. ZISMAN 
DEPARTMENT OF GEOLOGY AND GEOGRAPHY, HARVARD UNIVERSITY 


Communicated May 20, 1933 


Purpose of the Investigation.—That remarkably active branch of geo- 
physics, seismology, is dependent for theory and practice on accurate 
information regarding those physical properties of rocks that are most 
intimately related to the propagation of elastic waves in the earth. Most 
important of all for use in diagnosing the composition of hidden rocks is 
volume compressibility. Due to recent improvements in technique by 
Bridgman and others, it is now possible to measure the linear and cubic 
compressibility of a rock under hydrostatic pressures as high as 10,000 
to 30,000 kilograms per square centimeter. This is the range of pressure 
at depths of roughly 35 to 100 kilometers in the earth. While Bridgman 
has supplied some data on the compressibility of rocks at high pressure, 
the greater part of the available information comes from the work of Adams 
and Williamson at the Geophysical Laboratory of the Carnegie Institution 
of Washington.' 

Adams and Williamson worked in a pressure range of 2000 to 10,000 
“megabars’”’ (now sometimes called ‘‘bars’’), a megabar being one million 
dynes per square centimeter and equal to about 1.02 kilograms per square 
centimeter or about ninety-nine per cent of one atmosphere. Their 
measurements were made under two different conditions. In the one 
case the rock specimens were exposed directly to the fluid transmitting 
the pressure; in the other case the air-dried specimens were covered with 
thin, impermeable tin-foil which prevented the direct access of the fluid. 
Adams and Williamson pointed out that the measured compressibility 
should be the greater in the second case, provided the rock carries cracks 
or pores in communication. For clearly, when the pressure goes on, the 
cracks and pores should be more easily closed if the specimen is ‘‘covered”’ 
than if the transmitting fluid is allowed to enter cracks and pores and 
so add to the resistance to closure. This theoretical expectation agreed 
with the results of the many experiments, except in the case of three rocks: 
marble, Westerly granite and Stone Mountain granite. These anomalies 
were not explained. During the experiments reported in the present 
paper, twenty-eight different specimens of rock were tested in both the 
covered and uncovered states; in every case the behavior was that fore- 
told by Adams and Williamson. 

Those authors showed that the compressibility of most rocks decreases 
rapidly as the pressure rises from an initial amount of three hundred 
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atmospheres. This fact was discovered when they compared their values 
of compressibility at pressures over two thousand atmospheres with those 
obtained by Adams and Coker, who calculated the compressibilities of their 
rocks from measurements of Young’s modulus and Poisson’s ratio at a 
pressure range of about seventy to six hundred atmospheres.? From 
the two sets of results Adams and Williamson constructed curves of com- 
pressibility, reaching back to zero pressure, but their method for the lower 
pressures was by frank extrapolation. For seismological problems relating 
to rock formations at and close to the earth’s surface there has been 
evident need of a direct experimental study of rocks in the pressure field 
below that occupied by Adams and Coker. The following paper presents 
the results of such a study. 

Moreover, the geophysicist needs to know how far rocks, especially 
those at shallow depths, depart from isotropy. That the departure may be 
great has been shown by recent measurements of Young’s modulus and 
Poisson’s ratio at pressures reaching about fifty-six atmospheres.* The 
same fact has become evident in the course of the present research, when 
cubic compressibilities were determined by first measuring linear com- 
pressibility at a pressure range of one to about eight hundred atmospheres. 
The amount of the anisotropy in standard rocks can be gauged by study 
of the tables of data to be seen on later pages. It is clear that the seismolo- 
gist must consider anisotropy when he is dealing with waves passing 
along the surface layer of the globe. 

Description of Piezometer—The piezometer designed for this work is 
on the same principle as Bridgman’s well-known lever-piezometer used 
at pressures ranging up to twelve thousand atmospheres.‘ For research 
in that field of pressure the instrument could manifestly accommodate 
only short columns of the solids to be investigated. To secure reliable 
results at low pressures it was necessary to use columns six or eight times 
as long, the length actually favored being eight inches (twenty centimeters). 
The width was regularly five-eighths of an inch (1.6 centimeters). Such 
great lengths and actual thicknesses secured against errors due to the 
excessive influence of specially large crystals in the rocks, and also assured 
high sensitivity for the apparatus. 

The principle of the Bridgman instrument, built with specifications 
to match the new kind of investigation, may be briefly reviewed (see 
Fig. 1). One end of the rock cylinder rests on the piezometer case at the 
point A. The other end is in contact at B with the short arm CB of a lever 
ECB whose fulcrum is at C, the end of the pointed screw DC. The arm 
EC is about fourteen times longer than arm CB. A spring JI maintains 
CB in firm contact with the upper end B of the specimen. At EH is a stiff 
straight piece of manganin resistance wire about one centimeter long. 
This is in gentle, sliding contact with a blade of manganin, G, which is 
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supported but electrically insulated from the main body of the piezometer 

M;M.2M,. At F and H are soldered wire connections for the electrical 

circuit. When the hydrostatic pressure is increased on the liquid in the 

steel cylinder into which the piezometer is screwed, the specimen shortens 
in all its dimensions. Point A is fixed, but point B moves downward. 

Then the lever arm ECB rocks about the fulcrum at C and point E moves 

to the right. Therefore, the distance FG decreases as the result of the sliding 

of the manganin wire over the fixed 
oe. 2 knife-edge contact atG. By measur- 

Z ’ 4 ing the consequent change in elec- 

Flat Spring trical resistance of FG, we can, 

from our knowledge of the resis- 

tance per unit length of the wire, 
calculate the decrease in length of 
the rock under pressure. 

Ms In order to eliminate troubles 
from the contact resistance at G, 
the resistance of FG is measured by 
a potentiometer bridge method de- 
scribed in Bridgman’s papers. For 
this method F and G serve as po- 
tential and E and H as current con- 

Piezometer nections. 

— It is apparent that such an instru- 
ment measures only the difference 
in the hydrostatic linear compressi- 
bilities of the rock and the steel of 
the piezometer case, for when the 
length of the specimen shortens so 
does the distance M2M3, due to the 
compressibility of the steel. The 
true compressibility of the rock 

a ee is obtained by allowing for that of 

Mm, A Me the steel, for which Bridgman’s 

SCHEMATIC DIAGRAM OF PIEZOMETER oe 2.00 X 10-1 cm.*/dyne, was 
used. 

In order to test the reliability of the instrument, a cylinder of soft steel 
of the same kind as in the piezometer jacket was used in place of the rock 
specimen. No change in resistance under pressure was observed that 
could be due to relative motion of the specimen and the jacket, such as 
would occur if distortion took place in the instrument when the high 
pressure permitted any internal stresses in the steel to exhibit their pres- 
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If a specimen of some homogeneous substance free from internal strains 
and cavities is used, the graph of shortening of the specimen as a function 
of pressure will be a smooth curve free from hysteresis. This is the most 
reliable test for the existence of slipping at the bearing points B and C. 
A cylinder of pure cast lead was used for such a test. Several independent 
runs gave the same result and proved the absence of hysteresis. At thirty 
degrees centigrade the linear compressibility of lead was found to be 8.40 X 
10-3 cm.?/dyne. Bridgman® obtained 8.07 X 10-1* cm.?/dyne. This 
difference may be due to the lack of perfect annealing in my specimen. 

The hydrostatic pressure was measured with a manganin resistance 
coil by inserting a 250-ohm coil of well seasoned manganin wire in the 
pressure chamber and measuring its change of resistance under pressure 
with a Carey-Foster bridge. The pressure calibration of this type of 
secondary gauge is simple. The resistance of manganin was shown by 
Bridgman to change linearly with pressure in the range 0-10,000 kg./cm.?. 
Therefore it is necessary only to measure the resistance of the coil at 
atmospheric pressure and again at some higher known pressure, such as 
that at which mercury freezes at zero centigrade (7640 kg./cm.”). There 
was no change in the calibration of the coil used during the year needed for 
these measurements. 

Preparation of Specimens.—The rock cylinders were core-drilled from 
the same blocks as those from which were cut the two-inch-diameter 
cylinders used for the measurement of Young’s modulus and Poisson’s 
ratio described in a previous paper.* The same method of indicating the 
directions of coring of the cylinder for compressibility tests was used. 
The cores were ground carefully to form accurate right cylinders. 

The compressibility of the specimen was measured when immersed in 
kerosene, the fluid transmitting the pressure. This will be referred to as 
the compressibility of the ‘“‘uncovered” rock. In contrast is the compressi- 
bility of the ‘‘covered’”’ rock, the cylinder in this case having been enclosed 
in an impermeable sheath of 0.002-inch annealed copper foil. This cover- 
ing was strong enough to prevent entrance of the kerosene into the rock 
and flexible enough to transmit the pressure freely to the rock within. 

The method of preparing these covers may be indicated. The rock 
cylinder was rolled in a sheet of foil, much as the tobacco of a cigarette 
is rolled in its paper cover. The end of the rolled sheet of foil then formed 
a seam running the length of the cylinder. This seam was soft-soldered, 
and to facilitate the operation the foil edge was tinned before rolling. 
Thin caps were made for covering the ends of the cylinder, by spinning a 
sheet of copper into a cap whose inside diameter was five-eighths of an 
inch and whose length was from a quarter to three-eighths of aninch. The 
cap fitted snugly over the end of the specimen and the rim was soft- 
soldered to the copper covering beneath. When hydrostatic pressure was 
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applied, the excess foil present was forced into the surface cavities in the 
rock, making an extremely close fit like that of the skin of an onion. 

After being covered, the rock was exposed to a pressure of 1000 kg./cm.? 
for five minutes. Then a small hole was made in the cover. If there had 
been a leak, it was detected by the issuing of small bubbles of kerosene 
from this hole. After some practice in the technique it was found that few 
covers leaked. After this test, the small hole was sealed with soft solder, 
and the specimen was ready for a compressibility measurement. Upon 


40r 


35Fr 


30F 


nN 
wo 


\ 
° 
T 


a 
T 


5 
T 


FIG 2 


Shortening of Cylinder Under Pressure 


uw 
T 








- : 3 i tg 

Pressure —> 
completion of an experiment the foil was stripped from the rock to permit 
detection of further evidence of a leak. 

A blank run on the compressibility of a cylinder of soft steel covered 
with foil was made, and it was thus possible to prove that there was no 
contribution to the observed compressibility attributable to the existence 
of air spaces between the foil and the flat ends of the rock. When the 
compressibility of the covered rock was measured, the smallness of the 
permanent set found after each pressure run was considered good reason 
for believing the observed pressure effect to be due to the change in length 
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of the rock rather than to the squeezing of the foil covering against the 
ends of the specimen. This conclusion was strengthened by the fact that 
successive runs on the same rock using entirely new covers gave the same 
compressibility within a few per cent. Moreover, the compressibility of 
covered and uncovered obsidian glass was found to agree, as it should, 
since such material is non-porous. . 

Description of Results.—Figure 2 is the graph of the shortening of a 
cylinder of covered Rockport granite against pressure. There is a per- 
manent set amounting to about three per cent of the total shortening. 
The uncovered rocks show somewhat smaller amounts of permanent set. 
From such a graph it is not difficult to obtain measurements of the com- 
pressibility as a function of pressure by reading off the slope for, say, eight 
values of the pressure. 


FIG 3 
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Just before taking a compressibility measurement to 800 kg./cm.? a 
ten-minute pressure seasoning to 1000 atmospheres was given the appa- 
ratus and specimen. A repetition of the compressibility measurement was 
found unnecessary. Even when the pressure seasoning lasted as long as 
five hours, the compressibility was unchanged. 

The permanent set was found greatest for the rocks exhibiting the 
largest change in compressibility with pressure, that is, for the most 
porous rocks. On the other hand, pronounced hysteresis, as indicated 
by a great width for the hysteresis loop, was not always easily linked with 
high porosity. It is quite possible that most of the hysteresis was due to 
the presence of internal strains in the rocks, such as there must always be 
in a rock essentially composed of non-cubic crystals. 
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In the following tables the linear (a) and cubic compressibilities (8) 
at thirty degrees centigrade are given as a function of pressure for both the 
covered and uncovered rocks. The cubic compressibility was computed 
by tripling the mean value of the linear compressibilities. In figure 3 is 
presented the cubic compressibility of granite as a function of pressure. 

The difference in behavior of each rock with and without cover is strik- 
ing. The values for the uncovered rocks should be considered as the more 
accurately determined, chiefly because the small change in compressi- 
bility with pressure then observed made it easier to obtain from the graphs 
accurate values of shortening with pressure. The values of the com- 
pressibilities for the compact rocks are good to two per cent, while for the 
rocks showing rapid changes in compressibility with pressure, the error 
may rise as high as five per cent for the initial compressibility and two to 
three per cent at 840 kg./cin.’. 

Throughout the whole series of measurements the compressibility of 
each covered rock is greater than its value when uncovered. The difference 
between the two becomes smaller as the pressure is raised. This is what 
one would expect, for at high pressures the looseness of texture in the rock 
due to the presence of cavities and cracks is substantially removed. 
Thus we find the two values almost equal at 840 kg./cm.?. 

The change in the linear compressibility with direction parallels the 
vectorial change of Eando. This is a natural relation, the cylinders having 
been cut respectively from the same blocks of rock and in the same direc- 
tions. The results for French Creek norite will serve as an example. 
Of the three directions studied, the greatest linear compressibility occurs in 
specimen No. 2, and the least in No. 1. The measurements of E and o 
show number two to have the smallest value and number one the greatest. 
Thus, as we should expect, a large E and o accompany a small linear 
compressibility. 

Description of Rocks Used.—Most of the rocks used in these experiments 
were described in a previous publication.* The following descriptions of 
additional rocks studied here are due to the kindness of H. J. Fraser, E. B. 
Dane, Jr., and C. S. Hurlbut, of the Department of Geology of Harvard 
University. 

Obsidian, Lipari Islands—This well-known material is essentially an 
unstrained rhyolitic glass, containing thread-like crystallites of iron oxide 
and minute feldspars totaling about two per cent of the volume. Align- 
ment of the crystallites indicates a definite flow structure, which is macro- 
scopically not conspicuous. A chemical analysis of this glass is to be found 
in Washington’s Tables. 

Peridotite (Dunite), North Carolina.—A large block of this astonishingly 
fresh, ultrabasic rock was generously supplied by Mr. A. W. Hyatt of 
Balsam (Balsam Gap), North Carolina. A sugary, holocrystalline, allo- 
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TABLE 1 
LINEAR (a) AND CUBIC COMPRESSIBILITIES (8) AT 30°C. (Cm.?/Dyne) 
(Cylinders numbered, mutually at right angles to one another) 


PRESSURE a X 1018 a X 1018 a X 1018 B X 1013 
(Ke. /cm.?) cyL. 1 CYL. 2 cyL. 3 MEAN 


Quincy granite at surface (un- 


covered) 0 6.92 9.42 8.32 24.7 
60 6.76 9.05 8.12 24.0 
120 6.62 8.69 7.93 23.2 
240 6.30 8.27 7.62 22.2 
360 6.10 7.95 7.40 21.4 
480 5.93 7.78 7.25 21.0 
600 5.89 7.74 7.15 20.8 
720 7.73 7.09 20.7 
840 7.72 7.07 
Quincy granite, 100 feet deep 
(covered) 0 20.9 31.1 23.6 75.6 
24 18.9 25.3 20.4 64.5 
48 17.1 21.1 17.4 55.6 
72 15.8 17.4 15.5 48.7 
96 14.7 15.1 14.2 43.9 
120 13.6 13.7 13.0 40.2 
144 12.7 12.7 12.1 37.6 
180 11.5 11.6 11.2 34.4 
300 9.52 9.81 9.51 28.8 
600 8.26 8.54 8.44 25.3 
720 8.18 8.22 24.6 
Quincy granite, 100 feet deep 
(uncovered) 0 6.50 6.23 19.2 
60 6.35 6.15 18.8 
120 6.20 6.07 18.5 
240 5.88 5.93 42:7 
360 5.72 5.82 17.3 
480 5.58 5.72 17.0 
600 5.51 5.65 16.7 
720 5.49 5.63 16.7 
840 5.48 5.61 16.6 
Quincy granite, 235 feet deep 
(covered) 0 30.0 89.9 
24 26.3 78.8 
48 23.2 69.7 
7 20.8 62.3 
96 18.7 56.0 
120 16.9 50.8 
144 15.5 46.5 
180 13.8 41.3 
300 10.4 31.2 
600 8.32 25.0 
720 7.84 23.6 
Quincy granite, 235 feet deep 
(uncovered) . 0 7.68 6.35 31.13 
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Rockport granite, 100 feet deep 
(covered) 


Rockport granite, 100 feet deep 
(uncovered) 


French Creek norite (covered) 


French Creek norite (uncovered) 





TABLE 1 (Continued) 


PRESSURE 
(Ko. /cm.?) 
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180 
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360 
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600 
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48 
72 
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120 
144 
180 
300 
600 
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120 
240 
360 
480 
720 
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300 
600 
720 
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19.1 
17.4 
15.9 
14.6 
13.5 
12.5 
Sg 
10.6 
7.82 
5.40 
4.93 
4.43 
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4.38 
4.33 
4.30 
4.26 
4.24 
4.22 
4.20 
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TABLE 1 (Continued) 


PRESSURE 
(ko./cm.?) 


Sudbury norite (covered) 0 
24 

48 

72 

96 

120 

144 

180 

300 


720 
Sudbury norite (uncovered) 0 
60 
120 
240 
360 
480 
600 
Olivine diabase (covered) 0 
120 
300 
480 
720 
Olivine diabase (uncovered) 0-200 
400-720 
Peridotite (uncovered) 0 
60 
120 
240 
360 
480 
600 
720 
Orthogneiss (covered; cylinder 
No. 1 cut with axis at right 
angle to plane of schistosity) 0 
120 
300 
600 
720 
Orthogneiss (uncovered) 0 
60 
120 
240 


360 
480 
600 
720 


* No sensible change of @ in range of pressure. 
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CYL. 2 
a X 1018 
10.8 
10.16 
9.41 
8.80 
8.22 
7.75 
7.40 
6.93 
5.96 
5.57 
5.56 
5.42* 
5.42* 
5.42* 
5.42* 
5.42* 
5.42* 
5.42* 


2.98* 
2.98* 
2.98* 
2.98* 
2.98* 
2.98* 
2.98* 
2.98* 


6.32 
6.17 
6.03 
5.82 
5.65 
5.57 
5.52 
5.48 


cyL. 3 
a X 101% 


48.5 
27.4 
14.3 
9.61 
9.06 
.18* 
18* 
18* 
18* 
18* 
18* 
18* 
18? 
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BX 1018 
MEAN 


31. 
29. 
27. 
25. 
24. 


. 9 . . . 
RNIOCOCHWORAWANATAEHROHAAQN 


20.7 
20.0 
19.3 
18.4 
Laws 
17.2 
16.8 
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TABLE 1 (Continued) 


PRESSURE cYL. 1 
(kc./cm.?) a@ X 1018 aX 


Quartzitic sandstone (covered; 
cylinder No. 1 cut with axis 


at right angle to bedding) 0 19.7 23 
24 18.3 21 

48 i722 19 

72 16.3 18 

96 15.6 17: 

120 15.0 16 

144 14.5 15 

180 13.8 15 

300 12.1 13. 

600 10.9 11 

720 10.5 10 

Quartzitic sandstone (uncovered) 0 8.39* 9 
60 8.39* 9 

120 8.39* 9 

240 8.39* 9 

360 8.39* 8. 

480 8.39* 8 

600 8.39* 8 

720 8.39* 8 


Vermont marble (covered cylin- 
der No. 2 cut with axis at 


right angle to bedding) 0 58.8 75. 
24 42.0 47 
48 24.0 26. 
72 18.0 18 
96 14.2 15 
120 ai. 7 12 
144 10.2 10. 
180 8.79 9. 
300 6.67 ‘fe 
600 4.29 5. 
720 he 5. 
Vermont marble (uncovered) 0 4.17 5. 
120 4.06 5. 
240 3.90 5. 
360 3.79 5. 
480 3.69 5. 
600 3.61 5. 
720 3.53 5. 
840 3.46 5. 
Limestone (covered; cylinder 
No. 1 cut with axis at right 
angle to bedding) 0 14.5 6 
60 13.9 6. 
120 13.4 6 
240 12.4 5 


* No sensible change of a in the range of pressure. 





CYL. 2 


1018 


Proc. N. A. S. 
CYL. 3 B X 1013 
a X 1018 MEAN 
15.7 58.7 
14.3 54.1 
13.3 50.2 
12.4 47.1 
11.8 44.8 
11.3 42.8 
10.9 41.3 
10.4 39.2 
9.48 34.8 
8.92 30.9 
8.82 29.9 
8.86* 26.7 
8.86* 26.6 
8.86* 26.5 
8.86* 26.3 
8.86* 26.2 
8.86* 26.1 
8.86* 26.0 
8.86* 26.0 
46.5 180.0 
26.2 115.0 
18.0 68.2 
13.6 50.3 
10.9 40.9 
9.12 33.1 
8.32 29.2 
7.95 25.9 
6.59 20.5 
5.10 15.0 
4.53 13.9 
4.47 13.8 
4.40 13.5 
4.28 13.3 
4.09 12.9 
3.82 12.6 
3.50 12.2 
8.49 29.2 
8.29 28.2 
8.09 27.5 
7.75 26.1 
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TABLE 1 (Continued) 


PRESSURE crn. 1 CYL. 2 cy. 3 BX 1018 
(xc./cm.2) a X 1018 a X 1018 a X 1013 MEAN 
360 11.6 5.79 7.44 24.9 
480 11.0 5.73 7.25 23.9 
600 10.7 5.68 7.10 23.5 
720 10.4 5.65 6.94 23.0 
Limestone (uncovered) 0 LEG" 5.88* 7.22 24.7 
60 12:6? 5.88* 7.13 24.6 
120 11.6* 5.88* 7.06 24.5 
240 11:6" 5.88* 6.96 24.5 
360 11.6" 5.88* 6.81 24.2 
480 UD bg 5.88* 6.71 24.1 
600 15-6" 5.88* 6.58 24.1 
720 114;6* 5.88* 6.45 24.0 
840 116° 5.88* 6.32 23.7 
Dolomite (massive; covered) 0 14.2 10.60 37.1 
24 13.0 10.2 34.8 
. 60 10.6 9.63 30.4 
120 8.28 8.65 25.4 
240 6.62 6.73 20.1 
360 5.88 5.37 16.9 
480 5.40 4.91 15.5 
600 5.18 4.70 14.8 
720 4.99 4.44 14.2 
Dolomite (uncovered) _ 0 4.11 3.81 11.9 
360 4.11 3.81 11.9 
720 4.11 3.81 11.9 


* No sensible change of a in range of pressure. 


triomorphic mass of medium grain, it is composed of olivine (92 per cent 
MgeSiO, and 8 per cent Fe,SiQ,, as measured optically), 97 per cent; 
serpentine, two per cent; with chromite, carbonate and hornblende 
totaling one per cent. The olivines reach 4.4 mm. in diameter, the average 
being 0.75 mm. The chemical composition of this dunite must approxi- 
mate that of the olivine alone, as found by Gonyer.® 


Dolomite, Pennsylyania.—This material of Cambrian age, obtained 


through the courtesy of Professor B. L. Miller of Lehigh University, came 
from the Bethlehem Mines Corporation quarry just east of Bethlehem, 
Pennsylvania (about 40° 37’ N. Lat., 75° 22’ W. Long., Allentown topo- 
graphic sheet). The rock is massive, with obscure bedding in the block 
actually studied. The mineral data are as follows: 


MAXIMUM AVERAGE 
VOLUME, DIAMETER, DIAMETER, 
PER CENT MM, MM. 
Calcite 98 0.25 0.01 
Quartz 1 0.25 0.1 
Feldspar . 1 0.2 0.1 
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TABLE 2 
SUMMARY OF COMPRESSIBILITY DATA 
COMPRESSIBILITY 
PRESSURE 8 X 1038 
(KG./CM.2) UNCOVERED COVERED DENSITY 
Quincy granite, at surface 0 24.7 2.61 
720 20.7 ss 
Quincy granite, from depth of 100 feet 0 19.2 76.6 2.59 
720 16.7 24.6 
Quincy granite, from depth of 235 feet 0 21.1 89.8 2.64 
720 17.8} 23.6 
Rockport granite 0 19.5 91.7 2.63 
720 17.3 26.61 
French Creek norite 0 14.0 59.0 3.05 
720 12.5 15.2 
Sudbury norite 0 16.5 31.5 2.85 
720 15.7} 16.3 
Olivine diabase 0 14.6? iW Be | 2.96 
720 13.33 12.6 
Peridotite 0 1.2 °° fi 3.27 
720 9.37 oe 
Orthogneiss 0 20.7 2 2.64 
720 16.6 Se 
Obsidian 0 29.6 29.5 2.34 
720 29.6 29.5 
Quartzitic sandstone 0 26.7 58.7 2.64 
720 26.0 29.9 
Marble 0 13.9 180.0 2.71 
720 12.2! 15.0 
Limestone 0 24.7 29.2 2.69 
720 24.0 23.0 
Dolomite 0 11.9 37.1 2.82 
720 11.9 14.2 
Soft steel 0 5.99 7.79 


1 Value at 600 kg./cm.?. 
2 Mean for range of 0-200 kg./cm.?. 
3 Mean for range of 400-720 kg./cm.?. 


Numerous cracks carry a limonitic stain and along them the calcite has 
been locally recrystallized, with an increase in the size of the crystals. 
Summary.—The compressibilities of a number of rocks of geophysical 
interest have been measured at thirty degrees centigrade for pressures 
up to 840 kg./cm.”. The prediction of Adams and Williamson of a rapid 
rise in compressibility as the pressure is decreased has been verified. 
It is a remarkable effect in the case of covered rocks, and a very much 
smaller one for uncovered rocks. The variation of the compressibility 
with pressure is greater for the more porous rocks. Thus, in order of 
decreasing magnitude of this effect, we have marble, granite, limestone, 
Sudbury norite and Vinal Haven diabase. The compressibility of the 
uncovered rock is always found to be less than that of the same rock when 
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covered. At pressures up to 300 kg./cm.? the difference may be large, 
but it rapidly approaches zero at higher pressures. In many cases the 
difference is already unimportant at 840 kg./cm.?. Simple explanations 
of these pressure effects can be given, based upon the existence of many 
easily compressed cracks and cavities between mineral grains. The linear 
compressibility often shows a large variation with direction in the rock; 
it is naturally most striking in layered or bedded rocks. 

Acknowledgments.—These experiments, suggested by the Committee on 
Geophysics of Harvard University, were performed in the Research 
Laboratory of Physics, through the courtesy of the Department of Physics 
of the same institution. I wish to thank also Professor P. W. Bridgman 
and Professor R. A. Daly for their active interest and advice. 
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COMPARISON OF THE STATICALLY AND SEISMOLOGICALLY 
DETERMINED ELASTIC CONSTANTS OF ROCKS* 


By W. A. ZISMAN 
DEPARTMENT OF GEOLOGY AND GEOGRAPHY, HARVARD UNIVERSITY 


Communicated May 20, 1933 


The speeds, V and v, of longitudinal and transverse waves of sound 
in an infinite, homogeneous and isotropic medium obeying Hooke’s law 
are related to the other elastic constants in a simple manner. If Young’s 
modulus, the cubic compressibility, Poisson’s ratio and the density are 
E, 8, o and p, respectively, then it has been shown by the mathematical 
theory of elasticity that these six constants are related as follows: 














eee E(1 — a) 

ig Va + 6) (l — 20) () 
: * * pees 2) 
¢-S (3) 


In short, we can regard such a medium as having its elastic behavior 
completely described when we assign the density and any two of the 
constants FE, o, B, V and v. 

It is fundamental for the theory of the seismological analysis of the 
structure and composition of the earth’s interior that we investigate experi- 
mentally the validity of assuming that rocks do behave elastically as homo- 
geneous elastic media obeying Hooke’s law. 

There are a number of methods of studying the general elastic behavior 
of rocks. We can measure the elastic constants E and o with a static 
testing machine. This will disclose any deviations from Hooke’s law. 
If rock specimens of different orientations are cut from a massive block 
of freshly quarried rock, we can ascertain the amount of elastic anisotropy. 
With the static constants so obtained the compressibility can be calcu- 
lated from equation (3). Comparison of this computed value with the 
compressibility measured directly for the same solid by subjecting it to 
hydrostatic pressure furnishes a sensitive test of any departure from homo- 
geneity and isotropy. Again, we can measure V and v seismologically 
in the quarry where the rock specimens for the static tests were obtained, 
and from equations (1) to (3) compute EZ, ¢ and 8. These may be appro- 
priately called the dynamical elastic constants. A comparison of these 
with the static values is a manifest necessity, if the seismologist is to use 
the elasticity as a criterion of the nature of rocks below the earth’s surface. 
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With the support of the Harvard Committee on Geophysics, the author 
commenced in 1931 a laboratory investigation of the elastic properties 
of rocks, including some in which the seismic-wave velocities at and near 
the surface have been measured. 

The purpose of this paper is to interpret the related results of all these 
experimental studies. 

The laboratory measurements have already been published."? Leet 
and Ewing have published dynamical elastic constants for Quincy granite’® 
and a report on those for Rockport granite and Sudbury norite is in 
preparation. 

The speeds observed at Quincy were V = 4.96 km./sec. and v = 2.48 
km./sec. From these, the dynamical constants were found to be 


E = 4.3 X 10" dynes/cm.? 
o = 0.33 
B = 22.7 X 10-" cm.?/dyne 


with probable errors of about 2 per cent. 

The statical measurements of E and o and the linear compressibility 
measurements on the granites from Quincy, Westerly and Rockport all 
demonstrated their anisotropy. This became particularly evident in the 
case of the Quincy granite of which many samples were studied. In all 
instances the values of E for even almost perfectly fresh surface rock were 
found to be smaller than those for fresh rock taken from depths of 100 and 
235 feet. The measurements of o show little difference with samples or 
locality. The mean values for E and o in Quincy were 3.5 X 10"! dynes/ 
cm.” and 0.10, respectively. Thus the statical values are definitely smaller 
than the dynamical. 

In the author’s report on the measurement of the hydrostatic com- 
pressibilities of rocks, two values of the compressibility were distinguished : 
the compressibility of the covered rock and that of the uncovered rock. 
The former was found by studying the change in volume of the rock when 
covered by an impermeable cover of metal foil. The latter was found 
by studying the changes in volume of the rock when exposed directly 
to the hydrostatic fluid so that the liquid was able to penetrate into the 
pores of the rock. The former value always exceeds the latter. The 
difference between the two values of compressibility is least for the com- 
pact rocks and greatest for the most porous. As the hydrostatic pressure 
increases, the two values of the compressibility approach each other, the 
difference usually becoming small at pressures around 800 kg./cm.’. 

This contrasted behavior is due to the existence of cracks and cavities 
among the crystals making up the rocks. 

The mean value of the initial compressibility of covered Quincy granite 
is 90 X 10-' cm.?/dyne. The mean initial value for the uncovered 
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granite is 21.6 X 10~-'* cm.?/dyne. Evidently, the compressibility of the 
uncovered rock agrees well with the seismologically determined value. 

Summarizing the results for granite, we find the statical constants E 
and o are considerably smaller than the dynamical; the statical com- 
pressibility of the covered rock is much greater than the dynamical value; 
and the statical compressibility of the uncovered rock is in good agreement 
with the dynamical value. 

In the norite at Sudbury, Ontario, Leet obtained accurate measure- 
ments on elastic wave speeds. The results were V = 6.22 km./sec. and v 
= 3.49km./sec. The density being 2.85, these give: 


E = 8.82 X 10" dynes/cm.? 
og = 0.27 
B = 15.6 X 10- cm.?/dyne 


The mean statical values are: 


E = 8.36 X 10" dynes/cm.? 
a = 0.22 


The initial compressibility of the covered rock is 28.6 X 10~!* cm.?/dyne, 
and of the uncovered rock 16.4 X 107" cm.?/dyne. Thus the statically 
determined constants E and o in norite, though in better agreement 
with the dynamical values than was the case in granite, are again defi- 
nitely smaller than the dynamical. The compressibility of uncovered 
rock is in good agreement with the dynamically determined value. 

The Sudbury norite is more compact than Quincy granite; this is 
indicated by the much higher value of Young’s modulus for norite. 

These two comparisons of field and laboratory work suggest the following 
generalizations: The statically determined elastic constants E and oa 
are smaller than the respective seismological or dynamical values; the 
compressibility of uncovered rock is essentially the same as the seismically 
effective compressibility; the more compact the rock the more nearly 
do the statical and dynamical constants agree, and the nearer together 
are the respective compressibilities of covered and uncovered specimens. 

Recently Richards published the results of his measurements on the 
statical and seismological elastic constants of limestone cores cut from 
various depths in a large oil-bearing structure in southwest Persia.‘ As 
he was mainly concerned with seismic waves traveling in directions 
parallel to the bedding of the rock, he confined his statical measurements 
to those giving the elastic constants E and o in the plane of the bedding. 
The mean statical values of E and o were 5.7 X 10'! dynes/cm.? and 0.25, 
while the mean seismological values were 5.5 X 10"! dynes/cm.? and 
0.29, respectively. The average compressibility computed from his seismo- 
logical measurements is 22.8 X 10~'* cm.?/dyne. 
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The Persian limestone is closely similar to a compact Pennsylvania 
limestone studied at Harvard, for which, when covered, the observed 
compressibility is 29.0 X 10~'* cm.?/dyne and, when uncovered, 24.7 X 
10-13 cm.?/dyne. The small observed difference between the two types 
of compressibility clearly indicates great compactness of the rock. The 
measurements on the linear compressibility for directions perpendicular 
and parallel to the plane of bedding of the Pennsylvania limestone proved 
strong anisotropy. This was shown also by measurements of E and o on 
specimens cut in and at right angles to the plane of the bedding. The 
mean values of £ and o in the plane of bedding were 6.7 X 101! dynes/cm.’. 
and 0.25, respectively. Thus my statical results are in fair agreement 
with Richards’ statical values. 

As in the comparisons with granite and norite, we find bedded limestone 
exhibiting a lower statical value of o than that obtained seismologically, 
and here too the compressibility of the rock, when uncovered, is nearly 
that deduced from the seismic waves. 

The departure of the rocks from ideal elasticity may be tested in another 
way. Table 1 gives the results of computing from the values of E and o 
and from equation (3) the compressibilities at two different pressures, on 
the arbitrary assumption that all the rocks behave as homogeneous, 
isotropic media. In nearly all cases the computed compressibility is 
smaller than the observed compressibility, the difference increasing with 
the porosity of the material. The marble is the most porous; the Vinal 
Haven diabase the least. 


TABLE 1 
COMPRESSIBILITY AT COMPRESSIBILITY AT 
1.2 kc./cm.? 56.0 KG./cm.? 
CALCULATED OBSERVED CALCULATED OBSERVED 
ROCK CM.?2/DYNE cM.?/DYNE cM.?2/DYNE cM.?2/DYNE 
Limestone, Pennsylvania 32.6 X 10718 29.0 XK 10—-" 36.4 XK 107-!3 28.3 XK 10-"8 
Diabase, Vinal Haven 14.5 17.0 14.3 16.5 
Norite, French Creek 33.2 58.0 30.7 48.0 
Granite, Rockport 100 ft. 66.2 84.0 54.5 65.0 
Granite, Quincy 235 feet 70.6 86.0 47.5 67.0 
Marble, Vermont 71.0 140.0 41.4 59.0 
Quartzitic sandstone, Penn- 
sylvania 36.9 56.4 35.7 48.7 
COMPRESSIBILITY AT 
36.5 Kkc./cm.? 
CALCULATED OBSERVED 
cM.?/DYNE cM.2/DYNE 
Norite, Sudbury 20:2 KX 10-*'23:.6 < 10-* 


However, there is no theoretical justification for assuming anything 
like equivalence for the compressibilities, computed and observed, at 
such a pressure as 11.2 kg./cm.?. Possibly if the computed value for this 
pressure be compared with the observed hydrostatic compressibility at 
three times as much, or ‘33.6 kg./cm.’, a more significant result might be 
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obtained. Even then the respective differences are considerable and 
greater than can be explained by experimental error. It is further note- 
worthy that the differences are not due to lack of isctropy in the rocks, 
Vinal Haven diabase and Quincy granite at 235 feet, for example, being 
almost ideally isotropic. 

The analogous values of Poisson’s ratio may also be compared. Again 
assuming the validity of equation (3), and also assuming the correctness 
of the observed initial values for compressibility and Young’s modulus, 
we compute o for the rocks in both the covered and uncovered conditions, 
and at zero pressure. The results are shown in table 2. The values given 
in the right hand column are tending toward those expected if the rocks 
were free from cracks and pores. Although these were by no means all 
closed by the actual pressures, the computed values of o are not wildly 
unreasonable. In utter contrast are the values given in the first column 
of table 2. Manifestly the classic equations cannot be applied to most 
rocks of the earth’s crust except when these are at such depths as to ensure 
the closing of cracks and pores. 


TABLE 2 
POISSON’S RATIO POISSON’S RATIO 
ROCK USING B, COVERED USING 8, UNCOVERED 
Granite, Quincy 100 ft. —0.11 +0.31 
Granite, Rockport 100 ft. —0.06 +0.38 
Granite, Quincy 235 ft. +0.03 +0.39 
Diabase, Vinal Haven +0.21 +0. 25 
Norite, French Creek —0.14 +0.35 
Norite, Sudbury +0.06 +0.27 
Quartzitic sandstone, Pennsylvania —0.12 +0.22 
Limestone, Pennsylvania +0.24 +0.28 
Marble, Vermont —0.36 +0.43 


The probable cause of the difference between the statical and dynamical 
compressibilities of standard rocks at low pressure may now be considered. 

Let us assume that the dynamical methods, of which the seismological 
method is but one, give values of E, o and 6 for the rock when free from 
cracks and easily compressed cavities. The statical values of E and o are 
lower because part of the yielding of the rock is increased by deformation 
of cracks and cavities. The compressibility of the covered rock is high 
for the same reason. The compressibility of the uncovered rock is the same 
as the dynamically obtained compressibility, because the uncovered rock 
is soaked with the fluid transmitting the hydrostatic pressure and hence 
the measured effect is only that due to the compressibility of the solid 
material in the rock. 

It is desirable, if possible, to give a physical foundation for the generali- 
zation made above, which, after all, is an inference based on a limited 
number of comparisons of seismological and statical measurements. 
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A qualitative explanation can be based upon considerations of wave 
energy. A plane wave of sound is propagated in the direction indicated 
by the arrows in figure 1. The wave is supposed to pass over two cracks 
A and B, such as might mark the boundary between adjacent crystals 
in this rock. The portion of the wave striking A and B would suffer 
loss of energy, due to the two reflections at the rock-air and air-rock 
boundaries. A small amount of wave energy would be transmitted across 
the boundary and a large amount would be scattered around the cavities 
A and B. As the cavities in the rock will not in general form an orderly 
array, we realize that such scattered energy will be of small importance 
except as a cause of absorption for the wave propagation over long dis- 
tances in the rocks. Thus we are led to conclude that most of the wave 
energy transmitted through the rock masses will be made up of the beams 
transmitted between the crevices and cavities. But these beams will 
move with the velocity 
corresponding to the com- 
pact rock material. In other 
words the seismically effec- Pe CE g 
tive elastic constants will cavity cavity 
be the high E and o one 
would obtain in a statical 
measurement on a rock of { | t 
that mineralogical composi- 
tion, free from crevices and 
cavities. For waves of low 
frequency the elastic con- 
stants will be the statical . FIG. 1 
values, and for high fre- 
quencies they will be the high values due to a rock free from cracks 
and cavities. 

Recently, Dr. D. S. Muzzey, working in the Research Laboratory of 
Physics at Harvard University, has obtained dynamical measurements 
of E for various rocks by using a magnetostriction method to set the 
specimens vibrating. The same five-eighths-inch rock cylinders were used 
that were employed for my measurements of the hydrostatic compressi- 
bilities. Thanks to the kindness of Dr. Muzzey, it is possible to present 
here some of his results. A detailed description of method and results 
will be published separately by him. The following data in table 3 are for 
sound frequencies ranging from 9000 to 13,000 cycles per second. 

The experimental error is estimated to be plus or minus two per cent. 
Dr. Muzzey encountered much difficulty with the less compact rocks like 
granite and marble, due to absorption. This is just as might have been 
predicted. The mean of his values of E is 9.1 < 10"! for Sudbury norite. 
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This compares well with the seismological value of 8.82 X 10!! dynes/cm.?. 
Again, the very compact rock, Vinal Haven diabase, shows good agreement 
between the statical value of 10.2 X 10" dynes/cm.? and Dr. Muzzey’s 
dynamical value of 9.85  10!! dynes/cm.?. 


TABLE 3 
E (MUZZEY) E (STATICAL) 
ROCK DYNES/CM.? DYNES/CM.? 

Norite, Sudbury 1 9.13 X 10}! 8.07 X 10" 
Norite, Sudbury 2 8.85 7.96 
Norite, Sudbury 3 9.00 8.05 
Norite, Sudbury 4 9.48 9.36 
Diabase, Vinal Haven 9.85 10.2 
Obsidian, Lipari 6.82 


The discussion shows definitely that the comparison of field and statical 
laboratory measurements on elastic properties of rocks near the surface 
of the earth is a difficult matter for most rocks, because of the presence 
of numerous small cracks and cavities. The lack of agreement between 
dynamically and statically determined elastic constants is a real property 
of the rock, and is smaller the more compact the rock. 

The cracks and cavities in rocks will doubtless cause dispersion and 
absorption for elastic waves. If measurements of the elastic constants E 
and o are desired for comparison with data derived from seismological 
observations, they can be obtained by either statical or dynamical labo- 
ratory methods on highly compact rocks, and by dynamical methods 
alone on porous or loosely compacted rocks. The laboratory measurements 
of the compressibilities of the uncovered rocks do agree closely in every 
case with the seismologically measured values. Thus such laboratory 
measurements of the compressibility should prove of value to seismologists. 

For compact surface rocks, and for rocks at depths in the earth where 
the pressure due to the weight above is such as to ensure freedom from the 
presence of cracks and cavities, it is certainly correct to treat them as 
homogeneous elastic media obeying Hooke’s law and possessing the same 
statical and dynamical elastic constants. 

* Paper No. 5 published under the auspices of the Committee on Geophysical Re- 
search and of the Division of Geological Sciences at Harvard University. 
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THE GENERAL CASE OF NON-HOMOGENEOUS LINEAR 
DIFFERENTIAL EQUATIONS 


By W. J. TRJITZINSKY 
DEPARTMENT OF MATHEMATICS, NORTHWESTERN UNIVERSITY 


Communicated June 12, 1933 


1. Our present object is to establish the asymptotic nature of solutions 
of the equation L,(y) = b(x), where the differential operator L,(y) is of 
the form 


L,(y) = ao(x)y™ (x) + ar(x)y~? (@) +... + ea 9 (@) + 
An (x) y(x) [ao(x) = 0; On (x) # 0) (1) 


The coefficients a;(x)(i = 0, 1, ... m) and b(x) are analytic for | x | = po 
(|x| + ©) being representable by convergent series of the form 


B u—1 


1 im 
a,x? +a,-,%7? +... +an? +at+a,x ?+..., (2) 


where p (2 1) is an integer. More generally, these coefficients may be 
considered to be merely asymptotic, in certain regions, to series of the form 

’ (2).1 It is assumed that not all the coefficients of the series (2), associated 
with a(x) and a,(x), are zero. Thus, without any loss of generality, we 
take ao(x) = 1. No restrictions are made concerning the roots of the charac- 
teristic equation corresponding to the homogeneous problem L,(y) = 0. 

2. At the basis of the present treatment lies the asymptotic theory of 
the homogeneous problem, which the present author has completely de- 
veloped (with no restrictions on the roots of the characteristic equation) 
in a detailed paper to appear elsewhere. In the sequel this paper will be 
referred to as (JT). The development of the theory of the homogeneous 
equation (or linear system, as well), in (7), has been based on successive 
applications of a method of “‘Iterations’’ (inspired by a method so called 
and originally applied by G. D. Birkhoff to difference equations), fac- 
torization and suitable integrations. 

Those results of (J) which are necessary for the purposes at hand will 
be now stated. 

Let L(y) denote L,,(y) with the coefficients replaced by the correspond- 
ing series. As follows from a work of E. Fabry (Thése, 1885, Paris) the 
formal equation L*(y) = 0 is satisfied by a full set of formal series solutions 
of the form 


v=l—1 l-—v 


5; (x) = €%) xr d; (x), Q(x) = LD Gx*, 
v=0 


0;(x) = OP(x) + O}(x) logx +... + O7*(x) log (x), (8a) 
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ee ' 
OF (x) = OF + Oly & + OP A+... (3b) 
(u = 0,1, ... mm; 4 = 1,2, ... 2), 


where /;, m;, k; are integers [m; 2 0; k; = r;p; integer r; 2 1].? 
A series O(x) of the form (3a), (30) is called a 0-series. 
Branches, extending to infinity, along which 


RQ ;(x) = 0 [Q.5(%) = Qxr) — Q)] (4) 


for some 7 and j (i ~ j) will be termed Q curves.’ Let k be a suitable in- 
teger. The complete vicinity of infinity 


% S60 =argx S 2rk + H (5) 
(|x| =p =o >0) | 
is divided in a number of successively adjacent regions 
Ba, Ry... BR (5a) 


separated by Q curves; interior neither of these regions are there any Q curves. 
The region R; will have a certain Q curve in common with R,. Within 
any particular region R;.a certain ordering of the values RQ;(x) is main- 
tained. 

Let R; denote any particular one of the regions (5a), whose boundaries 
have different limiting directions at infinity. Let 


Oinig(®) = Oig(x) — Qj,(e) [# 0; H=1,...2] (6) 


constitute the totality of differences Q; ;(x), each of which has a non- 
negative real part in R; and is such that for some 8 


a} RO: 4.6 
|x Ble Cnn” —» 0 (6a) 


along both boundaries of R;. Suppose there are some such functions (6). 
Let R; denote then a subregion of R; with one boundary in common with R; 
and another boundary, interior R;, such that along it all the left members of 
(6a) increase indefinitely for every 8B (>0). We let R; denote an analogous 
subregion of R;, with the other boundary in common with R;. The limiting 
directions (at infinity) of the two boundaries of R; are correspondingly 
the same as those of R;; the same is true of R! (provided the regions are 
suitably specified, as is possible). 

The main result of (TJ) is embodied in the theorem. 

THE FUNDAMENTAL EXISTENCE THEOREM. Consider an equation 
L(y) = 0 of the type specified above. Let R; be any particular region of the 
set of regions (5a). If there exist no functions (6) such that (6a) holds, or if 
the boundaries of R; have the same limiting directions (at infinity), there exists 
a full set of analytic solutions y,(x) (v = 1,2, ... ) such that 
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Yo (x) sips €2rl#) 700, (x) sa 5y(x) (7) 
(y= 1, 2,.... # 2 i &,). 


The asymptotic relations here are with respect to the power series factors. 
If there exist functions (6) such that (6a) holds, while the boundaries of R; 
have different limiting directions, consider the two subregions R; and R! 
of R; (as specified in the italics following (6a)). In this case, also, there will 
exist a full set of solutions satisfying (7) in Rj, and another set of solutions 
satisfying (7) in Rj. The s;(x) are formal solutions of the formal equa- 
tion Lx (y) = 0. 

Analogous results hold for linear differential systems. 

3. Consider now the non-homogeneous problem. Let R denote any 
particular region for which existence of solutions (of the equation L,(y) = 
0), satisfying the asymptotic relations (7), had been asserted in the Funda- 
mental Existence Theorem. Thus, R may denote R; (a region of the set 
(5a)) or it may denote one of the regions R;, R;. This region R will be 
subdivided further as follows. Consider the branches, extending to in- 
finity, along which 


RQi(x) = 0 (8) 


for some i. In general, of course, there are several branches satisfying a 
particular equation (8). R will consist of a number of consecutively ad- 
jacent regions 


RR? G=et&:..D (9) 


separated by branches each satisfying an equation (8); moreover, no such 
curve is to be interior a region R®. It may happen that there is only one 
such region R®, necessarily coincident with R. 

Let R® be a particular region of the set (9) with the limiting directions 
(at infinity) of its boundaries distinct. Let 


—Qin(x) [#0; H=1,2,...g) (10) 


denote the totality of functions —Q;(x), each of which has a positive real 
part interior R and is such that for some 6 


|x~* |e *@in™ —» 9 (10a) 


along both boundaries of R®. It may occur that, given a particular re- 
gion R“, there exist no corresponding functions (10) such that (10a) 
holds. Suppose there are such functions (10). We let then R™? denote a 
subregion of R® with one boundary in common with R® and another 
boundary, interior R®, such that along it the left members in (10a) (H = 
1, 2, ... q) increase indefinitely for every 8B (> 0). On the other hand, R” 
will denote a similar subregion of R, with the other boundary in common 
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with R®. These regions can be so constructed that the limiting directions 
of their boundaries are correspondingly the same as those of R“. 

Let y;(x) [= y;, (x); i = 1, 2, ... nm] constitute a fundamental set of 
solutions of the equation L,(y) = 0, satisfying in R the asymptotic rela- 
tions (7) of the fundamental Existence Theorem. Let (y'4-?(«)) denote 
a matrix with yF—» (x) in the i-th row and j-th column (7,7 = 1, 2,... m). 
A solution of L,,(y) = b(x) would be given by 


v(x) = 3 mile) SC) aaledde, (11) 


=} 
where the 9; ;(x) are elements of the inverse of the matrix (y{4-  (x)); that is, 
Gis(x)) = Of"). 

The integrations in (11) are specified, for a region R” or R’ (or R™"), 
as the case may be, according to the integration methods used in (T). Ap- 
plication of these methods is possible in view of the conditions satisfied by 
the regions R®, R®’, R®”. 

On taking account of (7) and of the asymptotic form of 6(x) it will 
follow that, for x in R® or R®’ (or R®”), 


rf 
y(x) ~ x* (x), (12) 


where / and k, (> 0) are integers and 0(x) is a O-series. 

EXISTENCE THEOREM. Let L,(y) = b(x) be an equation of the type 
specified in (§1). Let R® be any particular region of the set of regions (9). 
If there exist no functions (10) such that (10a) holds, or if the boundaries of 
R® have the same limiting directions at infinity, there exists an analytic 
solution y(x) satisfying the asymptotic relation (12) in R. If there exist 
functions (10) such that (10a) holds, while the boundaries of R have differ- 
ent limiting directions, we consider two subregions R’ and R®" of R® (as 
specified in the italics following (10)). In this case, also, there exists an ana- 
lytic solution y(x) satisfying (12) in R“’, and another solution satisfying 
(12) in R". The right member in (12) constitutes a formal solution of the 
formal non-homogeneous equation corresponding to the equation L,(y) = 
b(x). 


1 Correspondingly, the results will be asserted for suitably restricted regions. 
2 Except with x and log x, superscripts here do not denote powers. 
3 RQ(x) denotes ‘‘real part of Q(x).” 
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A CANONICAL FORM FOR A SET OF VECTORS 
By M. S. KNEBELMAN 


DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 


Communicated April 22, 1933 


By a canonical form for a set of vectors we understand the form of the 
components of these vectors in a particular codérdinate system in which 
these components have some desired simplicity. In this note we propose 
to obtain a canonical form for a set of contravariant and for a set of co- 
variant vectors. Without loss of generality we may assume that the given 
set is independent and we shall also assume that the set is complete since 
the reduction of an incomplete set is accomplished by reducing the com- 
pleted set. 

1. Let &)(x) be a set of r contravariant vectors in an n-dimensional 


space, r < n. The vectors being independent implies that the matrix) , 


|| E% || is of rank 7. We shall suppose that the functions £(x) are analytic 
in the neighborhood of the origin. The set being complete, we have by 
definition 








J ° G 3 @ / a 
20 ne — & a = Ciy'x)EGy, (1.1) 


where Ci (x) are the structure {unctions of the set. The rank of the matrix 
\| € G || being 7, a transformation of coérdinates, which amounts to a mere 
renumbering, will insure the fact that | f;) | ~0,i,j=1,...,r. Then the 
system of differential equations 
de 
a, <P = 0 (1.2) 
fii) Sx" 
admits (because of (1.1)) ™ — r independent solutions ¢/,), o = r + 1, 
.., m. When we apply to our set of vectors a transformation of co- 


ordinates defined by x” = ¢@,), x' = x’, in the new codrdinate system! we 
shall have 


§G;) = 0 o am Fis (1.3) 
Moreover, any transformation of the form 2” = x(x! ... x") x’ = 


x’*(x’ +1, ..., x") will preserve the property (1.3). We have therefore 
reduced our set of vectors to a form in which all components after the r" 


are zero and such that 


re) ra) 
dk — tp) i Cie Eins (1.4) 


bin SF ox" 
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where the codrdinates x” +}, ..., x"—if they enter at all—may be treated 


as parameters. A transformation of the codrdinates x, ..., x” can always 
be found such that & = 6%). In this coérdinate system (1.4) becomes 


Eu 





a 2 
ag Cre En)» (1.5) 
Ebi) dij.) 
Ei) SF Fa g tin) Sa ia Cj. ki tiny, fu ky a, Ue 
From the first of these it follows that 
Ey (x! sig = £, (0, x2, ..., x")-Ab,(x! ... x"), (1.6) 
where A}, is a non-singular matrix defined by 
x! ¢1)3 
Aln = 8% + Cie (0, x? =)21+ Ch, y + Chass + ‘v9 


a subscript preceded by a comma indicating intrinsic covariant differentia- 
tion, i.e., if F’ is a set of scalar functions then 





Fi = &y a + PGi. (1.8) 
It can be easily verified that 
a“ = Ch Als (1.9) 
and by means of this to establish that, if A, is the inverse matrix of A, then 
= = Ch Alm (1.10) 
At, = of + Ci(0, x*, ..., 2") x1 + Ch-y = +... (1.12) 
where 
Fi = &j a + F*Cj,. (1.12) 


When (1.6) are substituted into (1.5) we get 


Oi, dé; i i 
thy iw Ea) res * | chahet, — tntio 





OAY 0Aq 
(at, — Ah > LY Lai, Ht. (1.18) 
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where ey are functions of x*, ..., x” at most. Since the left-hand side of 
(1.13) is free of x', so is the right-hand side; therefore its value is unaltered 
if we put x! = 0. Now from (1.7) and (1.11) we have A‘, = A‘, = a 
and 0A4,/dx’ = C4,(0, x%, ...) 5; for x1 = 0. Thus equations (1.13) 
become 


1 kb, 1 dk; i 
Ei ne 4 — eh, a = (Cha + Chitty) + Chieti» (1.14) 





where all functions involved are the values of the corresponding original 
functions for x! = 0. 

If we consider (1.14) for i > 1 we shall have (because ty = 0if7> 1) 
a complete set of r — 1 independent vectors since | e | = | eH, | ¥ 0, 
whose structure functions? are 


Crs (x?, ...,%") = Che, x2, ..., 0") + Chet + Chet. (1.18) 


We may now apply a transformation of coérdinates of the form x! = x}, 


x = x" (x2, ..., x"), x” = x”, such that #,, = dé, = é,) Hy) = af 
and as before we shall find from the equations of completeness that 
Ey = En (0, 0, x* ...) AD}, (@? ... x”), (1.16) 
where 
(x?)? 


Ad}, = OR + Cofi(0, x8, ...,2")x? + Coit y2 "dame 


Obviously this process may be repeated r — 1 times, since for one vector 
there are no equations for completeness. We now observe that #1) = 
i and because of (1.6) and (1.16) the second vector is such that for x! = 0 
its components are £2) (x? ... x”), 657 > 1. Similarly the third vector is 
such that for x1 = x? = 0 its components are £{g) (x, ...), £3) (x3, ...), 
53 i > 2 and so on. We may therefore combine these results into the 
following theorem: 


A set of s (S n) independent contravariant vectors in an n-dimensional 
space being given such that the completed set contains r (s< r S n) independent 
vectors, a codrdinate system may be found in which all components after the 
r” are zero throughout the space, the remaining components being 


t, : 64 throughout the space 
Ey > Ela)» i> 1forx! =0 
£(3) : Ely), a), ay = 2 for xi = x2 = 0 


E(s) ° Ey), PAS Te ee —_ 1 for x} =giz eines = (, 
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It may be shown further that the transformations involved in the above 
process have enough arbitrariness in them to enable one to choose the sub- 
space V,_;4, arbitrarily—except for certain restrictions—on which the s 
vectors are to have the above canonical form. 

In particular an ennuple of vectors and a curve being given, the curve 
not belonging to the congruence determined by the ennuple, a coérdinate 
system exists in which the curve has the equations x! = x? = ... = x""!= 
0, and the vectors have components such that £%) = 6g fora 2 8B. 

2. A set of r (< m) independent covariant vectors may be reduced to 
canonical form by a similar process although there are some essential differ- 
ences. Let nia) (x) be a set of independent vectors and without loss of gen- 
erality we may assume that | n{” |~ 0. The set being complete, we have 
by definition 


(i) Ong” 
ox? Ox™ 





= Cie(x) ne'np”, (2.1) 


the C’s being the structure functions of the set. It can be shown that the 
system of partial differential equations n{?0x*/dx” = 0, ¢ > r is com- 
pletely integrable* and that a codérdinate system therefore exists in which 
n? =00>r. Moreover, a transformation of the form 


x’ = x(x}, ..., x’) <r 
x =x" (x!,...,x") o>r 
will preserve this characteristic of each covariant vector. In the co- 
ordinate system determined above (2.1), being tensor equations, give for 
asrp>r 
(i) 
J 


Ox? 








= Cin} nf? (2.2) 


oo 
ox* 


Ox! 

From the first of these it follows that the components are free of gi th 

.., x" and from the second—since | nf” | ~ 0—that Ci, are free of the 
same variables. Hence 

A complete set of r (< n) independent covariant vectors being given, there 
exists a coérdinate system in which the components n\? = 0 ¢ > r and the 
remaining components are functions of x}, ..., x” at most. 

It may be observed that this theorem does not hold, in general, for a set 
of contravariant vectors. 

We now look for a coérdinate system in which the same single component 
of every vector except one shall be zero. Since | n{” | ¥ 0 not every 
(r — 1)-rowed minor is zero and without loss of generality we assume it to 
be the cofactor of n{”. Then the equations 
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Ox! 
ox! 


(i) 


nj = § 


admit a solution. Hence there exists a codrdinate system in which nf? = 
4; and when this is put into (2.2) for k = 1 we get 


anf! 





ox! = Ci ni”. (2.3) 
Hence, as in $1, 
nj(x!, ..., x") = nf(0, x2, ..., 2”) Bi,(x, ..., x’), (2.4) 
where 
(xt)? 
Bu, = 5 + Cin (0, x?, ..., xx! + Cissy — T +... GS 


Ciaz1 = —Ciaa (cf. (1.12)) 
=z = Cin Bin. (2.6) 


When (2.4) are substituted into (2.2) we get 





(m) on™ : aes OBi, OBi, 
ee Bm, (a) () m{ .(h a 
ax" ing CiyBioBi Brin} ne? — By; (11 3 rw nh Sx 2B), (2.7) 
where 7{” = nf” (0, x?, ..., x”) and BY; Bi,,=6. When i or j=1 (2.7) 
are identically satisfied and since the left side of (2.7) is free of x! so is the 
right. Setting x' = 0 we find 


(m) fo") , 
Reais SES ‘e r)(h) (l) » 
Dak a sa Cii(0, 2 i we x”) nj, Nk, Jis ky a ee (2.8) 
We observe that »{” is a set of r covariant vectors in a space of r — 1 
dimensions. Also since |n{”| # 0, the first vector is linearly dependent 
on the last y — 1 and therefore 


( ‘ 
ma? + on i = 0, (2.9) 


the y’s being determined by these equations. Hence (2.8) for m > 1 re- 
duce to 


(1) ) 
in One , 


a Ox! ss Chit ni nk, (2.10) 








where 
Cit, ris Cin (0, x’, sees x") + CE Phi + Con Fl: (2.11) 
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We are thus confronted with the same problem we had at the beginning of 
this section except that the number of vectors and of variables is diminished 
by one. As before we can find a coérdinate system in which (2.12) 


nf = 8 (2.12) 
and because of (2.10) 
a = nf (0, x? ... x”) Bi, (x? ... x); (2.13) 


where Bi,, are defined analogously to Bi, (cf. (2.5)). Hence nf = 6) 
throughout the space and for x1 = 0 nf? = of (x2, ..., x”) 7 = 6 
(4; > 1). Similarly the third components of the vectors are such that for 
x! = x? = 0 they reduce to n{? (x? ...), nf («* ...) and 62 (i, > 2). 
Hence 

A complete set of r (< n) independent covariant vectors n (x1, ..., x") 
being given, a codrdinate system can be found in which n® = 0 for o > 1, the 
first components are 6', throughout the space, the second reduce to n\ (x*, ..., 
x’), 62 i; > 1 for x! = 0, the third components reduce to n{?, n®, 63 ig > 2 
for x! = x? = 0, etc., the r™ components reducing to n”, (x"), ... 1,7} (xn), 
1 on the subspace V,_,+, defined byx! =x? =... =x7~'=0. 


1 That this transformation is non-singular is easily proved. 

2 The structure functions of a complete set of vectors besides being skew-symmetric 
satisfy Jacobi’s identities which may be written as C};, + Ch; + Chy = 0. It is 
easily verifiable that C’}',, satisfy identities of the same sort. 

3 This system is of a type similar to that studied by K6énig and its general solution 
depends on n arbitrary functions of m — r variables thus insuring the non-singularity of 
the transformation of codrdinates. 
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ON A THEOREM OF FEJER 
By W. C. RANDELS 


DEPARTMENT OF MATHEMATICS, BROWN UNIVERSITY 


Communicated May 25, 1933 


L. Fejér has shown! that, if a function f(x) is odd, and positive and convex 
upward in the interval (0, 7), the (C, 3) transforms of the partial sums of 
its Fourier series will be convex upward in the same interval, and, if the 
function is symmetric about the point x = 1/2, the (C, 2) transforms will 
also be convex in that interval. In this paper we propose to show that in 
both cases the indices obtained by Fejér are the lowest possible. 

Let us define the function 








x 
i O53 236 
f(x) = f(—x) = — f(x). 
"=— Xx 
) gee Oe Sls a 
xr—b 
Then 
= 2 ; 
= ——_——— sin nb sin nx. 
fe) = 2 ize-ie 
The nth (C, a) transform of the Fourier series of f(x) is given by 
Sn (x) 
n + a\> 
b(r — b) ( i ) 
where 


St (x) = > - rt a) sinh ain 


=1 a v 


The factor multiplying S% is positive so we may disregard it and determine 
the condition that S; be convex upward. This condition requires that 
D? S% (x) must be < 0. But D?.S% (0) = 0, and therefore we must have 
D3 Si (0) < 0. Differentiating we get 


f=-pDs0 => (7520 %) v sino. 


=] no F 
We now investigate 
& = (a+ 1)sinbd + 2sin 2b = (a + 1) sind + 4sin bcos bd. 


If this is to be non-negative we must have a 2 —4cosb — 1. Now let 
b —~> 7, then the expression on the right tends to 3 and hence we must 
have a 2 3. 
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If b = /2, f(x) is clearly symmetric about the point x = 7/2. It is 
readily found that 


5 
py . ee See ae es n+1 
i 2) ysiny 5 (—1) | 9 : 
Also ; 


bwls 


z n 
h= De =1, 9 +10 -; 


ess 


’ 


according as m = 1, 2, 3, 0 (mod 4). 
If B > a, it is well known that 
S@) => Poa 
Jas n—yv 
and 


Ge6e SS ee 
(n — v)! 





ef 


A linear combination, with positive coefficients, of functions which are 
convex upward is itself convex upward, hence if Sy (x) is convex upward 
for all n, S§ (x) will also be convex upward for all m. Therefore, since 
S} (x) is clearly not convex upward for all n, we must have a > 1. We 
will now suppose 1 < a < 2. Then 


f= a 
v=]1 
Let n = 4m, then 


a ee 


v=] ee »=1 


where the first sum is taken over integers of the form vy = 1 (mod 4) and 
the second over integers of the form v = 0 (mod 4). But 


(Corte te 2) _ (sorte 2) _ Gaeta). ey) 

n—v+2 n—v (n — v + 2)! 

a oe in SA PO, 7 Oe Oe ee 
(n — v)! 


y (n—v+a-— 2)...(a— 1) 
—(n—v+2)(n—v+1)} ae <0 











and the last term in the second sum is 


: 2 
4(@-2)(a@+1) <0. 
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Also 
~,{n-vta-2 "_{n-—-v+a-2 
Py | a= 8 VE i= P ) 
= (F271) 2 ory, 


This shows that /{,, will be negative for sufficiently large values of m, if 
a<2. 


1In the address delivered at the Chicago meeting of the American Mathematical 
Society, which is forthcoming in the Bull. Amer. Math. Soc. 


ARISTOGENESIS, THE OBSERVED ORDER OF 
BIOMECHANICAL EVOLUTION* 


By HENRY FAIRFIELD OSBORN 
THE AMERICAN MusEUM OF NATURAL HISTORY 


Read before the Academy, Tuesday, April 25, 1933 


The term “‘biomechanical’’! applies to organic as distinguished from 
inorganic processes;? for example, the evolution of the living invertebrates 
and vertebrates of the past as observed in paleontology is chiefly recorded 
in the hard parts of animals formerly full of life. For the observed order 
of biomechanical adaptation and evolution, which is the subject of the 
present communication, no fitter, no more apt term can be found than 
aristogenesis proposed but not defined in my contribution to the “Centenary 
of Evolution” symposium of the British Association for the Advancement 
of Science in 1931 and in my last contribution to the National Academy.* 

First, the term aristogenesis is derived from the Attic Greek “‘aristos’’ 
(apioros) signifying “‘best in its kind,” in Euripides referring to warriors, 
the ‘‘best,” the “‘most useful; in Homer referring to chariots, the ‘‘best’’ 
or ‘‘finest.”4 The Greek “‘genesis’’ (yeveois) from the Aryan base meaning 
“to produce, bring into being’’ means “‘origin, source, birth, descent; ...mode 
of formation or production;’’® also, in Plato, ‘‘an origin, source, productive 
cause; origination, making (even of common things).’’® 

Second, the term aristogenesis relates to an actually observed order of 
nature, a concrete process; it has no theoretic significance such as the 
“entelechy” of Aristotle, the ‘‘pangenesis’” of Darwin, the “‘vitalism’”’ of 

* This is the ninth contribution on the Origin of Species and the principles of bio- 
mechanical evolution as demonstrated in paleontology. The numerical references in 


the present communication are from the complete bibliography in Osborn’s “‘Fifty-Two 
Years of Research” which will be supplied on application to the Osborn Library. 
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Driesch, the ‘“‘bathmism”’ of Cope, the ““emergence’’ of Morgan, the ‘‘holism”’ 
of Smuts. It is rather to be compared with the “biogenesis” of Haeckel, 
a term signifying the recapitulation of ancestral evolution in embryonic 
and foetal development. 

Third, aristogenesis as a process of nature is not to be challenged as a 
concept, still less as an explanation, hypothesis or theory. At present 
we know how aristogenesis works, we have set forth (Osborn: 874, 901) 
a series of principles of its operation, but we have no explanation or hy- 
pothesis to offer as to why it works or what are its causes. In this sense 
aristogenesis can only be challenged by a reéxamination and reinterpreta- 
tion of the paleontological evidence on which it rests. One principle of this 
evidence was discovered in the invertebrates by Waagen in 1869, namely, 
the orderly “‘W. mutations’” of the shells of molluscs which arise in a 
most inconspicuous manner and develop into conspicuous organs. New 
and amplified principles have been discovered in vertebrates as fully 
enumerated in my last two.communications (Osborn: 874, 901) in this 
series of biological inductions from my own paleontological researches 
of the last half century. (See full bibliographic list below.) 

The evidence for aristogenesis or the direct origin from the geneplasm 
of biomechanical adaptations is now overwhelming; Waagen’s principle 
of W. mutations in invertebrates (see Osborn: 353) of 1869 was re-dis- 
covered in 1894 by Osborn (Osborn: 92) in the teeth of vertebrates. It 
has been confirmed and amplified by researches recorded in fifty-one 
contributions between the years 1894 and 1933, chiefly based upon thou- 
sands of independent observations on the evolution of the grinding teeth 
(Osborn: 33) in the primates, in the equines and proboscideans (teeth), 
rhinoceroses (teeth and horns), and titanotheres (teeth and horns). These 
observations are harmonious and concurrent. They have been made purely 
inductively with no other purpose than the direct recording of facts. 
Now they enable us to give a clear definition of aristogenesis as follows: 

Aristogenesis is the gradual, secular, continuous, direct, reactive, adaptive 
origin of new biomechanisms. It is a creative process from the geneplasm 
of entirely new germinal characters. It is the orderly creation of something 
better or more adaptive. Certain lines of descent are distinguished by the 
potentiality of creative origin from the geneplasm of new adaptive biomecha- 
nisms. Germinally predetermined origin of new characters tends toward 
betterment, arises independently, in widely separated geographic areas, at the 
same or different aristogenic rates. Aristogenesis is a secular genetic and 
adaptive reaction rather than an immediate adaptive reaction to new habit 
or environment as in Lamarck’s theory. Aristogenes are independent in origin 
and early development of the Natural Selection principle of Darwin. 

Fourth, since single new biomechanisms originate in the geneplasm 
(a term preferable to germ plasm) they may be called “‘aristogenes’’ as 
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distinguished from ‘‘somatogenes’” which originate in the somaplasm. 
Aristogene is also a better and more distinctive term than ‘‘rectigradation,”’ 
the term first applied (Osborn: 314) to these new biomechanisms. 

Fifth, aristogenesis, as the direct adaptive origin of new biomechanisms, 
is of two kinds: 





Aristagenic and Allometric Origin of Species. Osborn, 1933 
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FIGURE 1 


Aristogenic evolution from Phiomia osborni to Trilophodon macrognathus of the third 
inferior grinding tooth in the long-jawed mastodons migrating from northern Africa 
(Phiomia) to southern India (Trilophodon). 


(a) ARISTOGENES (rectigradations) are governed or predetermined by 
germinal potentiality in certain lines of racial, specific, generic, family 
and ordinal descent. The closer the descent or phyletic kinship the more 
similar are the aristogenes. In extremely close racial lifies like that of the 
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fourteen-million year longirostrine mastodons cited in my last communica- 
tion to the Academy (Osborn: 901), the aristogenes are closely similar 
no matter in what part of the world they arise. Thus from 10 aristogenes 
in the lower Oligocene mastodon (Phiomia osborni of the Faytim of 
northern Egypt) the number rises to 19 in Trilophodon cooperi also in 
Trilophodon pontileviensis of the basal Miocene of France and of India, 
to 20 in the lower Miocene Trilophodon paleindicus, to 25 in the Trilo- 
phodon chinjiensis of the Punjab and 37 in Trilophodon macrognathus, 
the two culminating members of this fourteen-million year aristogenic 
series. 

(b) ALLOMETRONS are not governed or predetermined by germinal 
potentiality in certain lines of racial, specific, generic, family and ordinal 
descent. On the contrary, within species and even within races, for ex- 
ample the modern species of man Homo sapiens, diverse allometrons or 
more or less profound changes of proportion, are independently arising. 

The term allometron (Greek 4)\dovos ‘‘of another sort or kind, different’’® 
and perpov ‘‘to measure in any way’’*) signifies changes of proportion or 
intensity which may be expressed in measurements and indices, for ex- 
ample the proportional changes of the head expressed in the Greek terms 
brachycephaly (broad headed), dolichocephaly (long headed), brachyopy 
(short faced), dolichopy (long faced), bathycephaly (deep skulled), hypsi- 
cephaly (high skulled), etc., or in the limb proportions brachydactyly 
(broad fingered), dolichodactyly (long fingered), leptodactyly (slender 
fingered), dolichopody (long footed), brachiopody (short footed), doli- 
chomely (long limbed), brachymely (short limbed), etc., or in tooth 
proportions brachyodonty (short toothed), hypsodonty (long toothed), 
etc. 

Allometrons are relatively rapid in development, or temporal, whereas 
aristogenes are relatively slow, or secular. The causes of geneplasmic 
allometrons are obscure. In some unknown way they seem to be connected 
both with hormones, and with the habit inheritance principle of Lamarck 
and the natural selection principle of Darwin. 

So far as the modes in which the “origin of species’ is defined by bio- 
mechanisms we are now on absolutely sure ground. This ground is contra- 
Lamarckian and contra-Darwinian. It is also contrary to the neo-Dar- 
winian evolutionary hypotheses of the leading biologists and geneticists 
of our day. 
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ELECTROSMOSIS AND ANOMALOUS OSMOSIS 
By Davin R. Briccs 


Oto S. A. SPRAGUE MEMORIAL INSTITUTE AND THE DEPARTMENT OF PATHOLOGY, 
UNIVERSITY OF CHICAGO 


Communicated June 12, 1933 


That electrosmosis can account for observed anomalies in osmotic trans- 
fer of liquid across a membrane has been held by Girard,! Loeb? and Bar- 
tell. Their measurements have been made upon rates and direction of 
flow of fluid through various membranes bathed on the two sides by elec- 
trolyte solutions of different concentrations. Since, in such cases, the 
membranes were slowly permeable to the electrolytes used, equilibrium 
measurements were not possible. These dynamic measurements, while 
indicating the magnitude of the resultant of the forces acting upon the fluid 
within the membrane at the time of observation, were not of such a nature 
that the magnitude of these forces could be individually evaluated. 

Certain prerequisites to the establishing of an electrosmotic pressure 
across a membrane are necessary. There must exist a difference in elec- 
trical potential across the membrane, there must be an electrokinetic 
potential existing between the wall of the pores of the membrane and the 
fluid filling the pore and a flow of current must take place as a result of 
which the charged, movable layer of liquid in the pore is moved, carrying 
along with it the liquid in the center of the pore. 

Numerous measurements have been made upon the osmotic pressure 
exhibited by solutions containing colloidal electrolytes. In such cases, the 
membrane is impermeable to the colloid ion and the diffusible ions in the 
system arrange themselves at equilibrium in a manner described in Don- 
nan’s* membrane equilibrium theory. Loeb® found that the osmotic pres- 
sure shown by casein in HC! solution was entirely due to the unequal dis- 
tribution of these diffusible ions, i.e., the colloidal micellae exerted no de- 
tectable osmotic pressure. Bjerrum® studied a chromium hydroxide solu- 
tion, and Adair’ studied hemoglobin. Both concluded that, within cer- 
tain ranges of salt, hydrion and colloid concentration the observed osmotic 
pressures were the sum of the pressures derived from the unequal distribu- 
tion of diffusible ions and that derived from the colloid micellae. Ham- 
marsten,® using sodium nucleate, and Samac,’ studying amylopectin and 
ligno sulphuric acid, found that the observed osmotic pressure was less than 
that calculated on the basis of differences in diffusible ion concentration 
and concluded that interionic forces, exerted by the large colloid ion, com- 
pletely or to a large extent destroyed the osmotic acitivity of the small ions 
(the Gegenionen) present (and ionized from the colloid). With such varied 
conclusions it seems probable that some unrecognized force may have been 
acting in most cases to modify the true osmotic pressures. 
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In all such experiments the conditions necessary to the establishment of 
an electrosmotic pressure were present. The membrane potential resulting 
from unequal distribution of diffusible ions on the two sides of the mem- 
brane could function as the electrical potential across the membrane needed 
to cause liquid to be transferred electrosmotically through the pores where, 
undoubtedly, there must exist a {-potential of one sign or the other. The 
circuit through which the current must flow in order that the water may be 
moved electrosmotically through the pore could be completed within the 
pore, i.e., in one direction the current would be carried by the ions making 
up the movable half of the double layer along the pore wall and in the other 
direction by ions within the body of the liquid in.the pore (see Bartell,* 
Freundlich”). At equilibrium, a pressure would be set up as a result of 
the electrosmotic forces acting which would cause the observed osmotic 
pressure to deviate from that pressure which should result from purely 
osmotic forces acting across the membrane. Depending upon the sign of 
the membrane potential and that of the ¢-potential, the observed osmotic 
pressure could be either higher or lower than that calculated on the basis 
of differences in concentration of the osmotically active components of the 
system on the two sides of the membrane. The difference between the 
observed osmotic pressure, P,, and that calculated from the difference in 
osmotic concentrations, P., should then be definable in terms of the elec- 
trosmotic equation— 


_ 2E¢D 


, 
ar? 


P, 


where P,, is the electrosmotic pressure, which, if the theory holds, would be 
equal to P, — P, at equilibrium; £ is the potential across the membrane 
and is equal to the membrane potential derived from unequal distribution 
of ions across the membrane, being calculable, upon the basis of the Don- 
nan theory, from the difference in active concentration of any one measur- 
able ionic species (as H* ions); 1 is the radius of cross-section of the pore 
of the effective membrane; ¢ is the electrokinetic potential at the solid- 
liquid interface within the pore of the effective membrane; and D is the 
dielectric constant of the liquid in the region of the interface. 

A series of measurements of the osmotic pressure of gum arabic indicates 
the validity of this theory. The colloid used was gum arabic which had 
been purified and electrodialyzed free of all ions other than hydrogen ion. 
This electrodialyzed colloid acts as a strong acid (Thomas and Murray)" 
and each gram required 85 X 10~-® equivalents of alkali to neutralize it. 
Osmotic pressure measurements were obtained by placing known amounts 
of the colloid (the Na salt) in small collodion sacs and allowing equilibrium 
to be attained against a constant pressure applied to the solution inside of 
the sac while it was immersed in an external solution. The final concen- 
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tration of colloid within the sac was obtained by weighing the equilibrium 
contents of the sac. The total volume of the water in the system (inside 
and outside the sac) was known in each case and its distribution at equi- 
librium was obtained. Also the total equivalents of each diffusible ion in 
the system were known. After equilibrium had been attained, the H* 
ion concentration was determined electrometrically on the solutions from 
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both sides of the membrane. From this determination the membrane po- 
tential, E, was obtained, as was also the ratio of distribution, R, of diffusible 
ions across the membrane. Then, on the basis of Donnan’s equilibrium 
theory, it was possible to calculate the active concentrations of all diffus- 
ible ions on each side of the membrane. All experiments were made at 
25°C., at which temperature one mole of dissolved material exerts an 
osmotic pressure of 25,280 cm. of HO. Neglecting any osmotic pressure 
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which might arise from the large colloid ion, the difference in the concen- 
trations of diffusible ions multiplied by 25,280 gave the calculated osmotic 
pressure, P., which would result from the unequal distribution of these ions. 
The observed osmotic pressure, P,, was equal to the constant applied pres- 
sure against which equilibrium had been attained. 

The experiments included one series in which the observed pressure was 
varied and the external solution was common to all the sacs, one in which 
the pressure was constant and the pH was varied from 2 to 7 while the 
total equivalents of Cl~ ion in the system was constant, and a third series 
in which the pressure was constant and the concentration of NaCl in the 
system was varied from a very low value up to 0.01 N. In all experiments 
the system was kept on the acid side by addition of HCI in order to elimi- 
nate CO, absorption. e 

In these experiments the only diffusible ions in the system were H+, Nat 
and Cl-. In every case the value of P, was greater than P,. The values 
of E, as derived from pH measurements, where E (in millivolts) = 59.1 
(pH, — pH;), were always such that the external solution was negative to 
the internal solution. Since the electrokinetic potential in the pore must 
be such that the wall was negatively charged to the liquid (whether the 
pore was that of the collodion or that of a gum arabic membrane the wall 
would be negatively charged since both are acidic and have no isoelectric 
point and NaCl or HCl cannot reverse the charge) it appeared that the 
electrosmotic pressure, P,, should act counter to the true osmotic pressure 
Pas the value of P, — P, indicated it did. 

In the electrosmotic equation given above, D could be taken as constant. 
In these and earlier experiments it was found that 7 could also be taken as 
constant. Whether the membrane sac was one of collodion, cellophane or 
parchment, made no difference in the final equilibrium concentration of 
colloid and H+ ion distribution. The pore size of the supporting membrane 
could be varied widely with no effect upon these values. It appears that 
the effective membrane in all cases consisted of a layer of gum arabic laid 
down on the inner wall of the supporting membrane. For this reason the 
value of 7 effective in the electrosmotic equation was that of the gum arabic 
membrane and was therefore constant for all these experiments. That the 
effective membrane consisted of the gum arabic made possible estimates of 
the value of £, also, as seen below. 

By collecting all constant terms of the equation on one side, the equation 


becomes 
ar 
( ee 


Of the variables, E is taken equal to the membrane potential and P, equal 
to P, — P,. An absolute value of ¢ was not obtainable from the data 
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derived from these experiments, but by regarding the effective pores as 
those of the gum arabic membrane and assuming that ¢ would be a func- 
tion of a, the degree of ionization of the primary salt valences of the colloid, 
it was possible to obtain comparative values of ¢. It was assumed, also, 
that the value of a in the bulk of the colloid solution would equal the value 
of a for that part of the colloid making up the effective membrane. The 
equivalent weight and final concentration of the colloid were known. So 
were the final concentrations of H+, Na* and Cl* in the colloid containing 
solution (from calculations described above). It was possible then to cal- 
culate the fraction of the colloid which was ionized, that is, a. 





TABLE 1 
Po P. P, of Ea Ea(ci-]?"" 
Cm. Cm. Cm. {(H*}; [Nath [Crk E — 470-211 
HO H:O H:O X07 X107 X107~ Mv. a P, [C17]; P, 
27.3 54.2 26.9 22950 24260 5000 28.7 0.1855 0.1980 0.2010 0.03980 
37.8 70.7 32.9 26420 27930 4343 32.4 0.1844 0.1815 0.1950 0.03540 
48.7 82.4 33.7 28830 30490 3980 34.7 0.1717 0.1767 0.1915 0.03380 
85.7 117.8 32.1 36050 38110 3184 40.3 0.1522 0.1909 0.1828 0.03490 
97.4 184.8 37.4 39430 41700 2914 42.6 0.1531 0.1744 0.1795 0.03130 
110.5 147.5 37.0 41960 44350 2735 44.2 0.1519 0.1815 0.1770 0.03210 
119.5 157.1 37.6 43830 46360 2617 45.3 0.1520 0.1832 0.1753 0.03210 
127.1 157.5 30.4 53320 47320 4580 39.6 0.1348 0.1756 0.1972 0.03460 
127.1 164.2 37.1 42170 61520 4505 40.1 0.1567 0.1694 0.1968 0.03330 
127.6 171.2 43.6 30280 76400 4398 40.9 0.1805 0.1694 0.1958 0.03315 
127.5 186.2 58.7 19630 93400 4191 42.0 0.2262 0.1620 0.1940 0.03140 
128.1 201.0 72.9 9930 109470 4039 43.3 0.3025 0.1796 0.1922 0.03450 
128.8 250.0 121.2 4062 136000 3558 46.5 0.4211 0.1618 0.1870 0.03025 
129.5 388.3 258.8 933.0 194500 2503 55.9 0.7170 0.1550 0.1736 0.02690 
129.0 472.0 343.0 16.1 229500 2200 59.2 0.9470 0.1637 0.1690 0.02765 
131.5 320.5 189.0 49.5 129800 17.4 114.4 0.822 0.4975 0.0624 0.03105 
130.8 344.1 213.3 62.2 139800 25.8 111.0 0.875 0.4561 0.0665 0.03031 
131.5 351.5 220.0 74.1 146300 92 94.8 0.859 0.3701 0.0875 0.03240 
131.2 381.9 250.7 171 165950 309 81.5 0.876 0.2850 0.1128 0.03215 
131.0 398.5 267.5 235 177000 568 74.4 0.873 0.2430 0.1282 0.03114 
130.4 278.4 148.0 1905 163200 5540 43.2 0.522 0.1523 0.2080 0.03168 
130.4 217.0 86.6 7672 190400 23130 27.2 0.349 0.1097 0.2792 0.03062 
130.1 195.2 65.1 12020 215300 39530 22.2 0.295 0.1005 0.3126 0.03141 


If the theory holds that the pressure, P, — P, = P,, could be defined by the 
E 

electrosmotic forces acting, then the value = would be a constant through- 

out all the experiments (disregarding any osmotic pressure derived from 


E 
the colloid micellae). If { was a function only of a then also the value = 


should be constant. 
In the table of results, it is seen that, where the [Cl—]; value was held 
Ea 


P was equal to a constant. When, however, 
x 


nearly constant, the value 
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the foreign salt concentration in the colloid solution was changed widely 


Ea 
the value P. varied in a manner such that at higher concentrations of salt 
x 


(designated as [Cl~]; = = became lower. That this effect was not due to 


a change in 7, the radius ee cross-section of the pore, as a result of dehydra- 
tion of the micellae making up the effective membrane, was ruled out by 
adding a nonelectrolyte, ethyl alcohol, in place of NaCl. The concentra- 
tion of alcohol was varied over a range 10 times as great as the range of 


concentration in which the effect of NaCl was observed and no variation in 


Ea Ea 
P. occurred. It became obvious that the observed effect on P. with in- 
x x 


creased salt concentration was an effect upon ¢. That is, ¢ was a function 
of a but also a function of the foreign salt concentration in the solution. 
This latter function of was presumably derived from oriented adsorption 
of the salt by the colloid and should therefore be definable in terms of the 
adsorption isotherm and ¢ would be a function of the amount, a, adsorbed. 
The value of (a) could not be derived from the data at hand but the manner 
in — (a) varied with the salt concentration could be obtained by plot- 


deg P, “ against [Cl—]; and a comparative value of (a) could be obtained in 


terms of [Cl~];. In the graph shown, wee P, is plotted against log [Cl—] 


and from the equation for the straight line thus obtained it appears that 
E a [Cl- If +211 

P, 
plotted in the figure and shown in the last column of the table.) 

It appears, then, that ¢ is a function of a, the degree of ionization of the 
primary salt valencies of the colloid, and of (a) the degree of oriented ad- 
sorption of the foreign salt by the colloid. In this latter function of ¢ 
appears the possibility of the lyotropic series playing a part in the electros- 
motic forces acting, because oriented adsorption is a function of the rela- 
tive adsorbability of both cation and anion of the salt and will not be the 
same for any two ions of equal sign and valence. 

Without direct measurement of ¢ but upon the assumption that ¢ is a 
function only of a and (a), as the calculations indicate, it becomes apparent 
that anomalous osmosis observed for systems containing the colloid, gum 
arabic, can be fully defined in terms of the electrosmotic forces existing 
therein. The osmotic pressure derived from the colloid itself is so small 
as to be undetectable. 

It is obvious that electrosmosis is an important osmotic factor which 
must be evaluated in every case where the osmotic properties of colloid 
solutions are being studied, and that its evaluation will furnish a satis- 


is a constant. (This value for the earlier experiments is also 
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factory explanation to most cases of anomalous osmosis in such systems. 

Experiments with the Ca salt of arabic acid and with increasing concen- 
tration of CaCl, indicate that at a low concentration of the divalent cation 
salt the sign of ¢ is reversed and negative osmosis occurs, i.e., the value of 
P, is greater than P,. This and lyotropic effects are being studied at the 
present time. 
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NOTE ON THE STURM-LIOUVILLE EIGENVALUE-EIGEN- 
FUNCTION PROBLEM WITH SINGULAR END-POINTS 


By Epwin C. KEMBLE 


RESEARCH LABORATORY OF PHYSICS, HARVARD UNIVERSITY 


Communicated June 7, 1933 


The general systematic theory of Sturm-Liouville (S.L.) eigenvalue- 
eigenfunction problems as discussed in standard texts such as Courant- 
Hilbert and Ince is restricted to cases involving regular end-points, finite 
fundamental intervals and homogeneous boundary conditions. The most 
important S.L. problems of mathematical physics lie outside this restricted 
class, however, and the question of the existence of discrete eigenvalues 
as well as that of the completeness of the system of corresponding eigen- 
functions is usually dealt with for each problem separately. The general 
case of an S.L. problem with one regular boundary point in the finite re- 
gion and a singular boundary point at infinity has been treated by Weyl! 
in a well-known paper whose primary objective was the development of an 
expansion theorem for the continuous spectrum. The class of problems 
directly treated in Weyl’s paper is not perfectly general and his technique 
is too elaborate to be readily accessible to physicists. Hence the writer has 
considered it worth while to develop a general treatment of the discrete 
spectra of S.L. problems with singular end-points which involves only ele- 
mentary methods. The procedure is only sketched in this note which the 
author hopes to follow with a more detailed publication elsewhere. 
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Consider the S.L. equation 


d d 
a p(x) of — g(x)y + Ap(x)y = 0 (1) 


in which ) is a parameter while p(x), q(x), p(x) are real continuous functions 
of x in the fundamental intervala <x < b. p(x) and p(x) are required to 
be positive at every interior point of this interval. The boundary points 
a, b are singular and either or both of these points may lie at infinity. The 
solution of this equation becomes an eigenvalue-eigenfunction problem if 
we impose suitable boundary conditions. Those defined by physical ap- 
plications vary from case to case and it is therefore convenient to introduce 
a set of standard mathematical boundary conditions to which the various 
physical conditions can usually be reduced. 
To this end we require that there shall exist a function g(x) such that 


Lim { p(x) | v(x) |? g-1 (x)} = 0; (2) 
t=—a 
p(x) | -y’(x) |2 g(x) < M; 0<x-a<B (3) 


where M and B are positive quantities independent of x. g(x) is to be in- 
dependent of the parameter \ whose eigenvalues we seek, and can usually 
be set equal to unity. This condition we call the singular point boundary 
condition (s.p.b.c.) for the end-point a. A corresponding condition is to 
be applied to 0. 

The s.p.b.c. has the following important properties. 


(a) If (x) and u(x) conform to the s.p.b.c. atx = a 
Lim {p(x)muf’} = Lim {p(x)ujus} = 0. (4) 
*—>a ‘i—>a 


(The asterisk denotes the complex conjugate.) 

(b) If u(x) and u(x) conform to the s.p.b.c. at x = a, any linear com- 
bination, say, w(x) = au; + Bus, will conform to the same condition. This 
is a consequence of the inequality 


| w |? <2 fl coms |? + | Bue |*}. 


(c) Let y:(x) and ye(x) denote two linearly independent integral curves 
of Eq. (1) constituting a fundamental system. Then both of these func- 
tions cannot conform to the s.p.b.c. at either end-point, for it follows from 
the linear independence that the Wronskian 


Wyiny2) = P(r, — nye) ¥ O, 


whereas if both y, and ye satisfy the s.p.b.c. at either end-point, it follows 
from (a) that W must vanish at the corresponding boundary point. 
An important corollary on (a) is that any two eigenfunctions y,(x), 
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ye(x) with different eigenvalues \, and 2 must satisfy the orthogonality 


relation 
b 
a pyrygdx = 0. 


From (c) we draw the further inference that there cannot be two linearly 
independent eigenfunctions with the same eigenvalue. 

We approach the question of the existence of eigenvalues by postulating 
the existence of an interval G of \-values, say, \’ £ A < X” such that for 
every \ in G there exists a pair of integral curves u(x), v, (x) which con- 
form to the s.p.b.c. at a and b, respectively, and which do not have an in- 
finite number of nodes in the neighborhood of these points. It follows from 
this assumption that there exist continuous functions of x and X, say, 
u(x,d), v(x,A) which for any given ) in G satisfy the conditions laid down for 
UM and 2. 

If the Wronskian W[u(x,d), v(x,A)] vanishes for any value of \ in G, 
say, Xi, “(x,A1) and v(x,di) are linearly dependent and conform to the 
s.p.b.c. at both boundary points. Hence , is an eigenvalue and u(x,\,), 
v(x,d1) are eigenfunctions. It is easy to prove by standard methods that as 
d increases the nodes of u(x,\) move continuously to the left while those of 
v(x,) move to the right. Whenever a node of u comes into coincidence 
with a node of v, they all coincide in pairs, the Wronskian vanishes and 
we have an eigenvalue. 

Let n” be the number of interior nodes of u(x,X”) [upper limit of G] 
and let ’ be the number of interior nodes of u(x,\’). For simplicity we 
assume that \” and 2’ are not eigenvalues., It follows readily from the 
above argument that there exist n” — mn’ eigenvalues in the interval G, 
each having an eigenfunction with one more interior node than the eigen- 
function of next lower eigenvalue. The case where n’ = 0 requires special 
discussion but forms no exception to the rule. If” > 0, and if the lower 
‘limit of G is at} = — ©, there is a minimum eigenvalue with a nodeless 
eigenfunction. If the upper limit of G is at \ = +, the sequence of 


ascending eigenvalues is denumerably infiniteand Lim A, = ©. 
ti > © 
The S.L. eigenvalue-eigenfunction problem is usually reducible to varia- 
tional form. One of these forms consists in finding stationary values of 


quadratic integral 


'b 
Tyra) =f {ply’ |? + a] y|? — dol y|?} ax, (5) 





the admissible comparison functions being subject to boundary conditions 
which insure that 


Lim {p[dyy*’ + dy*y’]} = Lim {pldyy*’ + dy*y’]} (6) 
ts—> <a x— > bd 
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and to the requirement that they yield convergent values of J. The 
Euler’s equation of (5) is (1), the s.p.b.c. take care of (6) and the eigen- 
functions of (1) yield the value zero for J. If they also yield convergent 


b 
values of the normalization integral N = f p | y |? dx, they must yield con- 


b 
vergent values of Q = f (p | y’ |2 +q | y |2)dev. In this case the eigenvalues 


of \ are stationary values of Q/N for comparison functions subject to the 
same restrictions as before plus the requirement that N shall converge. 
In applications to wave mechanics this latter requirement is usually one of 
the physical conditions to which the s.p.b.c. are assumed equivalent. 
Changing the dependent variable from y to w by the transformation 
w = p’*y converts the variational problem into one involving fixed end- 
points, since w(a) = w(b) = 0. It isa theorem of the calculus of variations 


#1 
that ify = U(x) isa nodeless extremal of an integral D[y] = f F(x,y,y')dx 
x0 


of the form (5) subject to fixed end-point conditions, and if F(x,y,8) > O 
for every point (x,y) in the neighborhood of U(x) and for every finite 8, 
then U actually minimizes D.* In proving this theorem the end-points 
Xo and x, are assumed regular. Were it not for this circumstance we could 
identify J[y,\,a,b] with D and show that when Eq. (1) has a nodeless ex- 
tremal yo yielding a convergent N[yo] the corresponding eigenvalue Xo is a 
minimum for Q/N. 

To overcome the difficulty regarding the singularities at the end-points 
it is necessary to introduce an infinite sequence of modified variational 
problems A;, As, ... A,, ... with regular end-points a,, >, inside the 
fundamental interval a < x < 5b and with fixed end-point conditions. 
This sequence can be so chosen that the lowest eigenvalue and eigenfunc- 
tion converge on \y and w = p’/*yy as m becomes infinite. By application 
of the above theorem to the members of this sequence we can prove that 
yo actually minimizes Q/N. Moreover, the minimum is absolute since for 
quadratic integrals of the type under consideration the increment due to 
any variation in the argument function is equal to the sum of the first and 
second variations of the integral and hence is positive or zero. 

The higher eigenvalues of Eq. (1) can be treated in the same way,’ each 
being proved to yield the minimum value of Q/N for all comparison func- 
tions which obey the s.p.b.c., yield convergent integrals and are orthogonal 
to all eigenfunctions of lower eigenvalue. When this point is established, 
one can adopt the procedure of Courant-Hilbert* to show that when the 
interval G extends from Xo (or a lower value of \) to ~, the eigenfunctions 
form a complete set. 

The question of the range of the discrete spectrum of eigenvalues of 
Eq. (1) and the completeness of its system of eigenfunctions is thus reduced 
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to the question of the validity of our postulates regarding the equivalence 
of the physical boundary conditions to the s.p.b.c. and the existence and 
range of the interval G. These points must be investigated separately 
for each individual case, using primarily the standard methods for the 
study of the behavior of solutions of a differential equation in the neighbor- 
hood of a singular point. 

The writer is greatly indebted to his colleague, Prof. Marston Morse, for 
numerous suggestions and much helpful discussion in connection with the 
proof that the eigenvalues of \ are minima of Q/N for suitably restricted 
classes of admissible comparison functions. 

1 Weyl, Hermann, Math. Ann., 68, 220 (1910). 

2 Bolza, O., Lectures on the Calculus of Variations, Chicago, p. 96 (1904). 

3 Kuen Sen Hu, Thesis, Theorem 10.3 [Contributions to the Calculus of Variations, 


Univ. of Chicago Press, 1933]. 
4 Courant, R., and Hilbert, D., Methoden der Math. Physik, Berlin, 1931, Kap. 6, § 3. 


STUDIES ON THE GROWTH HORMONE OF PLANTS. III. 
THE INHIBITING ACTION OF THE GROWTH SUBSTANCE 
ON BUD DEVELOPMENT 


By KENNETH V. THIMANN AND FOLKE SKOOG 


WrtitiamM G. KERCKHOFF LABORATORIES OF THE BIOLOGICAL SCIENCES, CALIFORNIA 
INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 


Communicated April 24, 1933 


It has long been known that when the growing point of a young di- 
cotyledonous plant is removed, the axillary buds on the stem below it 
begin to develop. As long as the terminal bud is present, the development 
of the axillary buds is inhibited. A lateral bud may also be inhibited 
by the rapid growth of another lateral above it, as we have found in our 
experiments, or opposite it, as in those of Dostél (1926). Furthermore, 
Snow (1929a) has shown that the inhibition, in Pisum, is principally due 
to the young leaves in the developing bud. The evidence indicates that 
this inhibition is probably caused by a special substance (see Snow, 1929b). 
There was reason to believe that this inhibiting substance is of the same 
nature as the growth-promoting substance of Avena coleoptiles, and 
the experiments to be described here confirm this belief. 

1. Methods——Young plants of Vicia faba, 4 to 6 weeks old, grown in 
the light, were used. The lowest lateral buds of these plants grow with 
great regularity. The plants selected were of approximately equal height 
and with as nearly as possible equal numbers of buds and leaves. The 
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growth substance used was obtained from the growth of Rhizopus suinus 
(Thimann and Dolk, 1933) and had an activity of about 2.10-* mg. per 
plant unit. The activity of this growth substance and of that obtained 
from the Vicia plants was tested on Avena coleoptiles under the standard 
conditions described by Went (1928), using the definition of units as given 
by Dolk and Thimann (1932). 

2. Production of Growth Substance by Vicia faba.—It was first found 
that the terminal bud produces growth substance in rather large quantities. 
From terminal buds of young plants 12 cm. high, 30 to 40 plant units 
diffused out into agar blocks in an hour. From the buds of older plants 
the amount obtainable was less. The undeveloped lateral buds produce 
practically no growth substance. When, however, the plant is decapitated 
and the lateral buds are undergoing rapid development, each bud produces 
almost one-half as much growth substance as the terminal bud of an 
intact plant. Small amounts of growth substance are also produced by 
the leaves, less than one plant unit per hour being obtainable from the 
oldest leaves, and about 5 plant units per hour from the youngest leaves 
under the conditions of the experiment. It is therefore clear that growth 
substance production is associated with the actively developing parts of 
the plants, and that rather large quantities of the substance are regularly 
passing into the stem. 

3. Inhibition of Bud Development by Growth Substance——In order to 
make quantitative determinations of the effect of applying growth sub- 
stance to the plant, it was necessary first to determine the time required 
for growth substance to enter the plant from an agar block. By plotting 
the amount of growth substance remaining in agar blocks after different 
periods of application to the stem, it was found that 6 hours were necessary 
for a complete transference of growth substance into the stem. 

Agar blocks containing growth substance were therefore applied to the 
tops of decapitated plants every 6 hours. Intact plants, and also plants 
to which blocks of plain agar were being applied, were used as controls. 
Measurements of the length of the lateral buds, and also of the main stem, 
were made daily. When the amount of growth substance applied was of 
the same order as that diffusing from the terminal bud into agar (160 
plant units every 6 hours), a slight but definite inhibition of the develop- 
ment of the lateral buds was observed. When, however, the amount of 
growth substance applied was larger than that diffusing from the terminal 
bud (1400-1700 plant units every 6 hours), the development of lateral 
buds was completely inhibited. Thus, in one experiment, the mean 
increase in length of the lateral buds in 8 days was 3.4 + 0.5 mm. in in- 
tact plants, 16.0 + 2.7 mm. in decapitated controls and 1.8 + 0.6 mm. 
in plants to which 1670 units of growth substance had been applied. 
Several such experiments were carried out, each comprising thirty to forty 
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buds, and the amount of inhibition was similar in each case. The inhibition 
also takes place in the dark. 

The necessity for applying larger amounts of growth substance than 
can be obtained from the terminal bud is fully justified on the ground that 
the application is generalized over the whole stem surface, while the normal 
supply from the tip is localized in the conducting tissue and therefore 
more effective in its action. Furthermore, there is evidence that wound 
substances may inactivate a part of the applied growth substance (Kisser, 
et al., 1931). The inhibition is not to be ascribed to any damage, since 
the plants remained in good condition throughout the experiment, as 
was shown by the rapid development of the lateral buds as soon as the 
application of growth substance was stopped. While it is possible that the 
inhibition is due, not to the growth substance, but to another compound 
of similar nature present in the active concentrates, the present experiments 
make it probable that it is the growth substance itself, produced in the 
terminal bud, which inhibits lateral bud development. 

Although it seems paradoxical that a substance promoting celi extension 
can also act as an inhibitor, this fact provides an explanation for much of 
the earlier work on inhibition, such as the experiments of Dostal, Snow 
(1929a) and Weiskopf (1927). The probable mechanism of the effect, 
together with a more detailed account of these experiments, will be pub- 
lished elsewhere. 


Dolk, H. E., and Thimann, K. V., Proc. Nat. Acad. Sct., 18, 30 (1932). 

Dostél, R., Acta. Sci. Nat. Morav., 3, 83 (1926). 

Kisser, J., Stasser, Kiffe and Géllner, Anz. Akad. Wiss. Wien, 275 (1981). 

Snow, R., (a) New Phytologist, 28, 345 (1929); (b) Annals of Botany, 43, 261 (1929). 
Thimann, K. V., and Dolk, H. E., Biol. Zentr., 53, 49 (1933). 

Weiskopf, B., Publ. biol. l’école hautes études vétérinaires Brno, 6, No. 4 (1927). 
Went, F. W., Rec. trav. bot. Néerl., 25, 1 (1928). 
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STUDIES ON THE GROWTH HORMONE OF PLANTS. IV. 
ON THE MECHANISM OF THE ACTION 


By JAMES BONNER 


Wn. G. KERCKHOFF LABORATORIES OF THE BIOLOGICAL SCIENCES, CALIFORNIA IN- 
STITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 


Communicated May 19, 1933 


Introduction.—-It has been known for some time that among the higher 
plants cell elongation as well as tropistic response is governed by a well- 
defined chemical substance which is secreted by young leaves and by the 
apical portions of stems.’? The mechanism by which this substance 
exerts its effect upon plant cells has, however, remained obscure and it 
is relatively recently that an attack upon this problem has been made. 

In an earlier paper it was shown that the amounts of this ‘‘growth 
substance” which enter the Avena coleoptile are far too small to bear 
any simple stoichiometrical relation to the substances which are formed 
in the cell wall during the resultant growth.* The growth substance is 
therefore a true hormone in that it is active in minute amounts and in 
some indirect manner. The work of Heyn,‘ of Heyn and van Overbeek,°® 
and of Séding* has shown that Avena coleoptiles which are artificially 
supplied with growth substance exhibit a marked increase in the plasticity 
of their cell walls. From this result it seemed possible that the growth 
substance itself might possess the property of being able, merely by its 
presence, to decrease the viscosity of the fluid in which the cellulose micelles 
of the cell wall are embedded. The cell wall could then be readily stretched 
by the force at hand, which would, under normal conditions, be the force 
due to the osmotic pressure of the cell contents. It was found in this 
laboratory, however, that if the osmotic force be replaced by the force 
due to a suitable weight, no increase in plasticity of the coleoptile results 
from the application of growth substance. It seemed, therefore, more 
probable that the hormone exerts its influence upon the cell wall through 
the mediation of the cell protoplasm. The present paper will describe 
briefly experiments which demonstrate that this is actually the case. 
A fuller account will appear elsewhere. 

Materials and Methods.—Avena plants of the pure line, ‘‘Sieges Hafer,”’ 
kindly supplied by Dr. Akermann of Svaléf, were used throughout. The 
plants were grown under the usual rigidly controlled conditions of tempera- 
ture and humidity.’ All experiments were carried out under red light, to 
avoid photropic stimulation of the plants. 

The conventional method of applying growth substance dissolved 
in agar blocks to the stumps of decapitated Avena coleoptiles is not well 
adapted to the study of matters concerned with metabolism. It was 
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found, however, that short sections (three mm. long) cut from the coleoptile 
and suspended in growth substance of suitable concentration elongated 
rapidly, even in the absence of their usual food source (i.e., the seed). 
The various substances whose effects upon the action of growth substance 
are to be investigated may be conveniently added to the growth substance 
solution in which the sections are suspended. In order to obtain sections 
which were initially relatively free of growth substance, the hormone- 
producing tip was removed from each plant two hours before the short 
sections were cut from the coleoptile. 

Experimental.—Avena coleoptiles have not previously been grown, as in 
the present case, suspended in solutions of growth substance. It is there- 
fore of interest to enumerate, as follows, some of the factors which were 
found to control elongation under these conditions: 

(a) Sections immersed in water or solution of growth substance do not 
exhibit ‘‘regeneration”’ of the capacity to produce growth substance which 
is characteristic of decapitated plants grown in the usual manner with only 
their roots in water. 

(b) Sections from near the apex of a previously decapitated coleoptile 
elongate more rapidly in growth substance of a given concentration than 
sections from near the base of the same plant. This indicates that cells 
which have already. undergone considerable elongation are less reactive 
to the hormone than those which have elongated only slightly. 

(c) The rate and extent of the elongation vary greatly with the growth 
substance concentration. There is an optimum concentration above 
which the hormone is toxic and may even bring about a loss of turgor. 

The growth of coleoptile sections under the influence of growth sub- 
stance was found to be completely inhibited by 10-* normal KCN, and 
by 0.05 per cent phenyl-urethane. This inhibition suggested a relation 
between the action of growth substance and respiratory activity. Coleop- 
tile sections were, therefore, immersed in a solution of growth substance 
under an atmosphere of nitrogen. The sections were not harmed by a 
prolonged lack of oxygen since they grew normally if placed in air. Never- 
theless, neither elongation nor any action of growth substance preliminary 
to elongation took place in nitrogen. This result might be due either to a 
failure of the hormone to penetrate into the cells in the absence of oxidative 
processes, or to the participation of growth substance itself in an oxidative 
process. In the first case it would hardly be expected that the respiration 
of a section should be effected by the addition of the hormone, since Steward 
has shown that, in the case of potato tissue at least, the uptake of ions, 
while dependent upon aerobic metabolism, nevertheless does not increase 
the production of carbon dioxide.” In the second case it would be difficult 
to predict a priori whether or not an increase in respiratory activity would 
occur upon the addition of growth substance. 
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A preliminary study of the effect of growth substance upon the respira- 
tion of coleoptile sections was made, using the standard Warburg tech- 
nique.* The results may be summarized as follows: 

(a) Growth substance in low concentrations stimulates the respiration 
of coleoptile sections by as much as 27 per cent during the first two hours 
after its addition. In high concentrations growth substance greatly in- 
hibits the respiration. There is apparently some discrepancy between the 
effects of hormone concentration on growth and on respiration, although 
the optimal concentrations are similar. 

(b) The respiration of sections is inhibited by the same concentrations 
of KCN and of phenyl-urethane which were found to inhibit elongation 
under the influence of growth substance. 

(c) The additional respiration observed upon the addition of growth 
substance exhibits the same sensitivity to KCN and to phenyl-urethane 
as does the normal respiration and is therefore probably identical with it. 

(d) Growth substance which has been inactivated for growth by 
oxidation with hydrogen peroxide does not stimulate respiration. 

Discussion.—From the foregoing results it seems possible that an in- 
crease in respiration is one member in the chain of processes by which 
growth substance brings about elongation. It is as yet unsafe to conclude 
that this is actually the case. In the first place, the hormone preparation 
used was, while of relatively high purity, still only approximately one 
per cent hormone. It is possible that an impurity in the preparation is 
responsible for the stimulation of respiration. In the second place, this 
stimulation may be a secondary phenomenon attending the presence of 
growth substance in the cell, and have in reality no bearing upon growth. 
These questions, as well as a possible relation of respiratory rate to cell 
wall plasticity, are now being investigated. 

The author wishes to express his appreciation to Dr. Robert Emerson 
for his many suggestions and for the use of his laboratory, in which the 
measurements of respiration were made, and to thank Dr. F. Went and Dr. 
Kenneth V. Thimann for their help and suggestions. 
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THE ENERGY OF UREA SYNTHESIS. II. THE EFFECT OF 
VARYING HYDROGEN ION CONCENTRATION WITH 
DIFFERENT METABOLITES 


By HENRY BorRSOOK AND GEOFFREY KEIGHLEY 


Witi1am G. KERCKHOFF LABORATORIES OF THE BIOLOGICAL SCIENCES, CALIFORNIA 
INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 


Communicated June 13, 1933 


In a previous communication! we have reported that an increased oxy- 
gen consumption can be observed accompanying the synthesis of urea by 
liver slices from ammonium bicarbonate. This observation is now 
confirmed, and the results of additional experiments are presented which 
elucidate to some extent some of the general features of the mechanism of 
the transfer of energy in this and probably other coupled reactions. 

The technique employed here was similar to that in the experiments re- 
ported previously, except that the vessels used were two and one-half times 
smaller, and the liver slices were placed in 3 cc. instead of 15 cc. of fluid. 
The composition of the Ringers solution was varied according to the pH 
desired. The hydrogen ion concentrations were measured colorimetrically 
before and after each experiment. The ornithine and the ammonium 
chloride, wherever it was used, were first brought to pH 7.4 with sodium 
hydroxide. Qo, is the number of cubic millimeters of O2 at standard 
temperature and pressure used per milligram of dry weight of tissue 
(dried at 100°C.) per hour. Qyre, is the number of cubic millimeters of 
CO, at standard temperature and pressure, liberated at the end of the 
experiment at pH 5.0 by urease per milligram of dry weight of tissue per 
hour. 

Some of the results obtained are given in the table. They show that with 
no metabolite added or with glucose, at pH 7.4 and 8.4 one molecule of 
additional oxygen is used for every molecule of urea synthesized. With 
dl lactate and with succinate at pH 8.4 the results were the same. On 
the other hand, at pH 7.4 with the latter two metabolites there was much 
less additional oxygen used in the synthesis of urea. These observations 
resolve the apparent discrepancy between our previous observations and 
those of Krebs and Henseleit.? 

The latter authors worked for the most part at pH 7.4. The one experi- 
ment whose protocoll on oxygen consumption is given was carried out at 
this pH with dl lactate as metabolite. The Ringers solution was the same 
as in experiments 1 and 2 above, and a gas mixture containing 5 per cent 
CO, and 95 per cent O. was used. The CO: was not absorbed and so main- 
tained the hydrogen ion concentration of the Ringers solution at pH 7.4. 
In our experiments the gas was pure O2, and the CO: was continually ab- 
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sorbed. As a result, in our previous experiments, a large fraction of the 
bicarbonate was decomposed, producing an alkaline reaction until finally 
a balance was attained between the CO, lost by the solution and that 
produced by the respiration of the tissue. We were, of course, aware of 
this induced alkalinity, but these conditions were preferred in those ex- 
ploratory experiments because they permitted us to observe the rate of 
respiration throughout the experiment; and because larger amounts of 
urea are synthesized under these conditions than when the pH was main- 
tained at 7.4 by a phosphate buffer. Warburg has shown that between 
pH 6.5 and 8.5 the respiration of liver slices is not affected by change in 
hydrogen ion concentration,* which our own experience confirmed. 

In order to obviate a lengthy discussion of the different hypotheses 
which might be invoked to explain these results, we shall describe briefly 
the results of some experiments in a different field which bear directly upon 
the question of the mechanism of the transfer of energy in biological coupled 
reactions.* In these experiments it was found that in order to obtain 
under anaerobic conditions the synthesis of lactate from pyruvate through 
the oxidation of formate by toluene treated bacteria, it was necessary to 
add a small amount of some substance (usually a dye) which is capable of 
oxidation or reduction by both the formic acid and lactic acid dehydro- 
genase mechanisms. In the absence of the dye the formic acid was not 
oxidized and the pyruvic acid was not reduced. Similar results were 
obtained with the formation of succinic acid from fumaric acid at the ex- 
pense of the oxidation of lactic acid to pyruvic acid. The transfer of 
energy here was effected by the dye, which passed in the reduced state 
from the locus where, for example, formic acid was oxidized and the dye 
reduced, to the center where, by oxidation of the dye, pyruvic acid was 
reduced to lactic acid. 

The simplest hypothesis accounting for all of the facts at present seems 
to us as follows: The coupled reaction in which urea is synthesized from 
ammonia is essentially the same with all the metabolites used, at pH 7.4, 
as well as at pH 8.4; because in all cases the synthesis of urea is acceler- 
ated by ornithine; and because with lactate and succinate at pH 8.4 the 
ratio of increased oxygen used to urea synthesized is the same as with glu- 
cose or no added metabolite at pH 7.4 and pH 8.4. Here the energy 
produced is nearly 10 times the amount necessary, and is derived from 
specific oxidations and not from the general respiration of the cell. The 
synthesis, i.e., the transfer of energy occurs only at certain specific centers. 
This is shown by the absence of any correlation between the total rate of 
metabolism, in the absence of urea synthesis, and the amount of urea syn- 
thesized when ammonia and ornithine are added. This is further sup- 

* These experiments were carried out with Mr. H. F. Schott in this laboratory, and 
will be submitted for publication shortly. 
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ported by the experiments on the reductive formation of succinic acid and 
lactic acid mentioned above; and by the demonstration by Krebs and 
Henseleit that urea synthesis is obtained only with intact cells; when the 
structure of the liver cells is destroyed by grinding, though respiration still 
occurs there is no longer any synthesis of urea. This evidence for the 
existence of specific centers also strengthens the first point that the essen- 
tial mechanism of the coupled reaction is the same under all conditions. 

It follows therefore that with lactate and succinate at pH 7.4 (but not at 
pH 8.4) other reactions requiring energy are in progress. In these the 
energy supplied to them is greatly in excess of that required. When 
ornithine and ammonia are added some of this surplus energy, which is 
possibly used in other reactions or which passes off as heat, is used in the 
synthesis of urea. The result is that under these conditions the increased 
oxygen consumption accompanying the synthesis of urea is less than under 
other conditions, e.g., at pH 8.4 or with glucose instead of lactate or suc- 
cinate, when the increased oxygen consumption is proportional to the 
amount of urea synthesized. 

The synthesis of glycogen suggests itself as an obvious possibility for this 
first-mentioned reaction requiring energy upon which the synthesis of urea 
is superimposed. It has been shown that both lactate and succinate may 
be converted to glycogen in the liver.“* In the case of the synthesis of 
glycogen from lactate it is well known that more energy is produced than 
the theoretical minimum necessary. It is probable that a similar state of 
affairs exists in the synthesis of glycogen from succinate. 

It is necessary to the hypothesis advanced above that the metabolite 
added to the Ringers solution, particularly the lactate and succinate, be 
used. The results show that, in the absence of ammonia, at all the hydro- 
gen ion concentrations tried, the rate of respiration with the latter two me- 
tabolites was much higher than when no metabolite, or when glucose was 
added. 

It is clear from the above data that the total energy change is not 
governed by the energy requirements of the reaction in which AF is 
positive. The total quantity of energy produced is much greater than the 
necessary minimum amount computed from solely thermodynamical 
considerations. It is probable that this is a very common, if not a uni- 
versal feature of biological coupled reactions. The production of a large 
amount of excess energy in the synthesis of glycogen from lactic acid may 
be cited as another example. The 1:1 ratio between additional oxygen 
used and urea synthesized, suggests that the factor governing the total 
energy change here is a stoichiometrical one, the consumption of a molecule 
of oxygen really representing the production of a molecule of CO, at a locus 
in the cell, and in such a state that it can be combined with ammonia and 
the other components necessary for the formation of urea. 
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Summary.—Experimental results are presented which indicate (1) that 
one molecule of oxygen is used for every molecule of urea synthesized; 
(2) that the total energy change is governed by the stoichiometrical 
rather than by the energy requirements of the reaction in which the 
free energy change (AF) is positive; (3) the existence of specific centers 
and specific mechanisms in the cell whereby energy is transferred from 
_ specific energy liberating reactions (i.e., not the general respiration of the 
cell) to the energy absorbing reaction; and (4) that this specificity does not 
exclude mechanisms whereby the energy produced by an energy liberating 
reaction may be shared by more than one energy absorbing reaction. 


1 Borsook, H., and Keighley, G. L., Proc. Nat. Acad. Sci., 19, 626-631 (1933). 

? Krebs, H. A., and Henseleit, K., Zeit. phystol. Chem., 210, 33 (1932). 

3 Warburg, O., Biochem. Z., 143, 317 (1923). 

4 Himwich, H. E., Koskoff, Y. D., and Nahum, L. H., J. Biol. Chem., 85, 571 (1930); 
Bayo, C. Pi-Sufier, and Pi, J. F., Biochem. Z., 242, 306 (1931). 

5 Ponsford, A. P., and Smedley-Maclean, I., Biochem. J., 26, 1340 (1932). 


ON THE TEMPERATURE-REGULATORY FUNCTION OF 
“SPONTANEOUS” ACTIVITY IN THE MOUSE* 


By T. J. B. Stier 
LABORATORY OF GENERAL PHYSIOLOGY, HARVARD UNIVERSITY 


Communicated June 9, 1933 


When every effort is made to eliminate excitation in the external en- 
vironment of a young mouse, movements still occur—the so-called ‘‘spon- 
taneous” movements. These movements do not occur haphazardly, but 
with a definite form, namely, in regularly recurring periods of movement 
alternating with periods of quiet. The periods of activity and of quiescence 
are of statistically constant length, having a mean probable error, as per- 
centage of the mean, of 5 per cent for activity, of 10 per cent for quies- 
cence, at 21 + 0.01°C. (Stier, 1930). The mice are still poikilothermic at 
the age of two days; the frequency of occurrence of activity cycles in- 
creases with temperature in the manner of a chemical reaction velocity 
(u, in the Arrhenius equation In K«y/RT, was found to be about 25,300 
calories; cf. Stier, 1930). The duration of a period of activity, however, 
decreases as the body temperature increases, e.g., at 16.5°C. the average 
duration of a single activity period was 0.7 min., at 29.5°C. the average 
duration was 0.3 min. 

The average duration of activity periods times the number of periods 
per hour yields a quantity which should be roughly proportional to the rate 
at which heat is produced by muscular movements of the young mouse. 
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Calculating this product for various body temperatures in the range 15° 
to 32.5°C. gives an indication of the heat contributed by the active muscles 
of the young mouse over the range of body temperatures tested (cf. Fig. 1). 
The internal temperatures were measured by a thermopile (cf. Stier, 1930; 
and Stier and Pincus, 1927-28). At the moment we do not know whether 
the rate of heat production at any one temperature remained constant 
during each period of activity. The kymograms of the spontaneous 
movements show a uniform vigor of movements during each part of a period 
of activity (Stier, 1930). We have also been unable, as yet, to find out 
how the “‘intensity”’ of muscular activity exhibited during single activity 
periods is related to temperature. Probably the frequency of leg move- 
ments and of body movements (cf. Stier, 1930), and the tonus of the skele- 
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FIGURE 1 





Calculated total “spontaneous” activity per hour of two-day old albino mice 
as a function of temperature. Ordinate: the product of the average duration 
of single periods of activity times the number of periods per hour. Abscissa: 
body temperature, °C. The observations are plotted as averages of data from 
several litter mates; each spoke on the circle refers to the number of animals 
included in each average. 














tal muscles, increases with temperature. The additional heat contributed 
by these sources should likewise change with temperature. The heat 
production, however, as indicated by thermopile measurements of the 
body temperatures of young albino mice (Pincus, Sterne and Enzmann, 
1933) of the same strain used by Stier (1930), increases as the body tem- 
perature is raised from 15° to 25°C.; and decreases progressively from 25° 
to 35°C. The diminution of heat production beyond 25° can be possibly 
explained by the decrease in duration of activity periods with temperature 
and by the peculiar relation of the frequency of occurrence of activity 
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periods to temperature. The number of activity periods per hour in- 
creases as the body temperature is progressively raised from 16° to 24°C. 
In the neighborhood of 25° a sudden diminution in the number of periods 
occurs (34 at 24°; 20 at 26°). Thereafter, the frequency again increases 
as the body temperature is increased (cf. Stier, 1930). It may be found in 
subsequent experiments that the ‘intensity’ of muscular activity follows 
a similar relation to temperature; or that it steadily increases from 15 to 
35. The net effect of the additional heat contributed by the “‘intensity”’ 
of spontaneous movements, no matter what the type of temperature rela- 
tionship, appears to be of minor importance in accounting for the heat pro- 
duction above 25°. It has already been found (Pincus, Sterne and Enz- 
mann, 1933) that the difference between internal and external tempera- 
ture rises to a maximum at 25° and then falls off at higher body tempera- 
tures. Below 25° the heat contributed by the intensity of movements and 
by other heat producing processes (e.g., cellular metabolism) may account 
for the difference in location of the maxima in the graphs by Pincus, Sterne 
and Enzmann (1933) at about 25° and the peak in figure 1 at about 22°. 
A quantity referring to the intensity of muscular activity should be there- 
fore included in the product—duration of activity periods times frequency 
of occurrence— in order to account entirely for the rate of heat production 
by the active muscles. A measure of this quantity will be sought in subse- 
quent experiments. 

The data of Pfliiger (Starling, 1930) on the O2 consumption of rabbits 
in the range of internal temperatures from 22° to 38.8°C. yield graphs 
quite similar to figure 1, the peak occurring at 37.5°C. In Martin’s ex- 
periments (cf. Starling, 1930) on the CO, production of rabbits and mea- 
surements of the O, consumption of mice by Benedict and Fox (1933), only 
the external temperature was recorded and so the data cannot be used for 
our present analysis. More experiments in which spontaneous activity 
is measured over a wide range of internal temperatures are needed to com- 
plete the study of temperature regulation. 

Heat regulation by the adult mammal can be upset by injuring the 
corpus striatum (Barbour, 1921). As yet it is not known whether the 
“center” controlling the periodicity of movements is identical with the 
heat regulatory center, or whether it is closely associated with it. Tracy 
(1926) has pointed out that the center controlling ‘‘spontaneous’’ move- 
ments may be closely associated with the respiratory center. 

On the basis of the information contained in figure 1 and the findings of 
Pincus, Sterne and Enzmann (1933), it is suggested that in the range of 
body temperatures above 25° (approximately above 23.7°C., air tempera- 
ture for mice of these ages) the two-day old mouse already exhibits a kind 
of regulation of its body temperature. If, for instance, a young albino 
mouse kept in a nest at, say, 30° is moved out into a room at 25° there 
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occurs (after one hour) a 13.5 per cent increase in the product: duration 
times number of activity periods per hour; if, however, an animal is 
taken from the same nest to an environmental temperature of 32°C. there 
occurs a 32.7 per cent decrease in the total amount of muscular movements 
(as measured by the units plotted in figure 1). There is thus a tendency to 
produce more heat at temperatures below that of the nest and less heat at 
temperatures higher than that of the nest. The extra heat produced at 
lower temperatures is of course rather rapidly lost in these young furless 
individuals. However, as the animal becomes older, thicker layers of in- 
sulating materials are acquired, heat is lost less rapidly when the tempera- 
ture of the environment falls, and consequently the fluctuations of internal 
temperature become narrowed in their extent. So in conjunction with the 
other heat-producing processes and the heat-conserving and -dissipating 
mechanisms, the body temperature of the adult animal is held constant at 
some fixed point. 

Experiments now in progress with older mice of the same inbred strain 
(the Bagg albinos) used in the above experiments (Stier, 1930) will test 
whether a similar mechanism controlling the heat production at different 
body temperatures is retained in adult animals. If the mechanism out- 
lined above for young mice exists in adult animals one can then understand 
why a decrease in the body temperature of the so-called homeothermic 
mouse brings about an increase in the amount of movement of the animal 
and so an increase in heat production; and an increase in the external tem- 
perature brings about a decrease in the amount of bodily activity. 

It is important in studying the nature of temperature regulation in 
mammals to ascertain the magnitude of the internal temperature main- 
tained over long periods of time and the actual heat production at many 
steps in a wide range of body temperatures, studying at the same time the 
development of the mechanisms which dissipate and retain the heat con- 
tributed by muscular activities and by other heat-producing processes. 
This procedure is being followed in experiments by Dr. G. Pincus and the 
writer upon several strains of inbred mice at various ages after birth. 

Summary.—The total amount of “spontaneous” muscular activity of 
two-day old albino mice increases as the animal’s internal temperature is 
raised from 16° to 22°C. and then steadily decreases at successive higher 
temperatures. Similar graphs have been obtained for the difference be- 
tween external and internal temperatures of young albino mice of the 
same strain in the range of body temperatures, 15° to 35°C. (cf. Pincus, 
Sterne and Enzmann, 1933). The maximum difference occurs at about 25°C. 

It is suggested that there is exhibited a kind of regulation of tempera- 
ture by the young mouse in the range of internal temperatures above 25°. 
In this upper temperature range a decrease in air temperature results in an 
increase in the amount of muscular activity and an increase in the heat 
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production; an increase in the temperature of the air is followed by a de- 
crease in the amount of muscular activity and a decrease in the production 
of heat. 

* The work in this paper was supported by a grant from the Milton Fund for Research, 
Harvard University. 
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THE DEVELOPMENT OF TEMPERATURE REGULATION 
IN THE MOUSE* 


By G. Pincus, G. DEROO STERNE AND E. ENZMANN 
LABORATORY OF GENERAL PHYSIOLOGY, HARVARD UNIVERSITY 
Communicated June 9, 1933 


In experiments dealing with the temperature relations of certain vital 
processes in young mice (Stier and Pincus, 1927-28; Stier, 1930; Pincus, 
1931) it was noted that the body temperature of even one-day old mice 
was ordinarily not quite that of the environment, despite their complete 
lack of fur and their small size, but, over the range 15° to 35°C., some- 
what above that of the environment. The difference between external 
and internal temperatures seemed to be roughly a function of the age of 
the animal, such that the older the mouse the greater the difference. It 
was thought, therefore, that by measuring this difference in a series of 
mice from birth to maturity a quantitative measure might thereby be had 
of the relation of the ability for temperature regulation to age. The simple 
procedure would seem to be to employ one or more standard constant 
external temperatures and to ascertain the difference between these tem- 
peratures and the internal temperature of mice of various ages, and then 
by plotting this difference against age a presumably rising curve would be 
had indicative of the nature of the age-regulatory relation (cf. Baldwin and 
Kendeigh, 1932, p. 108). Such a procedure was found to be quite inade- 
quate, however, because of two definite difficulties. 

The first difficulty involves the exact estimation of the internal tempera- 
ture of a partially homeothermic mouse kept at a constant external tem- 
perature. The technique employed, described in detail by Stier and Pincus 
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(1927-28), involves the measurement of rectal temperatures by the use 
of a thermopile.' The rectal temperature of a young mouse taken from 
the nest and placed in a thermostat kept at a temperature of, say, 16°C. 
decreases fairly regularly with time. The shapes of the cooling curves of 
live mice differ markedly from those for dead mice of the same age. This 
is illustrated in figure 1, where the cooling curves of a live and of a dead 
mouse of approximately the same age and weight, and kept at approxi- 
mately the same temperature, are plotted logarithmically. The dead 
mouse exhibits a typical Newtonian cooling curve, as demonstrated by the 
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FIGURE 1 


Comparison of the cooling curves of live and dead mice. 
Ordinate: logarithm of the difference between internal 
and external temperatures. Abscissa: time in minutes. 
Open circles = 5-day old mouse at 15°32 + 0°02. Closed 
circles, 6-day old dead mouse at 15°95. Note difference 
in rate of cooling and in the shape of the cooling curves. 
See text. 


linear relation between elapsed time and the logarithm of difference be- 
tween external and internal temperature. For purposes of comparison 
the cooling curve of the dead mouse has been extrapolated to intersect the 
curve for the live mouse at a point where the difference between external 
and internal temperatures in both live and dead mice is 2°52C. (0.4014 in 
logarithmic units), and the comparative rates of cooling may be adjudged 
by the slopes of the respective lines. The rate of cooling is much slower 
in the live mouse, nor does it follow Newton’s law. Itisas if some internal 
mechanism had put a “‘brake”’ upon the cooling process. It is conceivable 











VoL. 19, 19833 PHYSIOLOGY: PINCUS, STERNE AND ENZMANN 731 


that a measurement of this “‘braking’’ would give a measurement of the 
way the temperature-regulating mechanism responds to temperature 
changes. 

When a dead mouse is warmed, the rate of attainment of the higher tem- 
perature is also logarithmic. When a live mouse is placed at a higher tem- 
perature, however, quite a different situation is observed. The rate at 
which the live mouse attains the higher temperature is now much greater 
than the rate of warming of a dead mouse of the same age and weight, and 
the live mouse, within a certain temperature range at least, extremely soon 
attains a body temperature above that of the environment. So rapidly 
does the body temperature of a live mouse rise when the animal is placed 
at a higher temperature that it is often difficult to obtain readings for the 
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FIGURE 2 
Comparison of the warming curves of live and dead mice. Ordinate and 
abscissa as in figure 2. Open circles = 2-day old live mouse at 32°81 + 
0°02. Closed circles = 2-day old dead mouse at 30°70 = 0°02. Note 
early rapid attainment of external temperature by live mouse and subsequent 
rise of internal temperature. Observe also periodic fluctuations in internal 
temperature of live mouse. 





rate of warming. The course of events is illustrated in figure 2. The first 
five points on the curve for the live mouse represent the period when the 
body temperature of the animal is below that of the environment. There- 
after the body temperature rises above that of the environment, but not 
continuously or smoothly. The peaks and valleys in the curve represent 
periods of rise and fall in the body temperature. We continually encounter 
such evidences of rhythmic change in internal temperature, and our esti- 
mation of the “normal” level of body temperature is made difficult because 
of the presence of these rhythms. Such rhythms are, of course, never 
encountered in dead animals. 

If we discount the spontaneous rhythmic changes, and proceed to esti- 
mate the ‘“‘normal’’ body temperature of a series of young mice at various 
constant external temperatures, our second and most serious difficulty 
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is encountered. This is revealed when we ascertain the difference between 
external and internal temperature as a function of temperature. Such data 
are plotted in figure 3 for mice of three different ages. It is apparent that 
the curves relating temperature difference to temperature are cusped in 
each case. Over the range 25°+ to 35°+ this difference decreases (cf. 
Stier, 1933). Over the lower temperature range we observe the ordinary 
relation of chemical processes to temperature (the temperature coefficient 
is positive and high); over the higher range we have a negative tempera- 
ture coefficient. The latter relation is what one would expect in an animal 
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FIGURE 3 


Calculated ‘‘normal”’ difference between internal and external tempera- 
tures for live mice of various ages as a function of temperature. Ordinate: 
difference between internal and external temperatures in °C. Abscissa: in- 
ternal temperatures in °C. Closed circles = 2-day old mouse. Half circles 
= 8-day old mouse. Opencircles = 10-day old mouse. See text. 


capable of temperature regulation, i.e., the higher the temperature the less 
the difference between external and internal temperature. It will be noted, 
too, that the slope of the curves is apparently steeper over the lower tem- 
perature range. The rising curve occurring over the lower temperature 
range we have designated the “‘metabolic’’ response, the second curve the 
“regulatory” response. Stier (1930) has noted a similar critical tempera- 
ture at 25°C. concerned with certain aspects of spontaneous activity of one- 
day old mice. He employed mice of the same inbred strain (Bagg-Little 
albinos). 

It becomes difficult to ascertain the development of the temperature- 
regulating mechanism with age, because the most elementary comparison 
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would require knowing whether the constant external temperature used 
results in the attainment of an internal level proper to the ‘‘metabolic’’ 
phase or to the “regulatory” phase. Whether a given level falls into one or 

the other category is apparently conditioned not merely by the external 
’ temperature but also by the internal temperature. Furthermore, we have 
some evidence that in older animals the ‘‘regulatory”’ response may occur 
at temperatures below 25°C. Experiments now in progress are being con- 
ducted to ascertain exactly the nature of these two phases. 

The data of figure 3 roughly illustrate that with increasing age the differ- 
ence between external and internal temperature increases in both the 
“metabolic” and the “regulatory’’ phases. The oldest mouse examined 
in the present series was 24 days old, and even at that age perfect regula- 
tion, i.e., the ability to maintain the body temperature at a constant high 
level (e.g., 36°8C.) at any outer temperature, was not observed. The 
body temperature of adult mice of the same strain is normally 36°8+, 
but it is possible to lower this temperature by various means. Details of 
experiments with adult mice will be given elsewhere. 

Summary.—tThe rate of cooling of living mice aged 2 to 24 days is always 
slower than that of dead mice of the same age and weight. Conversely, 
the rate of warming of live mice is always faster than that of dead mice of 
the same age and weight. At any constant external temperature over the 
range 15° to 35°C. the rectal temperature of living mice is always 
somewhat above the outer temperature. The equilibrium difference be- 
tween internal and external temperature is not constant, but fluctuates 
rhythmically about a mean. When the difference between external and 
internal temperature, for mice of the ages used, is plotted as a function of 
temperature over the range 15° to 35°C., it is found that this difference 
increases to a maximum at 25°C.+ and then falls. The ascending limb 
of the curve has been designated the “‘metabolic’’ phase, the descending 
limb the “‘regulatory” phase. No exact description of the relation of the 
development of temperature-regulating ability to age can be given without 
the resolution of the réle of these two phases at all ages. 


* The work summarized in this paper was in part supported by the Milton Research 
Fund, Harvard University. 

1 In the experiments described in this paper the mice used were albinos of the 42nd to 
45th brother-sister generation of the Bagg-Little strain. 
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THE INFLUENCE OF USE OR FREQUENCY OF OCCURRENCE 
UPON THE STRENGTH OF MENTAL CONNECTIONS 


By EDWARD L. THORNDIKE 
INSTITUTE OF EDUCATIONAL RESEARCH, TEACHERS COLLEGE 


Read before the Academy, Tuesday, April 25, 1933 


Until recently the doctrine in mental dynamics which has been accepted 
most widely and used oftenest as an explanation of learning is the doctrine 
that the mere occurrence of a connection strengthens it. It has seemed 
almost axiomatic, but almost all the evidence showing that use strengthens 
connections is contaminated by the presence of satisfying after-effects, 
so that it is possible to maintain with some reason that the occurrence or 
use of a connection does nothing but give an opportunity for other forces 
to act upon it. The Gestaltists have also challenged the doctrine on 
general grounds. 

In The Fundamentals of Learning evidence was reported showing that 
the mere occurrence of a connection does strengthen it, in the shape of 
the strengthening found where after-effects of any sort were at a minimum, 
and also in the shape of the strengthening of connections which were 
punished. The latter evidence was not treated fully or critically in that 
volume. It has been extended by Lorge and the writer, and will be the 
topic of the present report. 

In experiments with human multiple-choice learning with reward and 
punishment, a connection which occurs and is punished often is strength- 
ened more by the occurring than it is weakened by the punishment. 
Evidence of this has been presented by Thorndike and others [’32, p. 108 ff. 
and p. 276 ff.] and by Lorge [’33]. We shall determine the extent to which 
this happens and the amount of the balance of strengthening, and observe 
the characteristics of individuals in respect of it. Experiments made to 
measure the “‘spread’’ or “‘scatter’’ of the influence of rewards and the 
influence of the relevance of rewards provide the necessary data.! 

Consider first the facts from an experiment (Q) in learning one response 
out of six for each of 960 situations, the response being the numbers 1 to 6 
and the situations 960 words, divided into 24 sets of 40. A word was read 
by the experimenter; the subject said 1 or 2 or 3 or 4 or 5 or 6; the experi- 
menter said “Right” or ‘“Wrong,”’ and then read the next word, and so on 
through a series of 40, and through the same series four times more. 
Each “Right” signified a money payment which was of real importance 
to the learner. There were 24 subjects. 

The influence of one rewarded occurrence in Experiment Q is such 
that the connection which occurs and is rewarded occurs again in the 
following trial or round in 4029 out of 9751 cases, or in 43.6 per cent.? 
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The influence of one punished occurrence is such that the connection which 
occurs and is punished occurs in the-.following trial or round in 9868 out 
of 43,915 cases, or in 22.5 per cent. 

If we assume that, apart from the occurrences, rewards and punish- 
ments, the subject would at any occasion be as likely to choose any one 
number as any other, the rewarded occurrence has increased the proba- 
bility of reoccurrence by 0.269 (from 0.167 to 0.436) and the punished 
occurrence has increased it by 0.058 (from 0.167 to 0.225). 

Some of the latter increase is, however, probably due to the ‘‘spread’”’ 
or ‘‘scatter” action of rewards, and some of it may be due to favoritism in 
the choice of numbers. We will consider the latter first. 

The percentage of repetitions which there would have been if the punish- 
ment and occurrence had had a net influence of zero would be 100 divided 
by the number of different responses (here 6), if the subject was as likely 
to select one response as another. But if he favored certain responses 
it would not. In any trial (call it II) the probability of repetition of a 
response made in the preceding trial (call it I) apart from all influence 
from the learning itself, is al all + bI bII + cl clI....... nI nII where 
al, bl, cl, etc., are the probabilities of occurrence in trial I and all, bII, 
cII, etc., are the probabilities of occurrence in trial II. Favoritism for 
certain numbers can lower the probability of repetition to 0 or raise it to 
1.00. For example, in a case where the subject must respond by any 
number from 1 to 6, we could have such cases as A, B, C and D below 
(the occurrences total 1000 in each trial). 


Case A: 
Trial I I, 400:.- 2, 300; <3, 200; 4; FOO): 6, 05° 6; 0 
Trial II 1, 0;.. 2, OF > 3; ie 0; 5, 400; 6, 600 
Probability of repetition, 0.000. 

Case B: 
Trial I 1, 1000; 2; 0: 3; Of 4, 0: - 5, 6; -6; 0. 
Trial II 1, O00; = 2, ee Oe > Ay OF Os: ~6, 0. 
Probability of repetition, 1.00. 

Case C: 
Trial I 4 @OOs «De 8BOe .3¢. 200.4, -1OBy- B, 0;::..6, 0. 
Trial II 1; 40022; .300;:: 3, 2003: 45. 2003. .6, @:. 6, 0. 
Probability of repetition, 0.300. 

Case D: 
Trial I Tt, . 366s. 2 260; 3, - BOG; 4, 1003: Sd, 50; 6, 50. 
Trial II 1. 260s Be Os a ee a 100.8; COs 6; 6. 


Probability of repetition, 0.238. 


Favoritism of a similar sort in both trials will raise the probability of 
repetition above 0.1667. Reversal of favoritism in the second trial will 
lower it below 0.1667. - 
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We therefore proceed to investigate the favoritism, if any, shown 
by each of the twenty-four subjects.in choosing among 1, 2, 3, 4, 5 and 6 
as responses. The responses of each are tabulated in successive sets of 40 
corresponding to the successive trials of the various series of words, then 
in successive sets of 200, each corresponding to the 200 responses for the 
five trials of one series, and finally by successive sets of 400.* First, 
each record is inspected to detect any tendencies to reverse favoritism in 
any two consecutive trials. No such were found. 

Secondly, each record is inspected to detect any systematic shifts 
in favoritism (such as a steady change from favoring 1 and 2 much to 
favoring them little, or from favoring 2, 3, 4 and 5 to favoring them little 
or vice versa). There were few such shifts. 

Thirdly, each record is inspected to detect any clear cases of short- 
lived changes in favoritism. It is, of course, hard to draw the line be- 
tween such and chance variations in choices. 

If there is no evidence of systematic shifts, or of genuine short-lived 
tendencies, the average influence of favoritism operating at any one time 
is estimated as follows: a? + b? + c? + d? + e? + f? / n? is computed for 
each successive 400 (or for so many 400’s as we have) and the average of 
the determinations is taken. This permits undetected genuine short-lived 
tendencies to favoritism to count somewhat, but still does not interpret 
chance variations in choice as genuine favoritisms. If there are any special 
shifts, the 400’s are so chosen as not to swamp the influence of the favorit- 
ism, and are checked by using sets of 200 each. 

As a result of much labor expended in the ways stated, a probability 
of repetition of the same number in the responses to the same word in two 
successive trials has been estimated for each of the 24 subjects.‘ They 
range from 0.170 to 0.220, twenty being between 0.174 and 0.205. These 
are shown in table 1 together with the excess frequencies of repetition for 
punished occurrences at various degrees of remoteness from any rewarded 
connection. 

Even at a remoteness of five or more steps from any rewarded connec- 
tion the frequency of repetition in the next trial for one punished occurrence 
of a connection is clearly in excess of what would be expected by chance 
plus number favoritism. It is + in 19 cases out of 24, and the median 
(16 permilles) is about three times its probable error. A punished connec- 
tion gains more strength by occurring once than it loses by being punished 
once. 

Why this should be so is an interesting problem. Some facts and 
possibilities are as follows: (1) Lorge and Thorndike [’33] have shown 
that in man a connection left ambiguous and neutral by having no announce- 
ment or reward or punishment of any sort attached to it gains still more. 
(2) The net gains shown in table 1 may then be reductions from still 
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larger gains due to the occurrences pure and simple. (Lorge does not, 
however, regard this as probable.) (3) The strengthening by the mere 
occurrence may be reduced not at all by the punishment in these experi- 


TABLE 1 


THE STRENGTHENING DvE TO ONE (First) OCCURRENCE OF A PUNISHED CONNECTION. 
EXPERIMENT Q. 


EXCESS OVER THE EXPECTED FREQUENCY FOR THE ACTUAL FREQUENCIES OF REPETITION 
IN THE FOLLOWING TRIAL FOR FIRST OCCURRENCES OF PUNISHED CONNECTIONS. 
PERMILLES 


a a 8 g n 
Be ch ee 
5} z % ao a5 Ms Be Daz 
+5 ae $9 #8 880s He, «68 oSe 6888 
a a BE ha Ba ate onS SEES Op 
OD 22 yaa g BOBO fy Rous Aaa 
CE A 
mks ra ae gs OoOSs ae aeee Ores 
i} a2 B & x az BORZ 9 
44 oe 11 et ee ee 
Anz 28 Amz Zw z ae ez RRO Raoo mem Zz 
1 2 3 4 5 6 7 8 
1 200 109 145 107 157 107 119 091 
2 174 033 ~=008 038 063 053 039 026 
3 170 027 071 057 008 009 024 002 
4 189 —003 O71 005 081 —051 —Ol11 005 
5 185 010 041 017 093 019 011 003 
6 196 007 057 013 035 025 =—023 024 
iT 210 030 = 167 036 —130 026 030 —017 
8 199 031 126 023 —045 025 017 050 
9 175 064 075 075 044 081 082 019 
10 185 071 081 073 194 090 070 077 
11 174 020 072 050 039 «4-006 -—031 —036 
12 208 035 109 045 023 049 026 —009 
13 180 039 #8146 044 106 —005 027 030 
14 180 025 107 035 ° 120 037 023 002 
15 180 050 060 039 —013 046 040 069 
16 174 045 022 081 159 023 038 042 
17 177 004 058 055 109 -001 -—027 -—020 
18 182 042 032 079 100 012 057 004 
19 205 014 032 022 045 007 028 —019 
20 189 111 061 094 144 078 144 107 
21 183 048 117 075 140 039 035 011 
22 190 052 174 081 041 013 038 042 
23 220 033 ~=—(018 014 165 049 044 042 
24 180 00s «sill 000 062 —004 —008 013 
Median 184 032 071.5 044.5 072 025 029 016 
Average 188 038 082 048 073 030 033 023 


ments, so that the 016 is the strengthening due to one occurrence pure 
and simple. (4) The spread or scatter effect may extend out further 
than past information has led us to expect. (5) There may be a strengthen- 
ing due to the satisfyingness of making any choice and response, which 
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is reduced by the annoyingness of finding it to be wrong. (6) Certain 
words may more or less persistently suggest certain numbers to certain 
subjects by reason of their meanings or associations. We shall return to 
this problem after reporting the facts for other experiments. 

Consider now the facts for Experiment HI and Ha (described in ‘‘An 
Experimental Study of Rewards,” Section II) in which the learning was to 
connect a number from 1 to 10 with each of 480 words (H1) and 160 words 
(Ha). Two readings of each series of 40 words and the right numbers to 
the subject before his first trial with that series put the chance for a right 
choice in trial 1 well above 0.10, but that does not concern us here. The 
percentage of repetition in the trial following a first occurrence that was 
punished was far in excess of 0.10, being over 0.24, and was 0.21 when the 
punished connection was five steps or more removed from any reward. 
We have to discover how much of this excess was due to number favoritism. 

Experiment HI really comprises ten different sub-experiments, each 
involving the responses to 40 words in each of four trials. What we really 
need to know, in order to allow perfectly for the influence of favoritism 
for numbers, is the status of favoritism before each response of each of the 
sub-experiments, or, since that is impracticable, at the beginning of each 
trial of each sub-experiment. That is, we should like to know the nature 
and amount of favoritism at any stage in the learning, regardless of how 
it was caused, so as to allow for its action point by point, or at least trial 
by trial, in the course of the learning. But even the favoritism trial by 
trial would be exceedingly difficult to estimate. The size of the samples 
(40) is too small for reliable measurements, and we should have to work 
out correction factors. Until the need for allowance trial by trial becomes 
acute nobody will, or should, undertake the labor. We may, however, 
reasonably undertake the labor of measuring the favoritism in the first 
two sub-experiments, the second two sub-experiments, etc., since that 
will give us N’s of 320, which should not increase the obtained favoritisms 
greatly above the real tendencies which they sample. We have done this 
with the facts for the first 960 choices of Experiment Hl, and the 640 
choices of Experiment Ha, with the results shown in table 2. 

If the tendency of a person in the 40 choices of any trial in Hl was to 
favor numbers as he did in the entire 960, the probability that, apart from 
any other influences than his number favoritism, he would give the same 
number for any given word in any two successive trials would be 0.113 
for Bra, 0.104 for Bre and 0.117, 0.105, 0.110, 0.103, 0.115, 0.120, 0.111 
and 0.106 for the other subjects in order. 

In so far as short-lived habits or systematic shifts of favoritism existed 
these figures are too low. They are not much too low. The division 
of the 960 into three 320’s in order of occurrence and the computation of 
the average of the same probabilities from the three 320’s gives results 
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TABLE 2 
FAVORITISM FOR NUMBERS. THE FREQUENCIES WITH WHICH BRA, BRE, ETC., USED 
1, 2,3, ...10 1n SUB-EXPERIMENTS 1 AND 2, 3 AND 4, AND 5 AND 6 OF Exp. HL AND IN 
1 AND 4 AND 7 AND 10 or Exp. HA 
BRA BRE BL 
12 3,4 56 1,4 7,10 iz 34° 668° 447,10 1,2 3,4 56 1,4 7,10 
1 48 35 31 27 25 27 47 38 33 33 31 24 16 27 15 
2 33 47 25 36 37 23 38 28 10 33 35 51 31 46 31 
3 37 40 35 26 30 47 32 31 40 30 42 41 40 32 32 
4 60 48 56 54 42 39 41 50 32 36 16 8 27 35 46 
5 18 30 28 33 32 35 19 21 29 36 36 22 48 53 55 
6 37 41 37 37 49 17 44 41 38 24 33 41 28 42 39 
7 30 19 33 17 25 38 34 40 25 59 44 22 37 20 9 
8 14 16 21 22 29 34 34 18 55 33 35 33 32 26 19 
9 43 30 29 31 32 35 14 25 14 22 31 34 19 38 49 
10 0 14 24 37 19 25 17.278. .0 14 17 44 42 1 25 
A. BU. L. M. 
12 3,4 56 1,4 7,10 12°34 66 1,4°7,10 13 $4: 66: 24°7;20 
1 23 28 31 39 28 7 11 26 24 36 42 42 26 22 38 
2 65 40 51 47 56 28 32 42 31 26 14 22 41 28 50 
3 41 39 40 37 36 53 38 41 53 34 46 16 28 28 27 
4 50 54 46 45 54 46 48 48 30 38 30 48 50 61 25 
5 is 12. .8 36-36 37 29 26 61 41 41 27 18 28 48 
6 8 21 26 22 27 55 40 30 34 51 26 29 35 52 30 
7 29 32 27 23 36 14 31 22 29 17 35 34 32 42 48 
8 39 42 49 38 21 28 21 40 21 24 oT 3dr Bt 16: St 
9 23 19 27 26 36 34 57 25 30 24 23 31 20 21 16 
10 29 33 18 10 18 1816.25: -¢ 20 36 34 48 22 22 
PEI RO. E. T. 
12 54: 66.1% 750 12 3,4 56 1,4 7,10 12 3,4 56 1,4 7,10 
1 29 30 30 14 19 29 41 46 29 39 30 18 28 22 25 
2 39 30 48 11 22 58 57 51 64 46 27 32 34 31 25 
3 46 38 25 31 33 29 13 18 15 9Q 45 27 42 39 46 
4 22 438 41 56 55 21 17:10. 2B. 14 39 51 44 31 25 
5 63 56 61 30 63 37 35 45 51 54 29 36 45 34 35 
6 41 46 31 33 31 64 27 48 40 34 53 57 52 55 33 
7 28 22 17 53 48 23 54 48 39 62 26 19 20 20 48 
8 11 138 25 45 25 2. a1 TT 19S 20 20 23 28 37 
9 14 14 26 33 10 1 Ee ae 26 37 17 45 36 
10 31.8: 16° 14. 14 24 40 46 19 41 25 23 20 15 15 
R. FT. 
19 84 86 14 7,10 
1 30 30 36 16 20 
2 23 29 47 44 42 
3 43 39 338 45 41 
4 40 42 32 35 38 
5 16 22 31 34 33 
6 47 41 36 38 39 
7 26 15 18 38 27 
8 39 46 38 19 25 
9 34 32 32 42 32 
10 22 24 17 9 238 

















PSYCHOLOGY: E. L. THORNDIKE Proc. N. A. S. 





740 


only 0.0044 higher. Inspection of the choices by successive 160’s shows 
littie or nothing striking in the way of short-lived habits or systematic 
shifts. Moreover, these subjects in the similar experiment (Ha) over a 
week before Experiment Hl, showed much the same sorts of favoritisms 
as they showed in Hl. 

We have computed the sums of the differences between the number 
of choices of the same number in the various 320’s. For neighboring 320’s 
the sums are 922, 856 and 998 (averaging 925). For the first 320 and 
third 320 of the 960, the sum is 1029. For any 320 of Hl and any 320 of 
Ha the sums are 1281, 1162, 1171, 1109, 1113 and 997 (averaging 1139). 
We shall therefore not be far wrong if we use the average of the separate 
probabilities for the five 320’s of Experiment Ha and Hl as the base from 
which to measure the influence of one punished occurrence in Experiments 
Ha and Hl. 

Table 3 presents, for each subject in Experiment Ha and Hl, observed 
frequencies of repetition in the following trial of a punished connection 
which at its first occurrence was five or more steps removed from any 
rewarded connection, along with this base (and also the five probabilities 
of which it is the average). There is in all subjects an excess of repetition 
above that explicable by favoritism for numbers. The average excess 
(092) is 6!/: times its probable error. The individual differences are large, 


as they were in Experiment Q. 

TABLE 3 
THE OBSERVED FREQUENCY OF REPETITION IN THE FOLLOWING TRIAL OF A PUNISHED 
CONNECTION FIVE OR MORE STEPS REMOVED FROM ANY REWARDED CONNECTION AT 
Irs First OcCURRENCE (OBS.), AND THE APPROXIMATE FREQUENCY OF REPETITION 
TO BE EXPECTED FROM NUMBER FAVORITISM ALONE (FAV.), IN PERMILLES. ALSO THE 
PROBABILITIES FOR EacH COLUMN OF TABLE 2. THE DATA ARE FROM Exp. HL AND HA 


EXCESS OF 


OBS. OVER EXP. Hl EXP. Ha 

OBS. FAV. FAV. 1,2 3,4 5,6 1,4 7,10 
Bra 132 113 019 127 114 109 109 107 
Bre 124 111 013 107 112 109 116 112 
Bl 298 114 184 107 115 109 119 121 
A. Bu 171 118 053 126 114 119 112 119 
L 167 116 051 123 119 106 121 109 
M 298 112 186 109 108 109 119 113 
P 235 121 114 122 116 118 122 128 
R 277 128 149 124 125 133 122 135 
E. T. 237 112 125 109 116 115 113 108 
RF. 135 109 026 109 109 107 114 106 
Median 084 
Average 092 (P.E. = 014) 


In the case of Experiment F (in which the choice was of a number from 
1 to 10) the individual records of actual repetitions of punished connections 
remote from a reward were from too few cases to deserve computation. 
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The totals for the group of twenty subjects were: 605 first occurrences 
five or more steps removed from any reward, of which 74 or 12.2% were 
repeated in the following trial. The probability of a repetition apart from 
occurrences, rewards and punishments computed from the frequency 
of choices in the first 4 trials varies from 0.109 to 0.141, averaging 0.120. 
When computed from the last two trials, it varies from 0.102 to 0.231, 
averaging 0.140. The two medians are 0.1175 and 0.136. The frequency 
of repetition by chance plus number favoritism may then be set as approxi- 
mately 0.130. There thus appears to have been no balance of strengthening 
effect from an occurrence of a punished connection in this experiment. 

For Experiment P (described in Lorge and Thorndike, p. 374) table 4 
gives the facts concerning the expected frequency of repetition for each 
individual by chance plus his number favoritisms and concerning,the actual 
repetition for punished connections 5 or more steps removed from any 
reward.> The number of observations for an individual is small (about 
85 in Part A and about 45 in Part B) so that the per cents in column 6 
and the differences in column 7 are not individually reliable. But the 


TABLE 4 


THE STRENGTHENING DUE TO ONE (FirST) OCCURRENCE OF A PUNISHED CONNECTION 
IN EXPERIMENT P. IN PERMILLES. 


EXPECTED FREQUENCY OF REPETITION BY CHANCE ACTUAL 
PLUS NUMBER FAVORITISM REPETITION OF 
PUNISHED ACTUAL MINUS 
TWO WITH CONNECTIONS EXPECTED BY 
MOST ESTIMATE 5 OR MORE CHANCE PLUS 
FIRST LAST FIRST FAVORIT- FOR EX- STEPS FROM NUMBER 
FOUR FOUR Two ISM PERIMENT A REWARDED FAVORITISM 
n = 512 n= 6508 nm» = 256 nm = 252 ASA WHOLE CONNECTION 5 
PART A 1 2 3 4 5 6 7 
B 290 258 286 304 280 353 073 
H 252 254 250 359 265 337 072 
L 264 279 265 282 270 333 063 
M 256 256 258 275 260 567 307 
Se 260 251 254 271 260 402 142 
Sh 260 257 263 264 260 216 —044 
< 255 251 255 271 260 628 368 
Ca. W. 286 286 295 304 290 254 — 036 
Ch. W. 261 255 257 266 260 255 —005 
PART B n= 352 n=352 n=176 n» = 176 
Bel 254 255 263 287 260 452 192 
D 251 254 252 261 255 209 —046 
La 258 250 266 295 260 429 169 
Al 304 442 330 487 370 400 030 
Fl 257 261 259 282 260 375 115 
Ri 256 255 257 262 255 295 040 
Ro 287 252 286 296 275 408 133 
5.0: 252 259 256 302 265 500 235 
7%, 259 259 264 267 260 385 125 


264 305 305 270 224 —046 
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general results (median excess over what chance and number favoritism 
would give, 073) is about four times its probable error. 

A similar treatment of the records in Experiment M results in the facts 
of table 5. The measurements for any one individual are from too few 
cases to be reliable, but the average (0.037 with a P. E. of 0.016) indicates 
that a punished occurrence has a strengthening influence upon the connec- 
tion in question in this experiment as in Experiments Q, H and P. 


TABLE 5 
THE STRENGTHENING DUE TO ONE PUNISHED OCCURRENCE IN Exp. M, COMPUTED BY 
Pact—Pavy. Pa, Ps, Pc AND Pp ARE THE PROBABILITIES OF A REOCURRENCE OF THE SAME 
RESPONSE IN Two CONSECUTIVE TRIALS IN SUB-EXPERIMENTS A, B, C anD D, APART 
FROM ANY INFLUENCE OF OCCURRENCES, REWARD AND PUNISHMENTS. Pay IS THE 
AVERAGE OF THESE Four PROBABILITIES, Pact IS THE ACTUAL FREQUENCY OF REPE- 
TITIONS IN THE FOLLOWING TRIAL OF PUNISHED CONNECTIONS FIVE OR More STEPS 
REMOVED FROM ANY REWARD (ONLY First OCCURRENCES BEING USED). 1 = THE 
NUMBER OF CASES USED IN DETERMINING Pact 


STRENGTHENING 
DUE TO ONE 
PUNISHED FIRST 


SUBJECT Pa PB Po Pp Pav Pact % OCCURRENCE 
A 0.171 0.226 0.307 0.218 0.231 0.284 67 0.053 
Bi O.170.. 0.170 30.170. 0.187 0.177 . 0.221 68 0.049 
Be 0.195 0.257 0.178 0.180 0.203 0.085 82 -—0.118 
B; 0.245 0.208 0.235 0.265 0.2388 0.346 78 0.108 
L 0.173 0.183 0.192 0.172 0.180 0.176 85 0.004 
M 0.181 0.198 0.180 0.174 0.183 0.226 62 0.043 
O 0.353 0.350 0.454 0.297 0.364 0.477 Ad 0.113 
N 0.186 0.209 0.182 0.189 0.192 0.240 129 0.048 


For the 20 subjects of Experiment J the facts were as shown in table 6.° 
There is an excess of repetition due to one punished occurrence 5 or more 
steps removed from any rewarded connection in every individual. The 
median is 98 permilles, which is about 10 times its probable error. 

For the 8 subjects of Experiment L, the median excess was 46 permilles.® 

In experiments upon the influence of varying amounts of money reward,’ 
Rock found the percentage of repetition of punished connections to be 
33.5 where pure chance would give 25.0 and where chance plus position 
favoritism would give (according to the facts in a similar experiment) 
about 27.0. Even after a liberal allowance for the influence of the spread 
or scatter, there will remain a substantial amount of strengthening due to 
the mere occurrence of the connection. 

With the single exception of Experiment F, a punished connection 
gains more strength from its first occurrence than it loses by being punished. 
In my opinion, if the situation vanishes so promptly that the person is not 
led by the punishment to change his response to it, i.e., to create some 
second and different connection, and if he is immediately engaged in some 
new task so that he does not create some second and different connection 
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TABLE 6 
THE STRENGTHENING DUE TO ONE OCCURRENCE OF A PUNISHED CONNECTION 5 OR 
More STEPS REMOVED FROM ANY REWARD. EXPERIMENT J. IN PERMILLES. 


EXCESS OF 
ACTUAL OVER 
REPETITION DUE TO CHANCE AND ACTUAL EXPECTATION 
NUMBER FAVORITISM REPETITION, BY CHANCE 
SET A SET B SET C SET D ESTIMATE PUN. CONN. AND NUMBER 
SUBJECT n= 245 mn = 245 n= 245 mn = 245 FoR ALL 50R>REMOVED FAVORITISM 
Ah 128 120 117 192 140 249 109 
An 119 119 117 121 120 151 31 
Bi 104 116 143 113 120 170 50 
Br 126 163 157 185 160 203 43 
Br 2 151 130 131 133 135 241 106 
Bro 111 112 106 112 110 187 77 
Burg 162 143 166 168 160 428 268 
Burl 129 116 119 121 120 202 82 
Burk *169 135 116 169 147 329 182 
Ca 108 120 116 109 115 131 16 
Co 137 122 146 148 135 236 101 
F 139 134 139 122 135 229 94 
G 151 119 161 136 142 315 173 
H 109 113 134 114 118 150 32 
Le 129 112 124 108 118 161 43 
Lip 120 121 135 123 125 259 134 
Ma 118 146 145 138 137 272 135 
Me 203 121 141 121 147 315 168 
N 121 126 134 122 126 187 61 
O 385 281 122 151 235 372 137 


with some memory or other equivalent of the situation, the punishment 
has zero influence upon the punished connection. 

If this is so, the strengthening observed in these experiments is a measure 
of the strengthening due to a mere occurrence, provided proper allowance 
is made for (a) spread or scatter beyond four steps, (b) the possible per- 
sistent suggestion of a particular response (here a number, a part of a 
diagram or a position in a line) by the situation, and (c) a possible occa- 
sional confirming reaction due to the satisfyingness of making any choice 
and response. 

Observation and experiment can, with enough labor, learn what allow- 
ances, if any, are needed for (a) and (b). We have evidence already that 
there is no spread or scatter of the confirming influence of a reward beyond 
four steps in many of the seventeen relevant experiments in which there is 
certainly a large excess over chance, and almost certainly an excess over 
chance plus number favoritism. There is also evidence that some persons 
who are entirely unaware of any persistent suggestions of particular 
responses of the sort described in (b) show large excesses over chance plus 
number favoritism. Though doubtless a real cause, b does not seem likely 
to be an adequate one. 
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I have not found it possible, within such multiple-choice experiments 
as these, to distinguish between a confirming reaction due to the satisfying- 
ness of making any choice and response and a strengthening influence due 
to the mere use of a connection. 

The former is conceivable; it would permit the explanation of all 
strengthening of connections by one physiological function or mechanism; 
it would explain the phylogenetic development of learning by one variation, 
the appearance of a confirming reaction, rather than by two; it would avoid 
the difficulty of explaining why the mere use of some connections does 
strengthen them whereas the mere use of others does not. But there is 
no direct evidence for it in the experiments. The subjects feel no observa- 
ble satisfaction at saying amy number or marking any part of a diagram. 
Their satisfaction seems to them limited to the connections that are re- 
warded (or that, in later stages of an experiment, they know will be re- 
warded, even before they make them). And, in general, the evidence here 
seems to be against it, as in the experiments of Chapter IV of the Funda- 
mentals of Learning, and as in the formation of conditional reflexes by 
Cason ['22], by Hudgins [’33] and by others who have used techniques 
in which the subjects gain nothing by the formation of such connections. 

There is still another possibility allied to (c). The mere occurrence of a 
connection may strengthen it and yet do so indirectly and occasionally 
by arousing the confirming reaction, not directly and uniformly by the 
action of the connection itself. I have elsewhere suggested, as a working 
hypothesis to explain the ‘‘spread’’ or “‘scatter’’ phenomenon, that the con- 
firming reaction, being a physiological rather than a logical influence, 
occasionally misses the connection which the reward follows and belongs 
to, and influences by mistake a punished connection coming a little earlier 
or later. In a somewhat similar way the confirming reaction itself, which 
is in such experiments probably in a state of extreme excitability or sensi- 
tiveness or readiness to act, may occasionally be set off, not by its normally 
adequate stimulus, a state of affairs satisfying to the mental set or purpose 
or controlling powers, but by mistake by stimuli that normally would be 
inadequate. This hypothesis would have the merits noted above for (c) 
and would perhaps be consistent with the facts of other experiments. 
But discussion of all these possibilities may well be deferred until further 
experiments bearing upon them have been made. 

1 Some of these are reported in An Experimental Study of Rewards by Thorndike 
and others, 1933; others have been or will be published in various journals. 

2 Some of this strengthening may be due to repetition by the subjects of the numbers 
after the announcements of right. There was little time for them to do this, and there 
is no surety that those who did it profited by it, but it was not forbidden and was often 
attempted. For the purposes of this article, we do not care just what caused the 


43.6 per cent of repetitions of a rewarded occurrence. 
3 For 16 of the 24 subjects, every one of the 4800 responses was tabulated. For the 
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other 8 subjects enough were tabulated (800 for two, 1200 for three, 2400 for one and 
3800 for two) to satisfy all requirements. In all cases, the first two sets of 200 each 
and the last two sets were tabulated. In the case of the two for whom only 800 were 
tabulated, there was close agreement among the four sets of 200. 

‘ The complete detailed records are kept on file in the office of the Division of Psy- 
chology of the Institute of Educational Research at Teachers College for use by any 
one who wishes to make independent estimates. ‘ 

5 First occurrences only are used; and, of course, the count of repetitions is for the 
trial or round following the occurrence. 

6 Experiments J, L and M are described in An Experimental Study of Rewards. 

7™In a study as yet unpublished. 
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THE EFFECTS OF THE PITUITARY HORMONES ON THE 
MELANOPHORES OF FISHES 


By J. M. ODIORNE 
ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated May 27, 1933 


In contrast to the uniform effect which adrenalin has on the melano- 
phores of many animals, the results obtained by the use of extracts from 
the pituitary body have frequently been contradictory. In view of this 
fact it seems appropriate to publish an account of experiments made 
within the last three years which put the situation in a somewhat different 
light. Most of the experiments have been carried out on Fundulus hetero- 
clitus and Ameiurus nebulosus, but occasionally other fishes have been 
tested. I am pleased to acknowledge the kind suggestions made by Pro- 
fessor G. H. Parker, under whose direction the work has been done. 

Spaeth (1918) reported that the melanophores on isolated scales of 
Fundulus contracted when placed in a pituitary solution. This was con- 
firmed by Wyman (1924), who found that injection of the solution into the 
body cavity brought about a paling of the fish; similar results were ob- 
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tained by Gilson (1926). Nevertheless, all experiments which I have 
made have been negative; the melanophores of isolated scales failed to 
contract in post-pituitary solutions and fishes injected with such solutions 
showed no pigmentary response. 

“Pituitrin,’’ an aqueous solution containing the water-soluble principles 
of the posterior lobe of the pituitary gland, was the first extract used. 
If scales from a white-adapted Fundulus are placed in distilled water the 
melanophores remain contracted, while the melanophores of scales placed 
in Ringer’s solution expand within a few minutes. Consequently, if 
pituitrin is added to the Ringer’s solution one would expect the melano- 
phores to contract. Spaeth reported that they contracted completely 
and that the time required for this response was a function of the concen- 
tration of the pituitary solution. In spite of this, in my experiments such 
a result was not obtained even when equal volumes of pituitrin and Ringer’s 
solution were used. Yet if a small amount of adrenalin was added to such 
a solution nearly all of the melanophores contracted within five minutes. 
When pituitrin was added to distilled water in which scales had been 
placed the melanophores remained unaltered, showing that it had no 
expanding effect. The melanophores of scales transferred from distilled 
water to Ringer’s solution expanded, which shows that the melanophores 
were not seriously harmed by distilled water. 

The results of experiments on isolated scales were confirmed by intra- 
peritoneal and intramuscular injections of pituitrin. Dark fishes from a 
black dish injected with pituitrin (in doses of from 0.1 to 0.4 cc.) and re- 
turned to the black dish remained dark. If transferred to a white dish 
they rapidly became light. Conversely, light fishes which were injected 
with pituitrin remained light if returned to a white dish, but rapidly 
darkened when transferred to a black dish. In other words, they behaved 
like normal uninjected fishes, and when in the same dish with controls 
could not be distinguished from them by their color. 

That the pituitrin was capable of producing pigmentary responses in 
suitable animals was demonstrated by experiments conducted simul- 
taneously on Ameiurus. The responses of these fishes to backgrounds are 
slow as compared to those of Fundulus, but they become light greenish 
on a white background and extremely black on a black background. 
Light fishes injected with pituitrin became darker, the degree varying 
with the size of the dose (usually 0.2 cc. produces maximal darkening of a 
fish six to eight inches long; the same effect can be produced in smaller 
fishes by 0.1 cc.). ‘The response was slow, requiring an hour or more for its 
completion, but the fishes remained dark for hours after the injection. 
As a check on experiments with Fundulus, pituitrin from the same ampoule 
was injected into Ameiurus with the result described, so that the efficacy 
of the preparations used in these tests was certain. 
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In addition to pituitrin, preparations of its two principles are available. 
“Pitressin,’’ a solution containing the pressor principle (beta-hypo- 
phamine), formerly known as vasopressin, produced results in Ameiurus 
like those brought about by pituitrin, and failed to produce any effect on 
Fundulus. ‘“‘Pitocin,’’ formerly called ‘oxytocin,’ represents the oxytocic 
principle (alpha-hypophamine); it produced no effect on Ameiurus, and 
was not tried on Fundulus. 

Experiments on eels and black-backed goldfishes gave results like those 
obtained on Ameiurus. The pigmentary responses of eels are even slower 
than those of Ameiurus, though goldfishes react at about the same rate. 

An extract of the anterior lobe of the pituitary gland is also available, 
and the effects of this preparation, which is known as ‘‘antuitrin,’’ have also 
been investigated with interesting results. In Fundulus the injection of 
0.2 cc. of antuitrin produced contraction of the melanophores of dark 
specimens kept on a black background. The response was prompt but 
not all of the melanophores were concerned in it, for the fishes presented 
a spotty appearance unlike the even light color which results from the 
injection of adrenalin. The pale condition did not persist for more than 
two hours. 

Experiments using isolated scales furnished confirmation. When 
antuitrin was added to a small amount of Ringer’s solution in which 
scales had been placed the expanded melanophores contracted within a few 
minutes. Thus it gave results similar to those obtained by Spaeth with 
pituitrin. 

It seems likely that the dissimilarity between the results obtained by 
Spaeth and Wyman and those described here can be traced to differences 
in the preparations used. Spaeth (1918) employed ‘‘pituitrin’’ (Parke, 
Davis and Co.) and a pituitary extract furnished by H. K. Mulford & Co.; 
as the title of his paper indicates, he believed that the melanophores of 
Fundulus could be used for the biological standardization of pituitary 
extracts and other drugs. The samples which he used had been standard- 
ized by other methods, yet he wrote that “‘commercial samples of pituitary 
extract which were supposedly of the same physiological strength were 
found, by the use of the melanophore method, to vary widely.’’ If we 
assume that these early preparations contained not only the principles 
of the posterior lobe but substances from the anterior lobe as well, then 
Spaeth’s results can be explained and we can understand why supposedly 
standardized samples should ‘‘vary widely.’’ Instead of checking the 
strength of the extracts he was measuring their impurity with regard to 
material from the anterior lobe. Such an assumption is not unwarranted 
in view of the anatomical relationships of the various parts of the pituitary 
and the diffusion of substances from one part to another which has been 
amply demonstrated. Wyman (1924) in his experiments used pituitary 
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preparations which were supplied to him by Spaeth and which were from 
the same lot used by Spaeth for his 1918 paper, so that it is not surprising 
that their results were similar. 

The first report that pituitary extract causes expansion of the melano- 
phores of fishes was published by Abolin (1925). Following injections of 
“‘infundin”’ he observed expansion of these cells in Phoxinus, though the 
response was modified by the activity of other chromatophores. He ob- 
tained similar results using Nemachilus barbatula, Esox lucius, Carassius 
vulgaris and Leuciscus rutilis. Hewer (1926) found, on the contrary, that 
extracts of the posterior lobe of the cod pituitary produced contraction 
of the melanophores in Phoxinus. Infundin gave rather inconclusive re- 
sults. Giersberg (1930) also studied Phoxinus and reported that infundin 
produces the effect claimed by Abolin; in the perch no such positive 
effect was obtained. In Gobius and Pleuronectes pituitrin and “hypo- 
physin’’ cause contraction of the melanophores, according to Meyer (1931). 
Lundstrom and Bard (1932), working on the dogfish Mustelis canis, added 
the elasmobranchs to the list of animals which respond to pituitary sub- 
stances. The pale condition resulting from extirpation of the entire 
hypophysis or of the posterior (neuro-intermediate) lobe is replaced by the 
normal dark color following injections of suspensions of the dogfish pos- 
terior lobe, or of infundin, pituitrin and pitressin, but pitocin has no effect. 

From this we can at least conclude that different kinds of fishes respond 
in very different ways to pituitary preparations. Of the various fishes 
mentioned some do not respond to pituitrin at all, in the majority the 
melanophores expand and the fish becomes darker, in others (Gobius and 
Pleuronectes) the melanophores contract, and in the case of Phoxinus the 
effect is disputed. In Fundulus the result previously attributed to pituitrin 
seems to be due to some principle contained in antuitrin. Recently it has 
been estimated that twenty-two pituitary hormones have been reported 
by various investigators. In view of this it does not seem surprising that 
more uniform results have not been obtained. Very likely some fishes 
respond to one principle but not to another, and it is not certain that all 
fishes would respond in the same way to a single substance. 

Hogben (1924), in his monograph “The Pigmentary Effector System,” 
wrote as follows: ‘One other fact emerging from Spaeth’s work (1917) 
is of special interest, because it shows very clearly that the physiological 
properties of Fish and Amphibian melanophores are not to be regarded as 
identical. In the case of Amphibia, we have seen that pituitary extract 
exerts a specific and very characteristic action in promoting melanophore 
expansion. Spaeth, who has suggested a method of pharmacological 
assay based upon his observations, claims that pituitary extract induces 
contraction of the melanophores in fish scales. While there is reason to 
believe that the melanophores of Fishes and Amphibia respond in a different 
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manner to other stimuli, this appears to be a specially clear-cut instance of 
dissimilarity, since there is satisfactory evidence that the action of pituitary 
extract upon Amphibian melanophores is local.’’ Later he concludes: 
“Tf Spaeth’s observations on the pituitary reaction of Fish melanophores 
are correct, the mechanism of pigmentary response in Amphibia and Fishes 
must differ profoundly.” 

In view of the results obtained on various fishes since that time it seems 
clear that, although the melanophores of fishes and amphibians are perhaps 
not to be regarded as identical, they are not profoundly different. Indeed, 
the unanimity of response to pituitary extracts in the amphibians has 
been set aside by the work of Collins and Adolph (1926), who reported that 
the melanophores of the salamander Diemyctylus contract under the 
influence of pituitrin. Consequently the most that can be said is that not 
all fishes and amphibians respond to pituitary substances in the same way. 

Even though the difference between amphibians and fishes is not as great 
as was previously thought, it is nevertheless true that the amphibians as a 
group are more uniform than the fishes in this respect. In Fundulus the 
pigmentary responses have been shown to be to a great extent under the 
control of the nervous system. Though the nervous system also plays a 
part in the control of the color changes of Ameiurus and the eel, the slower 
responses of these forms suggest the possibility that here the nervous 
system is materially aided by hormones. Possibly the differences in the 
responses of fishes to pituitary substances are somehow related to the 
varying importance of the nervous system and the humoral system in 
these forms. 
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THE OCCURRENCE OF GUANOPHORES IN FUNDULUS 


By J. M. ODIORNE 
ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated May 27, 1933 


Though the changes in shade and color which occur in Fundulus hetero- 
clitus in response to backgrounds and other stimuli have hitherto been 
accounted for by changes in the melanophores and xanthophores, this fish 
also possesses chromatophores which apparently have been neglected if not 
entirely overlooked. Perhaps this is because of the distribution of these 
cells and their position in relation to the melanophores. The object of 
this report is to describe briefly some of their morphological and physio- 
logical characteristics, and to point out their significance in bringing about 
alterations in the appearance of the fish. 

If one of the large scales in the mid-dorsal region just behind the head is 
removed from an individual adapted to a black background, and the area 
thus exposed is observed with a microscope by reflected light, it will be 
seen to present the features illustrated by figure 1. The most conspicuous 
structures are the expanded melanophores; their processes partly conceal 
the xanthophores which lie among them, but deeper in the dermis. Im- 
bedded in many of the melanophores are the iridosomes which have been 


the subject of a recent paper by Foster (1933). According to the termi- 
nology of Ballowitz (1912) an iridosome is a group of “‘iridocytes,’”’ and 
when associated with a melanophore the combination becomes a “mel- 
aniridosome.”’ 


If instead of a black-adapted specimen a fish from a white background 
is examined in the same way, a much different condition is revealed 
(Fig. 2). The melanophores are contracted and appear to be surrounded by 
silvery halos, which are actually the radiating processes of chromatophores 
lying beneath the melanophores. These processes form an interlacing 
network in many places and hide the deeper xanthophores from view. 

If a piece of integument is removed and studied by transmitted light, 
the chromatophores, which always appear silvery (or luminously white) 
by reflected light, lose their brilliance and instead are a dull yellowish 
gray. They are never iridescent, though their brilliant whiteness is due 
to their ability to reflect incident light. If examined under oil immersion 
the substance which they contain is seen to be in the form of small particles. 
These may appear to be closely packed together or they may show Brown- 
ian movement. In either case it is difficult to determine whether they are 
merely granular or definitely crystalline, for they are smaller than the 
melanin granules. Some writers speak of amorphous guanin, while others 
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hold the opinion that it is always crystalline but, whatever its state, it 
seems likely that the substance is guanin. 

Consequently I shall refer to them as “guanophores,” a term used by 
Schmidt (1912) to indicate chromatophores (in amphibians) containing 
guanin. Ballowitz (1912) applied this term to similar chromatophores 
of fishes which had been called “‘iridocytes’’ by Pouchet (1876), but 
Ballowitz continued to use 
Pouchet’s term in his num- 
erous papers. Both Con- 
nolly (1925) and Fries (1931) 
have referred to the irido- 
somes of Fundulus as 
“guanophores or __ irido- 
cytes,” but in spite of the 
fact that an iridocyte is a 
type of guanophore it seems 
desirable to use the word 
“Sridocyte’” only when re- 
ferring to a guanophore 
which is iridescent. Bied- FIGURE 1 


ermann (1926) has already Photomicrograph of the area exposed by the re- 
pointed out that there are moval of a scale in the mid-dorsal line of a black- 
adapted Fundulus. The white spots in many of 
the melanophores are the iridosomes. The xantho- 
phores appear as small, dark gray bodies. 





guanophores which show no 
interference colors what- 
ever, and the guanophores 
of Fundulus which are the subject of this paper are of this type. 
Gilson (1926), in describing the iridocytes of Fundulus, said: ‘By 
reflected light the iridocytes appear shining and silvery, but by 
transmitted light they are brownish in appearance and, in a hasty exami- 
nation of tissue, may even be mistaken for melanophores.”” He did not 
describe their location, neither did he give illustrations, but from his 
description it would seem that he observed these structures though he 
did not recognize them as different from the true iridocytes. 

Ginsburg (1929) has described the development of the chromatophores 
in the guppy, Lebistes reticulatus, and he has recognized two types of 
“guaninophores,”’ as he called the guanin-containing cells. In addition 
to those which show interference colors there are others which are generally 
much larger, and which appear dull white or silvery white by incident 
light. Recently Sumner and Wells (1933) have described and illustrated 
structures in Lebistes which they have interpreted as ‘‘iridosomes,”’ 
following the terminology of Ballowitz. Apparently, they are the silvery 
white ‘‘guaninophores’? mentioned by Ginsburg. The guanophores of 
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Fundulus bear a very close resemblance to these “‘iridosomes”’ or “‘guanino- 
phores”’ of Lebistes. 

Ballowitz, in a series of contributions, has described chromatophores 
and chromatophore combinations in a number of different fishes, and an 
examination of these papers makes it clear that iridosomes may not only 
be dissimilar in different genera, but also that a single fish may have more 
than one type, as in Serranus. The melaniridosomes of Fundulus most 
closely resemble those of Perca and Acerina (Ballowitz, 1914, Taf. I. Figs. 
1-5). The guanophores of Fundulus are somewhat similar to the irido- 
somes of Gobius (Ballowitz, 1913, Taf. XI, Figs. 72, 73 and Taf. XII, 
Figs. 94a, 94b); to be sure the optical properties are not the same, but the 
configuration and position 
are comparable and both 
expand and contract.! Bal- 
lowitz did not investigate 
the effects of different back- 
grounds on the iridosomes 
of Gobius, but Meyer (1931) 
reported that they expanded 
in response to a white back- 
ground and contracted un- 
der the influence of a black 
one. 

The guanophores of Fun- 
dulus are most numerous 
on the dorsal surface though 
they extend downward over 





FIGURE 2 


Photomicrograph of an area similar to that 
shown in figure 1, after the application of adrenalin. 
The melanophores are contracted and the guano- _ the sides for some distance. 
phores are widely expanded, as in a fish adapted Infemales they are absent or 
to a white background. rare in the fins and tail but 

in males small guanophores, 
usually without overlying melanophores, are found in the dorsal, anal 
and caudal fins. They occur in areas which are free from melanophores 
(Fig. 3), frequently accompanied by iridescent patches, and intensify 
the color pattern which is primarily caused by the absence of the melano- 
phores. Isolated guanophores are occasionally found on the body as well. 
In either situation they are small and when expanded do not exhibit 
such extensive processes as do those below melanophores. In some fishes 
guanophores also occur below the melanophores which lie in the dermis 
external to the scales, and near their edges. If a scale bearing such guano- 
phores is removed and examined from the under side, they will be found 
to be small. They are located near the margin of the pigment mass (in a 
contracted melanophore) and usually do not cover its entire under surface. 
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The guanophores seem to be most numerous where the stratum argen- 
teum is poorly developed, though exception must be made of the fins, 
especially those of the females. Very likely the reason that they have 
escaped notice is that they are rare or absent in those parts of the body 
which lend themselves to the study of melanophores and xanthophores, 
namely, the isolated scales and the caudal fin. Besides they are not 
conspicuous when viewed by transmitted light. 

That the guanophores of Fundulus do have the ability to expand and 
contract is shown by the following experiment. If a few drops of adrenalin 
(1:1000) are put on an area such as that illustrated in figure 1 the migration 
of the melanin granules toward the centers of the melanophores can be 
plainly seen within a short time. Presently the outward migration of 
the guanin can be discerned 
and within about five min- 
utes the appearance of the 
area is similar to that shown 
in figure 2. From this it is 
clear that these structures 
are chromatophores and are 
not merely immobile de- 
posits of guanin beneath 
the melanophores which are 
exposed and concealed by 
the movements of the pig- 
ment within these cells. 
The isolated guanophores 
are favorable for observing 


these movements of the Photomicrograph of a portion of the dorsal fin 
of a male Fundulus, showing isolated guanophores. 





FIGURE 3 


guanin, because of the ab- 
sence of overlying melano- 
phores, but as already pointed out the processes are rather small. 

The experiment just described also demonstrates that the guanophores 
respond to adrenalin, but the response is opposite to that of the melano- 
phores. The same reaction occurs if 0.2 cc. of adrenalin (1: 1000) is in- 
jected into the body of a dark fish, instead of being applied to the naked 
dermis. 

‘“Antuitrin,” an extract of the anterior lobe of the pituitary body, causes 
contraction of the melanophores of Fundulus when injected into the body 
cavity or when applied to the naked dermis (see paper preceding this). 
The response is not as thorough as the response to adrenalin but is per- 
fectly definite. The guanophores expand under these conditions. 

In a fish from a yellow background the melanophores are contracted 
though not always maximally; at the same time the guanophores are 
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partially expanded but are only slightly visible. Consequently the xantho- 
phores are not concealed by either melanophores or guanophores and the 
animal assumes a decidedly yellowish green color. 

For the detailed study of the guanophores it is necessary to remove 
the overlying scales but this does not mean that their presence is not 
evident without this procedure. Different individuals vary with respect 
to the completeness with which they adapt their color to a white back- 
ground, and with respect to the time required to become maximally light, 
but the guanophores are certainly important factors in this adaptation, 
even though not all fishes assume the decided grayish color frequently 
seen. Situated as they are below the melanophores, the guanophores 
do not detract from the dark appearance of fishes in which the melano- 
phores are expanded; in light-adapted fishes the greatly expanded guano- 
phores increase the paleness caused by contraction of the melanophores. 
Furthermore, the reactions of the guanophores to various stimuli have 
so far been found to be just opposite to those of the melanophores, a fact 
which shows clearly that the chromatophores of this fish are far from being 
physiologically alike and that their behavior is decidedly ‘‘adaptive.”’ 


1 It has become customary to refer to the changes which take place in chromatophores 
as “‘expansion’’ and “‘contraction,”’ though this does not imply that the shape is actually 
altered. Instead, the pigment within the cell shifts, now being massed in the center 
(contraction), now distributed in the processes (expansion). 
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Ortmann (1911) describes a peculiar rhythmical movement of the mantle 
flaps of females of two members of the genus Lampsilis (L. ventricosa 
and L. multriradiata), fresh-water lamellibranchs. This movement was 
noted by the author in a female of Lampsilis nasuta. With the knowledge 
that the siphons, mantle edge and labial palps of bivalve molluscs often 
contain photoreceptors, it occurred to the writer to determine whether or 
not this rhythmical movement was influenced by light intensity. Such 
movement was being sought for in order to observe the effect of light 
intensity on frequency of muscular contractions where amplitude was not 
involved as it is in leg movement (Welsh, 1932). A few preliminary trials 
indicated that light did play an important rédle in determining the fre- 
quency of these rhythmical contractions. A series of determinations of 
the time required for ten movements, over a fairly wide range of intensities, 
is contained in table 1. The frequency of movement is shown plotted 


TABLE 1 


TIMES IN SECONDS FOR TEN MOVEMENTS OF THE MANTLE FLAPS OF LAMPSILIS WITH 
THEIR AVERAGES, AND THE FREQUENCY (NUMBER OF MOVEMENTS PER SECOND) AT 
EACH OF SEVERAL INTENSITIES OF ILLUMINATION. TEMP. 21.3°C. 


INTENSITY 
(FOOT-CANDLES) 0.20 0.41 0.83 1.66 3.6 12.0 23.0 180.0 689.0 


Time (secs.) for 
10 movements 65.0 48.0 28.6 24.5 22.5 21.2 21.0 20.1 20.0 
61.0:. 46:8 30.6268 22:3 2.8 -20:8. 2:1. G0 


62.0 49.2 32.2 27.4 22.2 21.4 20.8 20.1 20.2 

63.5 48.6 35.2 26.6 22.1 21.5 20.6 20.1 20.4 

61.8 46.7 83.0 26.4 - 21:9. 21:8 2.0: 24> 20:1 

32.5 27.6 22.6 21.5 20.9 20.2 19.7 

30.8 27.6 23.4 21.2 21.0 20.4 19.6 

35.0 27.4 23.6 21.5 20.6 19.9 19.8 

36.0 27.0 23.3 21.0 20.6 20.4 20.0 

38.0: °26.6% 23.6 + 21.:2>-. 21.2 -2:4» 1938 
Averages 62.66 47.86 33.18 26.59 22.75 21.32 20.84 20.19 19.96 
Frequency 0.159 0.210 0.301 0.376 0.489 0.469 0.481 0.495 0.502 


against the logarithm of the light intensity in figure 1. The curve is similar 
to that obtained in a study of frequency of leg movement, as a function 
of light intensity, in a water mite Unionicola (Welsh, 1932), except that the 
lower portion is incomplete as the rhythm was interrupted at low light 
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intensities, and ceased entirely after a short exposure to an illumination 
of about 0.2 foot-candles. Five consecutive readings at the two lower 
intensities were as many as could be obtained without interruption. 

It was impossible to carry out check experiments on other individuals 
because of their rareness locally, and the fact that the contractions occur 
only for a period of a few weeks when the marsupia are loaded with de- 
veloped glochidia and these are being discharged. 
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FIGURE 1 


Plot of data of table 1 showing frequency of movement (number of move- 
ments per second) of the mantle flaps of Lampsilis nasuia, plotted against the 
logarithm of the light intensity. Observations were begun at the highest 
illumination. 


Ortmann (1911) believed this movement to be for purposes of aeration 
of the glochidia during the breeding season. This seems to be the most 
logical explanation as the action of the mantle flaps greatly increases the 
amount of water passing through the mantle cavity; however, if this is 
its only function it is difficult to understand why it stops at low intensities 
of illumination or in darkness. In nature the movements must occur only 
during daylight hours, and during the night the circulation must be the 
result solely of ciliary activity. The phenomenon deserves further in- 
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vestigation as other factors such as CO, tension would doubtless be found 
to influence the frequency of movement as Redfield (1917) found in a 
study of rhythmic contractions in the mantle of certain lamellibranchs, 
a quite different process but one having a similar end result. 

The very precise regularity of the movement nearly equals the regularity 
of heart beat and provides a very convenient rhythm for use in an analysis 
of a number of phases of the physiological environment of such fresh-water 
organisms. It is proposed to investigate further this rhythmical move- 
ment when more breeding females are obtained. 
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